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Measurements of fluctuating surface pressure were made on a NACA 0015 airfoil 
immersed in grid generated turbulence. The airfoil model has a 2' chord and spans the 6' 
Virginia Tech Stability Wind Tunnel test section. Two grids were used to investigate the effects 
of turbulence length scale on the surface pressure response. A large grid which produced 
turbulence with an integral scale 13% of the chord and a smaller grid which produced turbulence 
with an integral scale 1.3% of the chord. Measurements were performed at angles of attack, a 
from 0° to 20°. An array of microphones mounted subsurface was used to measure the unsteady 
surface pressure. The goal of this measurement was to characterize the effects of angle of attack 
on the inviscid response. 

Lift spectra calculated from pressure measurements at each angle of attack revealed two 


distinct interaction regions; for co r 



< 10 a reduction in unsteady lift of up to 7 decibels 


(dB) occurs while an increase occurs for 0 ), > 10 as the angle of attack is increased. The 
reduction in unsteady lift at low co r with increasing angle of attack is a result that has never 
before been shown either experimentally or theoretically. The source of the reduction in lift 
spectral level appears to be closely related to the distortion of inflow turbulence based on 
analysis of surface pressure spanwise correlation length scales. Furthermore, while the distortion 
of the inflow appears to be critical in this experiment, this effect does not seem to be significant 
in larger integral scale (relative to the chord) flows based on the previous experimental work of 
McKeough (1976) suggesting the airfoils size relative to the inflow integral scale is critical in 
defining how the airfoil will respond under variation of angle of attack. 

A prediction scheme is developed that correctly accounts for the effects of distortion 
when the inflow integral scale is small relative to the airfoil chord. This scheme utilizes Rapid 
Distortion Theory to account for the distortion of the inflow with the distortion field modeled 
using a circular cylinder. 
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Chapter One 


1. Introduction 

When a body moves through an unsteady fluid the interaction of the two produces 
unsteady pressure fluctuations on the surface of the body which may radiate to the far-field (i.e. 
generates noise). The study of fluid/body interaction is of particular interest due the resulting 
noise generation (aerodynamic noise) as well as the unsteady loading produced which may 
induce undesirable vibration. 

A rather complex flow field develops as a body moves through an unsteady fluid. Near 
the surface of the body a boundary layer develops transitioning from laminar to turbulent. 
Associated with the turbulent boundary layer are tiny eddies which vary spatially and temporally 
and thus produce pressure fluctuations on the surface of the body. The airfoil also experiences 
fluctuating pressure as a result of the unsteadiness in the incoming fluid. These pressure 
fluctuations are associated with the inviscid response of the airfoil, or in other words, the 
imposition of the non-penetration condition. Both types of pressure fluctuations, those arising 
from the boundary layer and those from the inviscid response, are sources of noise and vibration. 
If the body is aerodynamic (such as an airfoil) and the characteristic scale of the unsteady motion 
(i.e. turbulence) on the order of the body dimensions the dominant source of noise arises from 
the inviscid response. 

The noise generated from turbo-machinery, helicopters, aircraft, and marine vehicles has 
made the study of fluid/body interaction quite important. However, understanding this complex 
interaction mathematically or experimentally has proven to be quite difficult. 

1.1 Background 

1.1.1 Pressure Fluctuations 

The study of aerodynamic sound amounts to understanding how rotational kinetic energy 
is converted to acoustic waves involving longitudinal vibration] . Lighthill (1952, 1954) 
developed the theory of aerodynamic sound through the acoustic analogy. Lighthill (1952, 1954) 
considers a region of fluctuating or turbulent fluid which occupies a larger volume of fluid at 
rest. He proposes that the equations describing the density variation in this fluid are equivalent 
to equations describing a uniform acoustic medium subject to an external fluctuating forcing 


1 Howe, M. S., Acoustics of Fluid Structure Interactions, Cambridge Press, 1998 
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function. Using this analysis Lighthill (1952, 1954) derives an inhomogeneous, exact wave 
equation for the propagation of sound generated aerodynamically and is given as, 


^-c;VV= dT * 

dt 0 


dxpXj 


where, Ty is the Lighthill stress tensor defined as, 


T ij =pu i u j +T ij +(p-c 2 0 p)S ij . 


( 1 . 1 ) 


( 1 . 2 ) 


The term pu j u j is the Reynolds stress tensor, r (/ is the viscous stress tensor, and p-c 2 0 p is the 


thermal stress. The Lighthill stress tensor represents the external stress system acting on the 
uniform acoustic medium and is the difference in the effective stresses in the flow and the 
stresses in the uniform acoustic medium at rest. 

The most general solution of an exact, inhomogeneous wave equation was first given by 
Stratton (1941) in his book Electromagnetic Theory and makes concessions for a control surface, 
S. If the surface is stationary and rigid, as in the case of an airfoil immersed in turbulent flow, 
the sound generated can be determined by implementing Curie’s (1955) solution to Lighthill’ s 
equation Here the time domain Greens function is used to derive the outgoing wave solution in 
integral form. In doing so, the notion of retarded time is developed to mean the time for a 
disturbance at location y to propagate to location x external to a boundary defined by surface 5. 
Curie’s equation is given as 


c 


2 

o 


(P-Po) 




(1.3) 


where quantities in [] brackets are evaluated at the retarded time r R = t - ^ % and 
P\ = {p-Po)* v -0 9 is the compressive stress tensor. Also, 8 is the Kroneker delta-function. 


Rearranging equation (1.3) allows the pressure in the far-field to be approximated as 

*>(*>*)” X %? h: J [/w)*ra+ Xi h ,2 jA[n jP '^ds (1.4) 

4^-c 0 \x\ ot J ol 4 nc 0 |x| ot J s 

where n is the unit normal on s directed into the fluid. In doing so, the aerodynamic sound is 
represented with sound produced by quadrupole (first term, volume integration) and dipole 
(second term, surface integration) sources. Dipole sources result from the net unsteady force 
exerted on the fluid by the airfoil or the inviscid response. Quadrupole sources result from the 
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turbulence stress distribution and are typically associated with the turbulent boundary layer. 
Under certain circumstances equation (1.4) can be simplified to include only the dipole term. 
This is evident through comparison of the order of sound power of each source. The dipole term 
generates sound power on the order ~A 1 p 0 u i M i , where A is the characteristic length of the 

turbulence. A quadrupole term produces sound power on the order ~A 1 p 0 u i M 5 . Evaluating the 

ratio of sound power of each, a factor of 1/M 2 persists. Therefore, the acoustic power in the 
dipole term exceeds that from the quadrupole by a factor of 1/M 2 . The radiation of sound is thus 
dominated by the dipole term when the Mach number is small {\JJc 0 « 1). This interaction 
situation is defined as acoustically compact and as such can be represented solely by dipole- 
sources. Equation (1.4) then becomes 

p{x,t)* - *' . 2 ^-j>[n jP ' v ]ds (1.5) 

4;rc 0 1*| at J s 

If the shear stress components of p'y are small compared to the local surface pressure fluctuation 
p\ equation (1.5) can be further simplified to give a simple relation between the acoustic and 
surface pressure 

p(x,t)& cos ^ —<b\p'\ds (1.6) 

v ’ 4 7rcxdt TL J 

o s 

where #is the angle between the outward normal and a vector which originates from the 
location of the outward normal and extents to the location where p is to be determined. 

1.1.2 Unsteady Thin Airfoil Theory 

Equation (1.6) clearly demonstrates if one wishes to obtain the acoustic pressure (far- 
field noise) an accurate estimate of the surface pressure is required. Thus, researchers embarked 
on a mission to develop techniques for computing the unsteady surface pressure of an airfoil 
interacting with turbulence. However, this is a forbidding task determining the unsteady surface 
pressure of an airfoil with thickness, at angle of attack, and immersed in a complex flow. 
Consequently, researchers worked to simplify the problem to a tractable nature. In doing so they 
developed unsteady thin airfoil theory. As the name suggests the airfoil is assumed to be thin, in 
fact having no actual thickness. Also, the airfoil is not at angle of attack and viscous effects are 
neglected (an implicit assumption in moving from equation 1.5 to equation 1.6). The airfoil can 
then be modeled as a vortex sheet with strength y (x,t) subject to a sinusoidal inflow with a wake 
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extending to infinity, y(x + tUf) (as shown in figure 1.1). The fluid motion is irrotational 
everywhere else with a perturbation velocity defined by a velocity potential which vanishes at 
infinity. Other conditions imposed include the Kutta condition on the trailing edge, a kinematic 
boundary condition on the airfoil surface and a dynamic boundary condition on the vortex wake. 
Lastly, since the airfoil was once at rest the total circulation about the airfoil and its wake must 
vanish. 

Consider figure 1.1; 




X 


Wake extending to infinity 


Figure 1.1: Unsteady Thin Airfoil Theory 

the total flow (a sinusoidal gust with wave front parallel to the airfoil leading edge) is written as 

U(r ,t) = U a0 +q(r ,t) 

( 1 . / 

= { U + U a0 )i + vj 


where q{r,t ) = ui +vj is the perturbation velocity and r =xi + yj is the position vector. It is 
convenient to write the v velocity component (or upwash) as w 0 sm(k x x + cot) or similarly as 
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Re|^H’ 0 e ,(i>A+ ® n J where w 0 is the upwash amplitude; in doing so any form of homogeneous 

turbulent upwash (in the 2-dimensional sense) can be built using Fourier analysis. Note that the 
frequency, co is related to the velocity as co = k x U v . Using this form of the upwash the 

perturbation velocity becomes 

q = ui +w Q sin (k x x + cot}j (1.8) 


Now, the flat plate with circulation is represented by a distribution of vorticity with strength y. 
As such, potential flow theory gives the following relationship between the perturbation velocity 
and vorticity, in complex number form 

oo 

u(z,t)-iv(z,t) = £ | ] ^jdz o (1.9) 

-c/2 

where / = V-l, z = x +iy and z a = x 0 +iy 0 . Applying the dynamic boundary condition on the wake 
gives 

y(x,t) = r(z 0 +Uj) c/2 < x < co (1.10) 

and decomposing equation (1.9) into its real and imaginary parts produces the velocity 
components in terms of vorticity 

oo 

u(r,t)= £ ( U1 ) 


c/2 oo 

-c/2 c/2 


r(Z+u*>Wx-Z) 

(x-if+y 2 


cU; 


( 1 . 12 ) 


The complete solution can now be constructed by imposing the kinematic boundary condition 
(no flow through the surface) on the vertical velocity (1.12) which gives 

c/2 oo 


-w a sin (k x x + cot) 


-i 

2tt 


I 

-c/2 


r(g,y,tXx-g) 

(x-g) 2 +y 2 



r(Z+uj)(x-Z) 

(x-t) 2 +y 2 




(1.13) 


Equation (1.13) governs the relationship between the inflow and the vorticity under the 
assumption prescribed above. Through inversion of (1.13) the vorticity can be obtained and 
utilized to compute the unsteady pressure on the airfoil surface. With the surface pressure 
known equation (1.6) can be used to determine the pressure in the far- field. 
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The relationship between the pressure and vorticity is now presented. First, defining the 
pressure difference between the upper and lower side of the airfoil 


A P = P V -P L (1-14) 

where Pu is the pressure on the upper side of the airfoil and Pl the lower side. By unsteady 
Bernoulli’s equation the velocity around the airfoil can be related to the surface pressure as 


AP = Pu-P l = 


^p {u ^ Uu f +p ^, 

dt 2 V UJ dt 


' L +£- 


(U M + u L f 


(1.15) 


with faj and fa are the velocity potential on the upper and lower surface and iiu and ul the 

perturbation velocity in the i direction on the upper and lower surface. Upon linearization of 
equation (1.15) the following relationship between unsteady pressure and velocity is obtained 

AP = + pU„(u L - Uu ) (1.16) 

This equation can now be cast in terms of the vorticity by noting that equation (1.11) can be 
simplified when evaluated on the upper and lower surface of the flat plate to give 

u(x,y = 0 ± ,t)=±^y(x,t) (1.17) 

where y = 0* indicates the upper and lower surface of the flat plate. Equation (1.16) can then be 
recast as 


A P = -pU aD y{x, t)-j- J y(x, t)dx (1.18) 

—oo 

A relationship has now been established (equation 1.13 with 1.18) between a sinusoidal 
inflow and the surface pressure. The crucial step in this solution is the ability to solve for the 
vorticity in (1.13). Considerable effort has been spent on extracting a solution to (1.13) in 
addition to extending Unsteady Thin Airfoil theory to the more realistic 3-dimensional case. The 
focus is now shifted to solutions to (1.13) and the search for rational representations of the airfoil 
response. 


1.2 Theoretical Unsteady Response of an Airfoil 

1.2.1 The Exact Solution to Unsteady Thin Airfoil Theory 

Von Karman and Sears (1938) followed by Sears (1941) were some of the first to pioneer 

the theory of turbulence/airfoil interaction. Using the Unsteady Thin Airfoil approach Sears 
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derived an explicit expression for lift and moment that involve a combination of Bessel functions 
which have become known as the Sears function. 

Liepmann (1952) extends Sears theory to consider the case of airfoil encountering 
turbulence of scale much greater than the span and utilizes statistical concepts in development of 
the lift response. He formulates the mean Cl squared in terms of an admittance function (Sears 
functions) and the power spectrum of a sinusoidal gust. This method is extended by Liepmann 
(1955) to allow the usage of the frequency spectrum of the incoming turbulence. Specifically, he 
introduces a two-wavenumber spectrum (spanwise and streamwise) to describe the turbulence 
which allows consideration of the more practical situation wherein turbulence scale is smaller 
than the span. The sinusoidal variation in angle of attack along the span is accounted for by 
taking the section lift at each spanwise location to be given by the local angle of attack and the 
two-dimensional Sears function. Ribner (1956) furthers this statistical approach to allow 
consideration of the full, 3 -wavenumber turbulence spectrum and three dimensional response 
function. 

1.2.2 Approximate Analytic and Numeric Models 

Unsteady Thin Airfoil Theory is a 2-dimensional theory and therefore does not easily 

allow for the more realistic case of a sinusoidal gust with wave front skewed relative to the 
airfoil leading edge. Additionally, Unsteady Thin Airfoil Theory does not account for 
compressibility which is quite important in the production of aerodynamic sound. Therefore 
several approximate closed form and numerical solutions have been developed to address these 
cases. 


1.2. 2.1 Skewed Incompressible Gust 

The skewed gust case was first considered by Filotas (1969) who found a closed form 
expression for the problem of an infinite, thin airfoil encountering a skewed, incompressible gust 
at zero angle of attack and derived approximate expressions for the lift and pressure distribution. 
His solution is based on Unsteady Thin Airfoil theory and utilizes an interpolation scheme to 
solve for the vorticity. His results were later demonstrated to be rather crude. Graham (1970a) 
also considered the skewed gust case and formulated a solution for the vorticity based on a 
Chebyshev series. The analysis is rather cumbersome in nature based on the consideration of 
both streamwise and spanwise shedding of vorticity due to the varying angle of attack along the 
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span. Graham provides the equations to numerically compute unsteady lift, chordwise pressure, 
and moment. 

Mugridge (1970) provides a closed form expression for the lift response of a thin airfoil 
of finite span immersed in a convected skewed sinusoidal upwash. He also showed the 
interesting effect that as spanwise wavelength of the gust becomes small compared to the span 
cancellation of surface pressure fluctuations occurs everywhere except at the extreme edges of 
the airfoil span. 

Amiet (1976a) developed a skewed gust incompressible solution valid at high reduced 
frequencies and large spanwise wave numbers from which the full compressible, skewed gust 
response can be derived. This work is extended by Amiet (1976b) to small spanwise 
wavenumbers utilizing Graham’s (1970b) similarity rules with Amiet’s (1974) compressible gust 
solution thus providing a solution to the complete spanwise wavenumber range. 

1.2. 2. 2 Compressible Parallel Gust 

The interaction of turbulence with a lifting body and the noise produced are tied closely 
to the compressible nature of the free stream fluid. It is therefore critical to consider the effects 
of compressibility in the turbulence/airfoil interaction problem. Early consideration of this 
important effect has been given by Possio (1938) and Reissner (1951). Reissner (1951) builds on 
Possio’s (1938) work deriving a solution in terms of the velocity potential using Mathieu 
functions. From this solution he formulates expressions for unsteady lift and moment. Amiet 
(1974) modifies Miles (1950) solution to consider the case of a compressible two-dimensional 
gust problem and develops a closed form expression valid at low frequencies. Amiet (1976a) 
later derived a closed form expression for the two-dimensional, compressible upwash problem 
valid at high reduced frequency. He employed the first two steps of an iterative procedure given 
by Landahl in which the leading and trailing edge boundary conditions are alternately applied. 
For this reason the theory only applies to high frequencies when the acoustic wavelengths are not 
large compared to the chord. 

Amiet (1975) also developed a relationship between the lift response function and 
turbulence spectrum which allowed prediction of far field acoustic radiation. In addition, he 
showed acoustic predictions can be made under certain conditions using only the unskewed gust 
response and limited information about the spanwise wavenumber content of the turbulence. 
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1.2. 2. 3 Skewed Compressible Gust 

Having developed solutions for the parallel, compressible and skewed, incompressible 
gusts researchers began to formulate approximate solutions to the more realistic case of a 
skewed, compressible gust. Graham (1970b) extends his work (although still computational in 
nature) to the case of a skewed, compressible gust encountering a flat plate at zero incidences. 
He transforms this problem into two sub-problems via a Prandtl-Glauert transformation 
depending on the spanwise trace speed of the gust along the airfoil leading edge. In the case of a 
subsonic trace speed the plate response is expressed in terms of the incompressible skewed gust 
problem which has been solved by Graham (1970a). For a supersonic trace speed the response is 
recast in terms of the unskewed, compressible gust which has been solved by Possio (1938). 

Adamczyk and Brand (1972) formulated a closed form solution for the skewed 
compressible gust which allows calculation of far field noise. Adamcyzk (1974) wishing to 
more clearly define the unsteady pressure, lift and moment response developed an approximate 
theory for the compressible, skewed gust which compared favorably to the exact theory of 
Graham (1970b). 

Chu and Widnall (1974) extend Amiet and Sears’ (1970) approximation to derive a 
similarity rule for the three-dimensional interaction of a two-dimensional wing and convected 
oblique sinusoidal gust in compressible flow valid at low reduce frequency. Osborne (1971) 
employed the GASP approximation to transform the Sears’ two-dimensional, incompressible 
solution into one for compressible flow in a straightforward analytic form. Although, as pointed 
out by Amiet (1974) this method cannot be utilized for the case of a two-dimensional airfoil 
since the solution involves integration over the infinite span. Therefore, Sing and Widnall 
develop a similarity rule which is valid in the presence of an infinite span wing. This method 
coupled with either Sears’ (1941) two-dimensional or Filotas’ (1969) three-dimensional 
incompressible flow solutions can be used as the bases for construction of analytic expressions 
for unsteady lift. 

A second order theory for the three-dimensional, compressible case was developed by 
Graham and Kullar (1977). These researchers used a general perturbation series to derive a 
second order solution for lift coefficient and compared their results with those of Amiet (1974) 
and other approximate solutions. They demonstrated insufficiencies in Amiet’ s solution as 
spanwise wavenumber tends to zero. 
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Martinez and Widnall (1980) unified aerodynamic/acoustic theory through the derivation 
of a closed form expression which couples the airfoil response with the acoustic field. Building 
on the work of Amiet (1976a) and Adamczyk (1974) they give a high frequency approximation 
to the airfoil response for a skewed, compressible gust from which both surface pressure 
fluctuations and far field noise can be determined. Their solution agrees with Adamczyk’ s 
(1974) for the case of zero sweep and collapses to Amiet (1976a) for the case of a parallel gust. 

1.3 Effects of Angle of Attack and Thickness on the Airfoil Response 

Unsteady Thin Airfoil theory provides well as a tool for understanding critical aspects of 

turbulence-airfoil interaction. However, real airfoils have thickness and are typically at non-zero 
angles of attack. If one is to produce superior aerodynamic devices (in terms of reducing 
aerodynamic noise) an improved representation of the airfoil and interaction is mandatory. 
Therefore, researchers aware of the deficiencies of Unsteady Thin Airfoil Theory began to 
investigate possible avenues for the incorporation of thickness and mean loading (i.e. angle of 
attack) effects. 

Perhaps the earliest attempt to consider angle of attack was that of Horlock (1968) and 
later extended by Morfey (1970). These researchers present a method for determining the lift 
fluctuation of an airfoil at angle of attack due to a gust parallel to the uniform flow field. 
Combining this method with established theory for transverse gusts (perpendicular to 
undisturbed flow) Horlock arrives at the fluctuating lift produced by both gust types. In his 
formulation Horlock argues that unsteady lift fluctuations arise not just by upwash fluctuations 
(producing an instantaneous change in angle of attack) but also by the streamwise velocity 
fluctuations changing the scaling of the lift coefficient. 

Atassi (1984) develops a theory for computing the unsteady lift of a thin airfoil with 
small camber and at small angle of attack. The unsteady lift is constructed by linear 
superposition to the Sears lift accounting separately for airfoil thickness, camber, and angle of 
attack. By developing a solution built on the Sears function Atassi’s (1984) theory is constricted 
to the case of a parallel gust and thus does not well represent actual fluid/body interaction. 
However, being a closed form solution does make this theory desirable and convenient for 
comparison with experimental data. 
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1.3.1 Effects of Distortion: Rapid Distortion Theory 

Goldstein and Atassi (1976) were conceivably the first to consider the distortion of 

incoming turbulence in development of the airfoil response function. They employed the 
technique of Rapid Distortion Theory (RDT) which considers the distortion of an incoming 
vortical gust (tilting and stretching of vorticity vector) by the steady-state potential flow around 
the airfoil. They formulate a solution for the case of a two-dimensional, incompressible gust 
providing separate methods for calculating the effects of angle of attack, camber, and thickness 
on the unsteady lift response in terms of a correction to the Sears formula. Interestingly, they 
show the effect modeled by Horlock’s theory is completely cancelled through a term arising 
from the distortion. 

McKeough and Graham (1980) extended the RDT based prediction to the case of a 
skewed gust which they found to be consistent with Goldstein and Atassi ’s result for zero 
spanwise wavenumber. They determined that first order corrections to the zero angle of attack 
case arise from both distortion and the streamwise component of the velocity fluctuation. They 
also show that in the case of an airfoil immersed in homogenous isotropic turbulence these 
corrections cancel. 

Atassi and Grzedzinski (1989) furthered the capability of RDT as an airfoil response 
prediction tool through the elimination of a key singularity arising from the infinite distortion of 
incoming vortical disturbances at the stagnation point. This singularity forms when applying 
RDT around a closed body, such as an airfoil and extends itself along the surface of the body at 
the exact location where the airfoil response is desired. Atassi and Grzedzinski showed that the 
irrotational field generated by the response of the airfoil must contain a canceling singular 
behavior and thus proposed a three way split of the fluctuating velocity field that accomplished 
this. Although greatly extending the capabilities of RDT in prediction of the airfoil response 
Atassi and Grzedzinski’ s method is a numerical procedure and their solution only considers the 
response at a few discrete frequencies. 

Meyers and Kerschen (1995) develop a theoretical model for the noise generated when a 
convected gust encounters an airfoil at non-zero angle of attack. The analysis uses RDT to 
consider high frequency gusts (i.e. noncompact— gust wavelength is small compared to the airfoil 
chord). They show the level of noise generated is increased by mean loading and the mean 
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loading effects scale on the airfoil total lift. Their analysis is somewhat limited, however, due to 
the high frequency nature of the theory. 

Scott and Atassi (1995) building on Scott (1990) extended the Atassi and Grezedzinski 
(1989) approach to a numerical method for solving subsonic flows with convected, three- 
dimensional vortical waves around a lifting airfoil. Their approach offers the computational 
efficiency of potential flow methods while accounting for the convection and distortion of 
incoming vortical disturbances by the mean flow field surrounding the airfoil. They compare 
their numerical solution of lift response to the Sears function and solutions of Possio’s integral 
equation investigating the effects of thickness, camber and angle of attack at particular 
wavenumber frequencies. In general, they find that thickness becomes important at higher 
frequencies. They also showed that for small angles of attack (~1°) a reduction in unsteady lift 
occurs at low frequencies. 

Reba and Kerschen (1996) develop a theory for the influence of angle of attack on the 
unsteady pressure distribution for a flat plate interacting with a high frequency gust. This work 
is extended in Reba and Kerschen (1997) to include the effects of airfoil camber. Their analysis 
utilizes RDT to consider the distortion of small amplitude vortical disturbances. The case of 
high aerodynamic reduced frequency, k is considered. The analysis shows that moderate levels 
of steady loading can significantly increase the amplitude of unsteady pressure fluctuations on 
the airfoil surface. 

1.4 Experimental Studies of the Airfoil Response 

Comparatively few experimental studies have complimented the development of airfoil 

response theories. Those that have been performed usually consist of measurements at a few 
discrete frequencies or measurement locations providing a crude representation of the airfoil 
response function. 

Jackson, Graham and Maull (1973) measured the lift spectrum on a NACA 0015 in grid 
generated turbulence (A/c = 0.417) at a Reynolds number (based on chord) of 16000 and 
compared the admittance function with two-dimensional Sears based on Liepmann’s (1952, 
1955) two-dimensional theory and his strip theory. They found good agreement with Graham’s 
(1970a) three-dimensional theory while two-dimensional and strip theory overpredicted by an 
order of magnitude and 50% respectively. Also, the measured admittance increases above the 
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three-dimensional theory at low frequencies the reason for which the authors believe is related to 
the distortion of incoming turbulence which acts to enhance lift. 

Commerford and Carta (1970, 1973) investigated unsteady pressure at 5 chordwise 
locations on a 5% thick circular arc airfoil i mm ersed in a circular cylinder wake at a chord 
Reynolds number of 400000. The cylinder wake produced a sinusoidal, two-dimensional 
upwash gust with a reduced frequency of 3.9. They found reasonable agreement at low angles of 
attack with the Sears function although noted that the phase angle variation of chordwise 
pressure was not well predicted. 

An experimental study of a circular arc airfoil immersed in a periodic wake generated by 
a pinwheel was performed by Fujita and Kovasznay (1974). Pressure measurements were made 
and compared to Meyer’s theory (thin airfoil theory) for varying angles of attack in the time 
domain. They show reasonable agreement between experiment and theory for unstalled angles 
of attack. 

Patterson and Amiet (1976a, 1976b) performed one of the few systematic unsteady 
pressure measurements which allow unambiguous comparison with theoretical developments. A 
NACA 0012 was immersed in grid generated turbulence of von Karman like spectrum (L/c = 
0.13) and the unsteady pressure measured between 15 and 70% chord through embedded 
transducers. They also measured radiated noise level with microphones located outside of the 
wind tunnel jet. One transducer on the airfoil could be moved in the spanwise direction allowing 
measurement of the spanwise correlation. Measurements were performed at 4 speeds between 
40 and 120m/s. They showed reasonable agreement between measured surface pressure spectra 
and correlations and Amiet’s (1976a) high frequency prediction at zero angle of attack. Amiet’s 
(1975) far field predictions also agreed well with far field measurements. These researchers also 
show an increase in radiated noise when the angle of attack is set to 8° and report (although do 
not show) an increase in surface pressure fluctuations which is small but measurable. 

McKeough and Graham (1980), using a configuration similar to Jackson et al. 
complimented their RDT prediction with unsteady lift measurements on a NACA 0015 at 0° and 
10° angle of attack immersed in grid generated turbulence (A/c =0.4, turbulence intensity = 
5.8%). Unsteady lift measurements were also made at 0° angle of attack in lower turbulence 
intensity flow ( A/c = 0.38, turbulence intensity = 3.7%). They showed as much as a factor of 2 
increase in mean-squared lift fluctuation at low reduced frequencies (<1) and argued that this 


13 




Chapter One 


must be a second order effect due to the canceling of first order effects arising from the distortion 
and the streamwise component of the velocity fluctuation in the case of isotropic, homogenous 
turbulence. 

Williams and Harris (1984) measured far field noise radiated by a helicopter rotor 
operating in grid and wake generated turbulence. These researchers uncovered an increase in 
peak SPL of broadband radiated noise with increasing mean loading. Also demonstrated was a 
greater rate of increase when the rotor is subject to lower turbulence levels which Williams and 
Harris feel is related to the distortion of the turbulence encountering the lifting airfoil. 

1.5 Motivation for Current Study 

The noise and vibration which result from fluid/body interaction is of great importance 
and has received considerable attention in the past 60 years. A wealth of theoretical models have 
been developed addressing the interaction of an airfoil airfoil with various representations of the 
inflow turbulence. However, relatively few experiments have been performed to validate the 
many theories which exist. In real life airfoils have thickness, camber, and often operate at 
varying angles of attack. Theories which consider these effects and are not rigorously validated 
against experimental data are as useful to design as flipping a coin. Thus, a significant need 
exists for comprehensive experimental data which can validate (or invalidate) the many theories 
in use today. 

Of particular interest in recent years has been the effect of angle of attack on the airfoil 
response. Having been postulated through various theoretical analyses, angle of attack effects 
have only recently been systematically documented. The authors presented a study of airfoil 
response to inflow turbulence in Mish and Devenport (2001, 2002). Surface pressure 
measurements were made on a 2 ft chord NACA 0015 at angles of attack varying from 0° to 20° 
in 4° increments. The airfoil was instrumented with 96 Sennheiser microphones mounted sub- 
surface (48 per side). Two inflows were considered each generated by a static grid; one with 
integral scale l/70 th of the chord and the other with integral scale l/8 th of the chord. 
Measurement of mean pressure was also performed which indicated the onset of significant 
separation at a = 16° and complete stall at a = 20°. One of the interesting results shown in this 
work was the variation in unsteady lift with angle of attack. In the presence of large scale 
turbulence (A/c = 1/8) unsteady lift was found to decrease at low reduced frequencies ( oj r < 10) 
while increasing at high reduced frequencies as the angle of attack increased. This result is 
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surprising based on the previous experimental work of Patterson and Amiet (1976a, b) and 
McKeough and Graham (1980) and is not supported by the previous theoretical work of Reba 
and Kerschen (1996, 1997). However, in small scale turbulence (A/c = 1/70) the unsteady lift 
increased across a similar reduced frequency range with increasing angle of attack. 

This work further investigates the effects of angle of attack on the airfoil response. New 
measurements of unsteady pressure in large grid flow made in an acoustically treated wind 
tunnel facility are presented. The quality of this data and Mish and Devenport (2001a, 2001b) 
and Mish, et al. (2002) data is further examined and the calculation of lift from unsteady pressure 
measurements is validated. Comparisons with Reba and Kerschens’s (1996) response theory, 
which incorporates angle of attack effects, are made, as well as additional comparisons with 
Amiet’s (1976a, b) zero angle of attack theory. This work also presents an improved model for 
the prediction of unsteady surface pressure which correctly accounts for the effects of angle of 
attack. This work is intended to be all inclusive and as such, presents some results which were 
originally shown in Mish and Devenport (2001b). 

1.6 Outline of Dissertation 

This work is presented in 6 chapters. Chapter 2 provides the background on 
measurement techniques and apparatus used to obtain unsteady surface pressure data for an 
airfoil in turbulence. An overview of airfoil response theories developed by Amiet (1976a, b) 
and Reba and Kerschen (1996) is given in chapter 3. The airfoil in turbulence set-up and data 
quality is discussed in chapter 4. Chapter 5 presents the analysis of airfoil in turbulence surface 
pressure measurements made on a NACA 0015 airfoil. Finally, conclusions and further 
considerations are given in chapter 6. 
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2. Apparatus and Instrumentation 

2.1 Single Airfoil in Turbulence 

2.1.1 Wind Tunnel 

The effects of angle of attack on the airfoil response were examined using the Virginia 
Tech Stability Wind Tunnel. The tunnel is a continuous, closed jet, single return, subsonic wind 
tunnel with a 6' x 6' test section 24 feet in length. A planform view is shown in Figure 2. 1 . 



Figure 2.1: Virginia Tech Stability Wind Tunnel 
The tunnel is powered by a 600hp DC motor driving a 14 foot propeller providing a maximum 
speed of about 60m/s. Tu nn el speed is regulated by a custom designed Emerson VIP ES-6600 
SCR Drive. Turbulence levels in the test section are extremely low, on the order of .05% or less, 
and flow in the empty test section is closely uniform (see Choi and Simpson, 1987). Although 
not originally designed as a low turbulence facility, the addition by the NACA of seven anti- 
turbulence screens, coupled with the other flow smoothing features of the tunnel, resulted in very 
low turbulence levels. The test section dynamic pressure is measured with a reference Pitot-static 
probe located 3.5" downstream of the test section entrance. 

2.1.2 Wind Tunnel Walls 

Measurements of unsteady surface pressure were made in two wind tunnel wall 
configurations; solid steel tunnel walls and acoustically treated walls. 

The standard tunnel test section walls are made of steel plate on three sides with the forth 
side made from Plexiglas (this configuration is referred to as solid walls). The acoustics in the 
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solid wall configuration are quite good as shown in figure 2.2 which is a plot of pressure sound 
level at a flow speed of 30m/s. The characteristics of this noise have been documented by 
Larssen and Devenport (1999). Acoustic noise levels were measured with a Bruel and Kjaer 
microphone fitted with a nose cone. The wind generated noise associated with this nose cone 
has not clearly been identified and therefore can not be subtracted out of the tunnel acoustic 
measurements. Thus, this measurement is a conservative estimate of the acoustic field within the 
test section and can be considered to represent the upper bound. 


PSL Vs Frequency 

Solid Walls; Microphone: Front Position on the Floor (2b) 



Figure 2.2: Pressure sound level at 30 m/s in solid wall configuration 
Despite the relatively acoustically quiet environment offered by the solid wall 
configuration, wall treatment is developed to further reduce ambient noise levels. This treatment 
consists of a series of acoustic absorbers (boxes) replacing the solid walls of the wind tunnel test 
section. The wind tunnel test section with boxes attached is shown schematically in figure 2.3. 
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Acoustic Treatment Modules 



Figure 2.3: Wind tunnel acoustic treatment schematic 
These boxes were designed to attenuate low frequency noise. The density of the material in the 
box and the box depth allow for control over the frequency range of attenuation. Thus, after 
analysis, a 20" nominal box depth is chosen which offers an expected attenuation as shown in 
figure 2.4. 



frequency [Hz] 


Figure 2.4: Expected level and frequency range of attenuation 
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The peak attenuation level of 33dB is expected to occur slightly above 100Hz with significant 
attenuation occurring from 10 to 200Hz. 

The mounting locations and preexisting tunnel configuration dictated that nearly every 
box have custom dimension. The nominal dimensions are 20" deep, 6' long, and 3' wide 
however, some boxes are smaller in depth and width due to constraints from equipment 
surrounding the test section. A cross-section of a typical box is shown in figure 2.5. 



Figure 2.5: Cross-section of typical acoustic box 
The box is made from 5/16" fiberboard, a material often referred to as MDF (medium density 
fiberboard). Each box is lightly stuffed with Owens Coming 701 insulation (1.5 lb/ft ) with a 
3 A" sheet of acoustical foam placed on top of the insulation. The acoustical foam serves as a 
barrier between the insulation and flow thus minimizing the amount of airborne insulation 
particles. A lightweight rope is strung from side to side over the acoustical foam to keep the 
insulation and foam in place. The boxes are attached to the tunnel test section with 2" flange 
aluminum angle. After the boxes are mounted to the test section sheets of aluminum screen are 
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attached which serve as the test section wall and flow surface. The aluminum screen is 1/8" 
thick with a 40% open area ratio. 

The effects of the acoustic treatment modules are documented through acoustic 
measurements with a Bruel and Kjaer microphone. The acoustic measurements of Larssen and 
Devenport (1999) made with the test section in solid wall configuration are compared with 
similar measurements in the acoustic treatment configuration. The pressure sound levels of both 
measurements are presented in figure 2.6 for a flow speed of 30m/s. 



1 10 100 1000 10000 100000 

Freq (Hz) 


Figure 2.6: Pressure sound level at 30m/s with and without acoustic treatment 
Considerable attenuation of the ambient sound field occurs between 10 and 500Hz with as much 
as 50dB of attenuation occurring between 10 and 20Hz. The spectra converge and crossover at 
500Hz with the acoustically treated spectrum rising up to lOdB above the solid wall spectrum. 
The spectra appear to converge again at 20kHz. The increase in spectral level above 500Hz with 
the acoustic treatment in place is likely the result of flow interaction with the screen. 
Additionally, the boundary layer is thicker and produces more high frequency noise with the 
screen present. 
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2.1.3 Turbulence Generating Grid 

Two grids have been developed for the Virginia Tech 6' x 6' Stability wind tunnel for the 
purpose of generating homogeneous isotropic turbulent flows for the study of unsteady airfoil 
response. The characteristics of each grid generated flow have been reported by Bereketab, et al. 
(2000) and are summarized here. The first, a square bi-planar grid with a 12" mesh size and an 
open area ratio of 69.4%, was mounted in the wind tunnel contraction 19.1 mesh sizes upstream 
of the location planned for the airfoil leading edge, at a point where the cross-sectional area was 
32% larger than that of the test section (figure 2.7). 



Figure 2. 7: Downstream view of large grid 


Acoustic foam wedges attached to the downstream side of the grid were used to minimize grid 
generated noise. The second grid, a metal weave with a 1.2" mesh size and an open area ratio of 
68.2% was mounted in the tunnel test section at a position 16.1 mesh sizes upstream of the 
planned leading edge location (figure 2.8). It should be noted that the small grid is located 
downstream of the tunnel reference Pitot-static probe which requires a correction of the pressure 
loss across the grid when used as a reference for mean pressure measurements on the surface of 
the test airfoil. 
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Figure 2. 8: View of small grid 


Table 2.1 summarizes flow speed and Reynolds number data for each grid as well as 
Reynolds number based on the chord of the NACA 0015 airfoil used in the subsequent surface 
pressure measurements. 


Table 2.1: Flow speed and Reynolds number 


Nominal Flow Speed 
(m/s) 

Small Grid Re 
(based on cell size) 

Large Grid Re 
(based on cell size) 

Equivalent airfoil 
chord Re 

30 

63000 

630000 

980000 


Three-component velocity and turbulence measurements were made with 4-sensor hot- 
wire probes over a 24"x24" cross section centered on the tunnel centerline in order to reveal the 
form of the grid generated flows. Velocity measurements were made in the plane of the airfoil 
leading edge (with airfoil removed). Turbulence intensities produced by the large grid are 
between 3% and 4% at all speeds and in all components as shown in figure 2.9. 
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Figure 2.9: Cross-sectional averaged turbulence intensities for large grid flow 


Turbulence intensities produced by the small grid are between 3.5% and 5% at all speeds and in 
all components as demonstrated in figure 2.10. 
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Figure 2.10: Cross-sectional averaged turbulence intensities for small grid flow 


The isotropy of the turbulence components in the large grid flow appears best at 30m/s 
where the difference between the smallest and largest components is less than 5% of the 
measured intensity as revealed by turbulence stress contours in figure 2.1 1. 
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Figure 2.11 (a): Large grid flow contours of Figure 2.11 (b): Large grid flow contours of 
u-component turbulence normal stress v-component turbulence normal stress 

normalized on Um normalized on Uj 


w 2 contours 




Figure 2.11 (c): Large grid flow contours of w-component turbulence normal stress 

2 

normalized on U x 

The variation in turbulence stress level appears to be random and within the uncertainty of the 
measurement. The turbulence stress field is closely homogeneous as well. 

The isotropy of the turbulence, as shown in figure 2.12 appears very good at 30m/s in 
small grid flow, the difference between the smallest and largest components at 30m/s being less 
than 1 part in 20. 
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Figure 2.12 (a): Small grid flow contours of Figure 2.12 (b): Small grid flow contours of 

u-component turbulence normal stress v-component turbulence normal stress 

2 2 

normalized on Uoo normalized on Uoo 


w 2 contours x iq ’ 3 
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Figure 2.12 (c): Small grid flow contours of w-component turbulence normal stress 

2 

normalized on Uoo 

The turbulence stress field shows some slight in-homogeneities at all flow speeds, attributed to a 
slight bow in this grid (about 0.4" at the grid center in the downstream direction), and the 
blockage produced by the frame used to support it. 

Table 2.2 presents the inf erred integral length scale for each grid at 30m/s. 

Table 2.2: Longitudinal integral scales (A) and velocity scales (u) that best describe the 2 grid 

turbulence flows. 



Large grid 

Small Grid 

u/U x (%) 

3.93% 

4.35% 

A (in.) 

3.22 

0.309 


As reported in Bereketab, et al (2000) for wavenumbers below the upper limit of the inertial 
subrange, the spectra and correlations measured with both grids can be represented using the von 
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Karman interpolation formula with a single velocity and length scale. The spectra may be 
accurately represented over the entire wavenumber range by a modification of the von Karman 
interpolation formula that includes the effects of dissipation. 


2.1.4 Airfoil Model 

Surface pressure measurements are made on an instrumented NACA 0015 with chord 
length of 2-ft and 6-ft span manufactured by Sandia National Labs. The airfoil trailing edge is 
round with a radius of curvature of 0.1075". Measurement of the airfoil geometry confirms the 
0015 shape which is calculated based on 
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Figure 2.13 presents the measured geometry plotted with the exact (analytically calculated) 
geometry. 


NACA 0015 Geometry 



Figure 2.13: NACA 0015 measured and calculated geometry 


The airfoil is manufactured as an extruded aluminum section. Minor imperfections are 
present on the surface of the airfoil consisting of scratches approximately 1/32" deep and 54 " to 
1 " long. These imperfections are located primarily out of the measurement region and cover less 
then 0.1% of the total surface. 
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The airfoil model is equipped with a removable access hatch which is 2-ft in span, the 
length of the chord and half the thickness of the airfoil. The hatch is positioned mid-span. 
Figure 2.14 reveals the airfoil with hatch removed. 



Figure 2.14: Wing with hatch removed 


The hatch is placed over the access and shimmed internally to ensure a smooth transition from 
the hatch surface to the wing. Four bolts are used to fasten the hatch to the airfoil. The hatch is 
then sealed with 0.00275" clear tape along both sides and the leading and trailing edge. Wax is 
then used to blend the tape into the surface of the airfoil. 

Mean pressure measurements are made with a 48 port scani-valve system. The scani- 
valve is connected and pressure measured with a Honeywell Model DRAL520GN pressure 
transducer interfaced with an IBM/ AT computer through a data translation DT2801-A A/D 
converter. A Dwyer series 427 Mark III handheld digital manometer (range 0-10" of H 2 O, 
accuracy ±0.01" of H 2 0) was also used to make mean pressure measurements at 23% chord. 
Tygon tubing with 3mm outer diameter is used to traverse the distance between the manometer 
and the interior of the wing. 0.075" OD copper tubing is used to transition the Tygon tubing 
from the inner surface of the wing to the outer surface. Measurement points are located 2" down 
from the hatch in a chordwise plane exactly perpendicular to the leading edge. Figure 2.15 shows 
the chordwise location of the pressure ports. 
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Figure 2.15: Chordwise distribution of mean pressure taps 
Mounting of the airfoil model is done to allow rotation of the model about the quarter 
chord for angle of attack adjustment. Two 4.6" outer diameter (3.75" ID) aluminum pipes 
protrude from both ends of the airfoil and are bolted to the models internal structure centered at 
the quarter chord. The pipes then slip through holes in the tunnel roof and floor plating. The 
pipe protruding through the tunnel roof then slides through a mounting block (Figure 2.16) 
which allows adjustment of the yaw and rolling angle. 



Figure 2.16: Roof mounting block 


Finally, the roof pipe is clamped to the mounting block (Figure 2.17) and the floor pipe is 
clamped to the floor plate. 
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Figure 2.1 7: Roof mounting structure 

Angle of attack adjustments are made by loosening the clamps and rotating the wing. 

Angle of attack is determined by setting the trailing edge to angle of attack markings on 
the wind tunnel floor. The angle of attack setting is also verified using the mean pressure data 
plotted with the potential flow solution. 

The airfoil model leading edge is located 229.2" (9.6 chords) downstream of the large 
grid as shown in figure 2.18. 



Figure 2.18: Schematic of the wind tunnel contraction and test section showing the 

location of the grid and airfoil model. 

Turbulence length scales at this location (with the airfoil removed) normalized on chord and 
leading edge radius are presented in table 2.3. 

Table 2.3: Turbulence length scales normalized on chord and leading edge radius 



Large grid 

Small Grid 

L/c 

0.134 

0.0129 

L/r c 

5.41 

0.519 
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2.1.5 Microphone System 

The airfoil model is instrumented with a unique multiple-microphone system. The 
system, illustrated schematically in figure 2.19, can accommodate up to 96 Sennheiser KE 4- 
211-2 microphones. 


1mm 









"V 




~Y“ 


96 sensors 24 4-channel amplifiers 9® signal lines 64 channel 16 bit 
embedded in model embedded in model from model acquisition 


Figure 2.19: Schematic of microphone system 

Measurements of unsteady surface pressure were made in the full 96 microphone configuration 
along with a complementary 40 microphone configuration. The 96 microphone configuration is 
shown in figure 2.20. 
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Figure 2.20: Ninety-six microphone layout 
The 40 microphone configuration is shown in figure 2.21. 



Figure 2.21: Forty microphone layout 

The Sennheiser microphones are sensitive (lOmV/Pa), provide low noise signals and a 
nominal response that is flat to within ldB from 40Hz to 10000Hz, calibratable over a larger 
range. The microphones are built with a 1mm pinhole (providing reasonable spatial resolution) 
which can be recessed a short distance (typically 0.5mm) below the airfoil surface without 
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significantly degrading the response. They can measure signals up to an SPL of 125dB and are 
insensitive to vibration. Their small physical size allows a minimum sensor spacing of about 
5mm (significantly less than the expected correlation length scales). We have built independent 
operating electronics and amplifier circuitry for each microphone using a compact modular 
approach that enables these systems to be embedded in the model close to the microphone 
locations. The system is powered by 6 custom-built quiet power supplies. It is operated, and 
measurements are made using a Hewlett Packard El 432 based data acquisition system. This 
system (which is expandable) provides simultaneous measurement of 64 channels with 16 bit 
accuracy at sampling rates up to 56kHz per channel, and has built in anti-alias filtering for each 
channel. It is important to note that this system is not limited to a 64-point measurement of the 
two-point surface pressure correlation function. Because of the spanwise homogeneity of a mean 
two-dimensional flow the 64-channel/96 microphone system can be positioned and manipulated 
to behave as a dense two dimensional array of 240 microphones. To understand how consider 
the full 96 microphone layout shown in figure 2.20. Coordinate x in this figure refers to distance 
downstream from the leading edge, and z to distance along the leading edge from the hatch edge. 
Microphones are located on both sides of the airfoil at the locations indicated. At first sight the 
microphone array appears very sparse in the spanwise direction, but this is deceptive. Since the 
flow is homogeneous in the z direction the absolute position of the sensors in this direction is 
irrelevant, and the two-point space-time correlation is only a function of their spanwise 
separation (the array defines many more spanwise separations than it does positions). The 
following is an example measurement strategy with this array 

• Simultaneously measure signals from all microphones in column 4, A1-A8, and Bl- 
B8 

• Repeat measurement except with row C, D, E 

• Because of spanwise homogeneity, this produces at least the same two-point 
correlation/spectral data as would be obtained from a simultaneous measurement 
using 240 microphones in the 2D array shown in figure 2.22 (i.e. correlation of every 
point pivoted about every point). Furthermore, the data set includes checks on the 
homogeneity assumed. 
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Figure 2.22: Effective microphone array 

A central component of any such measurement system is a reliable means for calibrating 
the microphones. For accurate measurement of the space-time correlation one needs not just 
amplitude, but full phase vs. frequency calibrations for each microphone sensor. For this purpose 
a calibrator was built, adapted from a NASA design, and is shown schematically in figure 2.23. 

To amplifier 



From compensated 
white-noise generator 
and audio amp 


Figure 2.23: Microphone Calibrator 
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The calibrator consists of a small wide-range loudspeaker mounted on one side of a shallow 
cavity built into the end of a Nylon housing. The Nylon housing is fixed into a heavy aluminum 
plate with the open end of the housing extending 2 inches above the plate. The perimeter of the 
open end of the housing is lined with a rubber gasket that seals the cavity formed by an 
aluminum cup, which slides over the Nylon housing. Each microphone is placed in a recessed 
hole on the top of the aluminum cup, which is of similar dimension as those in the wing. The 
speaker system is driven from an audio amplifier using a white noise signal, compensated for the 
response of the speaker/amplifier system, generated by the HPE1432A. The cross spectrum 
between the input signal and microphone signal is then compared with a similar measurement 
made using the B&K reference microphone. The calibration is shown schematically in figure 
2.24 where V\ (t) is the white noise signal generated by the HP1432 (and also, simultaneously 
measured by the HP 1432), P s (t ) is the pressure signal produced by the speaker signal, and Vi(t) is 
the voltage signal produced by the microphone which is feed back to the HP 1432 for 
measurement. 


White noise 

v At) 

source 

W 



m 


► 

HP 1432 Input 


Figure 2.24: Microphone calibration schematic 
The Sennheiser microphone response is obtained in the following manner: First, the speaker 
response, S p (co) [Pa/V] is determined using the Bruel and Kjaer reference microphone which has 
a known flat response, bk(a>) (0.803 V/Pa). This measurement in indicated as superscript a, 


s„h= e 


KH 

bk^co) 


■KH 



G a 

-£T [ pa/y ] 

m 


( 2 . 2 ) 


where E [ ] is the expected value. The Sennheiser microphone calibration is then found in a 
subsequent measurement (indicated as superscript b ) as 


M, (®) = 


1 


£[T(®K(®)] 's,W 



[V/Pa] 


(2.3) 


With this system a phase and amplitude calibration over the full frequency range of a 
microphone can be performed in less than 1 minute. A typical calibration performed in this 
manner is illustrated in figure 2.25. 
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Figure 2.25: Typical magnitude and phase calibrations 


The calibrator design was the result of an evolving process to achieve accurate and 
repeatable magnitude and phase calibrations. The initial calibrator design was handheld in 
nature allowing for calibration of each microphone while mounted in the wing. This design 
suffered from several flaws although primarily failing due to unreliable sealing between the 
calibrator and airfoil surface, particularly in the highly curved leading edge region, and the 
inability to consistently gather repeatable calibrations. Through further experimentation the 
Sennheiser microphone calibrations were shown to be unaffected by temperature variations 
within 70°F ±5°F. Furthermore, the calibrations show little (within the uncertainty of the 
calibration) drift with time. 

Each microphone is mounted from within the airfoil in a slightly larger diameter hole 
than the microphone. The hole is bored to within 0.5mm of the outer surface of the airfoil. A 
pin hole is drilled in the center of the bored hole and aligns with the microphone pin hole. A 
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typical cross-section through a microphone and microphone mounting hole is shown in figure 
2.26. 



Microphone signal cable 
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Figure 2.26: Cross-section through typical microphone and microphone mounting hole 


Figure 2.27 shows a microphone lying beside a mounting hole on the interior of the wing. 
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Figure 2.27: Microphone lying beside mounting hole on interior of wing 
The microphone is sealed and held in place with hot glue applied to the outer edge of the bored 
hole. 

As stated above, operating electronics are embedded in the wing below and beside the 
access hatch. Power is supplied to the electronics through 3 power cables which exit the airfoil 
through the lower (tunnel floor) mounting pipe. Data from each microphone exits the wing on 
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10 data cables via the upper (tunnel roof) mounting pipe. This configuration which ensures 
minimal data contamination from the power supplies and cables is shown in figure 2.28. 



Figure 2.28: Data and power cabling configuration 


2.1.6 Reference Acoustic System 

Measurement of the tunnel test section acoustic field was performed simultaneously with 
some measurements of unsteady pressure. A Bruel and Kjaer model 4138 (figure 2.29) 
microphone mounted on an aerodynamic stand with microphone tip in the same plane as the 
wing leading edge was used (Figure 2.30). 
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Figure 2.30: Overview of reference microphone 
mounted on stand next to wing 

The microphone was located 17.25" to the starboard side of the airfoil chord line as measured at 
0° angle of attack and 18.25" from the starboard test section wall at a height of 29.75". 

The reference microphone signal is conditioned and amplified and data then recorded by 
the HP 1432 while simultaneously recording unsteady pressure as measured by the Sennheiser 
microphones. As mentioned previously, this microphone is equipped with an aerodynamic wind 
screen which serves to protect the microphone diaphragm. The noise induced by this wind 
screen has not been undoubtedly identified and thus the acoustic levels measured will represent 
an upper threshold of the acoustic field present in the tunnel test section. 

2.1.7 Hot-wire Anemometry 

For some measurements a hot-wire probe was mounted just off the airfoil surface at the 
14% chord location. The goal of such measurement was to reveal the velocity/pressure 
correlation of the flow occurring around the airfoil and through interpretation infer the nature of 
the eddy convection and deformation. 

Velocity measurements were made using a single sensor hot-wire probe (figure 2.31). 
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Probe Prongs Stainless Steel Casing 

(length ~ 40mm) (length = 38. 1mm", diam. = 4mm) 


DUPONT 



Electrical Leads 
(length ~ 50mm) 


Figure 2.31: Diagram of hot-wire probe 

The tip of the probe consists of two stainless steel tapered prongs that hold the sensor wire 40 
mm upstream of the main body of the probe. The sensor wire is 0.8 mm in diameter and 5mm in 
length. The hot-wire is operated with a Dantec 56C17 bridge and a Dantec 56C01 constant 
temperature anemometer unit. The amplified signal is then recorded simultaneously with 
pressure measurements by the HP1432E. 

Velocity measurements were made 0.875 inches off the surface of the airfoil exactly 
above a microphone located at 14% chord and 33% span for all angles of attack including 
negative angles of attack. The mounting location was chosen to minimize interference associated 
with the probe and probe mount wake with unsteady pressure measurements. The hot-wire 
allows examination of pressure/velocity correlations on both the suction and pressure side of the 
blade. The probe is held with epoxy to a small aluminum mounting bracket. The mount bottom 
is bent at a 90° angle to form a flange which is then epoxied to the surface of the airfoil as 
displayed in figure 2.32. 



Figure 2.32: View looking down on mounted hot-wire 
The mounted probe is yawed 7° to minimize the possibility of the mount wake contaminating the 
pressure field over the remaining chordwise microphones. The mounted probe tip is in a plane 
approximately parallel to a chordwise tangent line drawn at the 14% chord location. 

The probe was operated at an overheat ratio of 1.7. Calibration is preformed a short 
period prior to tunnel entry using a jet calibrator. The probe is placed in the jet calibrator and 
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through comparison of the output voltage (Ei nv ) produced by the hot-wire in the jet flow and the 
cooling velocity sensed by the jet calibrator the hot-wires calibration is obtained. Kings 
law 

El = A l + B x U: ff (2.4) 

is used. The exponent n is taken as 0.45 and the constants A \ and B\ are determined via linear 
regression from 10 to 15 points. Flow temperatures within the test section were found to 
fluctuate ±5°F. Hot-wire signals are temperature corrected using the method of Bearman. 

2.1.8 Boundary Layer Trip and Leading Edge Bump 

For some measurements of unsteady pressure a boundary layer trip strip is used. The trip 

is made from 1" wide electrical tape cut in half with pinking shears. Two layers of tape are used 
with the top layer set 1/8" downstream from the bottom layer. The total thickness of the trip is 
0.08" with the bottom layer set 0.09% chord downstream from the leading edge. A close up 
picture of the trip strip attached to the airfoil is shown in figure 2.33. 



Figure 2.33: Close-up of boundary layer trip strip 
Measurement of unsteady pressure is also made with a small bump attached to the airfoil 
leading edge. The bump is made from yam covered with electrical tape and has a diameter of 
3/16" which extents the airfoil chord by 0.75%. 

2.1.9 Modified trailing edge 

Measurements of the unsteady pressure field were made with a modified sharp trailing 
edge to reveal trailing edge effects. This modification is incorporated using two l-l/2"wide, 
0.025" thick aluminum strips attached to the trailing edge along the entire span on both sides of 
the airfoil. The strips are attached to the airfoil with double sided tape and meet together to form 
a 0.050" thick trailing edge which extents the chord by 0.715". Figure 2.34 illustrates the airfoil 
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in test section after application of the trailing edge modification. A chordwise cross-section is 
shown in figure 2.35. 
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Figure 2.35: Cross section showing trailing edge 

extension 


Figure 2.34: Airfoil with 
modified trailing edge 


Figure 2.36 shows a dimensioned cross-section diagram of the airfoil with the trailing edge 


modification installed. 



Figure 2.36: Cross section showing dimensioned trailing edge extension 
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2.1.10 Flow Visualization 

Oil flow visualization is performed to document flow features on the airfoil surface at 
various angles of attack in large grid turbulence. Black construction paper covered with clear 
contact paper is wrapped around the airfoil over the region containing the microphones and 
attached to the airfoil with duct tape. Microphone locations are marked on the contact paper 
with a white paint pen. The oil mixture consists of 1 5 parts kerosene, 5 parts titanium dioxide 
(Ti02), and 1 part oleic acid. Figure 2.37 shows a sample oil flow visualization taken at 0° angle 
of attack. 



Figure 2.37: Example of flow visualization at 0° angle of attack 
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3. Single Airfoil in Turbulence Predictions 

Two theoretical solutions to the problem of an airfoil encountering turbulence are 
presented in this chapter. A linearized, inviscid theory developed by Amiet (1976a, b) for a flat 
plate with no mean loading is presented. This theory is currently in use by engineers and will be 
directly compared with experimental results. A more sophisticated theory, which incorporates 
the effects airfoil angle of attack is also presented here. The formulation is from Reba and 
Kerschen (1996) and uses Rapid Distortion Theory to predict the response of a flat plate in 
turbulence. Direct comparison of this theory with measured data is not appropriate; however, 
predicted angle of attack trends are used for discussion. 


3.1 Zero Mean Loading Predictions 


For the relatively easier scenario of a flat plate at zero angle of attack, Amiet (1976a, b) 
derived one of the most commonly used (likely as a result of its closed form nature) solutions to 
the problem. Amiet, reworks Filotas’ (1969) result into two solution regions, one which is valid 
at small spanwise wavenumbers ( k y ) and the other at high spanwise wavenumbers. In doing so 
Amiet was able to improve the accuracy of the complete solution to order k y . Amiet employed 
Graham’s (1970b) similarity rules which relate the general case of a skewed compressible gust to 
either a parallel compressible gust or a skewed incompressible gust, both of which are much 
easier to solve, depending on the value of the parameter <r a , where 


MK 

fk y 


(3.1) 


and 


p = Vl-M 2 . 


If <7 a < 1 the similarity is to the skewed incompressible gust, whereas, if 


(7 a > lthe similarity is to the parallel compressible gust case. 

With the two solution regions defined Amiet (1976b), noting that each of the similarity 
rules are valid outside of the regime for which they are primarily applicable (Graham 1970b), 
presents a solution for cr a < 1 derived from the parallel compressible gust case which is valid for 
small k y (< 0.25). Patterson and Amiet (1976a, b), based on Amiet (1976b) give the pressure 
difference (see appendix A for definition), with coordinate system origin at mid-span, in this case 
as 

AP(x, y,t,k x ,k y )= 2np 0 Uw 0 g(x, k x ,k y \ V) (3 .2) 
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where g(x,k x ,k y ) is the airfoil response function and defined as 

S e {K) 1 1-x ik *f°{%) 


g 


( X ’K’ k y) 


n 


1 + x 


(3.3) 


and S e (k x ) is the sears function and f a is given by 

f;~ \ 

V 2 


X(^) = (7i+i r -l)-f i ^-ln^VVrTi r -ln(l + Vl+T r )-ln2 


For large k y Amiet (1976a) derives the skewed incompressible gust solution using the 
procedure of Schwartzchild (1902) and Landahl (1961) from which the complete high frequency 
skewed compressible gust result can be obtained. The problem is formulated based on an 
infinitesimally thin airfoil representing a line discontinuity at zero angle of attack. The flow is 
assumed inviscid and incompressible with Kutta condition applied at the trailing edge. There is 
no condition specified at the leading edge. 

Due to the random nature of turbulence interacting with an airfoil it is necessary to work 
with statistical quantities such as the cross-power spectrum density of the pressure jump ( S qq ) at 

two points on the airfoil surface rather than with A P itself. The final form of S qq is obtained for 
a single sided spectrum (i.e. positive frequencies only) as a function of frequency /in Hertz as 

S^(xJ,t,J) = \6M(^]lg\*,K,X,)g(x,K,,k y )b„{K,,k y y*dk, (3.4) 


where rj = y — y',K x = 2rf fU, and is the energy spectrum of the turbulence, ft should be 
noted that a factor of n was left out of this expression in the authors masters thesis work (Mish 
2001). This has been corrected in all Amiet predictions presented in subsequent chapters. 

To compare with the measurements made using the microphone array a MATLAB code 
to predict the unsteady surface pressure field using the combined theoretical approach of Amiet 
(1976a,b) has been developed and utilizes the above formulations of S qq . Amiets (1976a, b) 
prediction scheme requires inputs related to form of the turbulence spectrum, its length scale and 
turbulence intensity, and half span of the modeled airfoil. All input parameters are summarized 
in table 3.1. 


Table 3.1: Summary of unsteady pressure prediction code input parameters 


Input Parameter 

Value 

Turbulence spectral form fO vvw ) 

von Karman as given by Patterson and Amiet (1976a) 

Free-stream velocity, Uoo 

30m/s 
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Large grid length scale, A: 

0.0818m 

Small grid length scale, A: 

0.0078m 

Large grid turbulence intensity: 

3.93% 

Small grid turbulence intensity: 

4.35% 

Airfoil half chord, b 

0.305m 


3.2 Incorporating the Effects of Turbulence Distortion 


3.2.1 Rapid Distortion Theory Background 

A real airfoil with thickness and at some angle of attack will have a mean potential field 

associated with its presence in the flow. As such, turbulence encountering the airfoil will be 
modified as a consequence of the mean field velocity gradients. The modification of the 
turbulence takes place through stretching, compressing, and tilting of the vorticity vector. Such 
interaction can be described through the vorticity transport equation which is derived by 
considering the curl of the Eulers equation and given here as, 


D Qj 

Dt 


dU i 

= Q J 


dx. 


(3.5) 


This equation is often referred to as the rapid distortion equation. This equation is, of course, an 
inviscid representation of the distortion. As such, it is valid only when the distortion occurs so 
rapidly that the interaction between turbulence and the influence of viscous dissipation can be 
neglected. 

This method of modeling the evolution of turbulence is referred to as Rapid Distortion 
Theory (RDT). Batchelor and Proudman (1954) performed some of the pioneering work in the 
area of RDT, developing the appropriate conditions for its applications and the techniques 
required to apply this theory to a real flow. When considering the response of an airfoil to a 
turbulent inflow it is important to consider the effects of the mean field induced distortion of the 
inflow turbulence. Therefore, several researchers have developed theories for predicting the 
response of an airfoil to turbulence using the RDT approach. 

The theory considered here was developed by Kerschen and Myers (1983), Myers and 
Kerschen (1995), and Reba and Kerschen (1996) for application to a flat plate airfoil subject to 
some mean loading (i.e. at some angle of attack). A solution is derived for the unsteady pressure 
distribution on the airfoil surface due to interaction with a three dimensional time harmonic gust. 
The airfoil is assumed to be a zero thickness flat plate at some mean angle of attack to a subsonic 
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inflow. The distortion of the inflow by the mean field velocity gradients is accounted for using 
RDT. Singular perturbation expansions are developed for the case of Mach number of order 1, 
wherein the acoustic wavelength is short compared to the airfoil chord and for the case of small 
Mach number, in which the acoustic wavelength is long compared to the airfoil chord. As such, 
the airfoil-gust interaction analysis is simplified in various asymptotic limits at high and low 
acoustic frequencies, where the acoustic frequency is given as 

k a = w a (o r (3.6) 


with w = 


P 2 


V 2 


,P = fi -Ml ,h is the spanwise wavenumber. 


Note, the notation used here is consistent with that presented by Reba and Kerschen (1996) and 
is used and defined only is this section. For the case of M x « 1 (the case for the present 
experiment, M x = 0.087), the low acoustic frequency limit is applied. In doing so the gust may 
be two or three dimensional; however, due to the relationship between the acoustic frequency 
and spanwise wavenumber a three dimensional gust is limited to only very small spanwise 
wavenumbers. Reba and Kerschen (1996) develop a FORTRAN computer routine the theory 
outlined above. 

Reba and Kerschen (1996) give the inflow as 

v' (co r , k n ,L)= eU„ ( A t , A n , A 3 y^Kv^-t ) (3.7) 


where U,£ is taken to be the upwash amplitude referred to as w 0 . The parameters of equation 3.3 
are given in table 3.2 


Table 3.2: Parameters used in Reba and Kerschen (1996) airfoil response theory 


Parameter 

Comment 

Non-dimensionalization 

£ 

Dimensionless small number 

- 

U X, 

Freestream velocity 

- 

C ^tt A n , A 3 ) 

Streamwise, normal to planform, and spanwise non- 
dimensional gust amplitudes 

- 

(D r 

Reduced frequency 

II 

"1 

b 

Half chord, c/2 

- 


Non-dimensional mean flow velocity potential 

(f)p/U cob 

¥ 

Non-dimensional mean flow stream function 

yfp f^oo /U cob 

■Ts 

Non-dimensional spanwise coordinate 

x p /b 

k 

Non-dimensional streamwise wavenumber 

k t b 

(Dr k n 

Non-dimensional normal to planform wavenumber 

- 
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CO r k 3 

Non-dimensional spanwise wavenumber 

- 

t 

Non-dimensional time 

tpb/Uoo 

Poo 

Prandtl-Glauert transformation 

(1-M 2 ) 1 ' 2 


A modified pressure, p is calculated by the code for the pressure and suction side of the airfoil 
which, for a given set of gust amplitudes, is a function of the wavenumber vector, chordwise 
position, Mach number, and angle of attack. The modified pressure is related to the physical 
pressure, p' on the surface of the airfoil by 

p\x,k n ,k 3i t) = -p(x,a) ri k ni k i ) pU x w o e ar{k,x, ~ t] (3.8) 

where w a is the upwash amplitude. The code also uses a FORTRAN IMSL library function for 
the evaluation of the complex error function. This function produces significant errors when the 
reduced frequency becomes large. Therefore, all analysis with this prediction scheme is limited 
to 1 < o)r <10 which seems to produce reasonable results. 

3.2.2 Derivation of Cross-spectral Density of Surface Pressure Fluctuations 

Reba and Kerschen (1996) do not present an explicit expression for the airfoil response 

function (i.e. transfer function). Additionally, their program is designed to compute the modified 
pressure for one set of inflow wave amplitudes and wavenumbers. Therefore, in order to 
estimate the surface pressure cross-spectral density using this code for the experimental inflows 
researched here it was necessary to implicitly determine the airfoil response function and 
integrate over all wavenumbers. This is accomplished by considering the response of the airfoil 
to a unit amplitude gust. As such, the transfer function is found in terms of the modified 
pressure, p ua as (keeping with the notation of Reba and Kerschen) 

g(x,G) r ,k n ,k 3 ) = -p ua (x,a> r ,k n ,K) P U r w „ ( 3 - 9 ) 

with W 0 = 1 . Summing over all wavenumbers, the relationship between the surface pressure, p' 
and the inflow can then be written in terms of the transfer function as 

OO 00 00 

p'(x,x 3 ,t)= | | ^ g[x,co r ,k n ,k i )w 0 e COr(kiX ^dk n dk i dco r (3.10) 

—00 —00 —00 

where W 0 is the upwash amplitude. The Fourier transform with respect to time can now be 

performed to give the surface pressure in terms of frequency such that 

00 00 00 00 

p'(x,x 3 ,co r )= J | | | g (x, co r , k n , k-, ) w 0 e a>r ^ Xi ~ t ] e ""’'dkplLd a> r dt (3.11) 
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Assuming the dependency of the spanwise wavenumber on reduced frequency is removed 1 and 
substituting equation 3.9 for the transfer function, g, equation 3.11 becomes 

OO 00 

p'(x,x i ,a r ) = U x pj J -p ua (x,a> r ,k n ,k 3 )w 0 e ihXi dk n dk 3 (3.12) 

—00 —00 

where k 3 = — . The pressure jump (or difference), A P is now formed by considering the 

difference between the pressure on the upper and lower surface of the airfoil 
A P(x,x 3 ,co r ) = p u '-p L ' 

”” . (3.13) 

= U ao P J J - L Puna (-T <o r ,k n ,k 3 ) - p Lua (x,(0 r ,k n , L ) J w 0 e‘ 3X3 dk n dL 

—00 —00 

where p Uua and p Lua are the modified pressures computed for a unit amplitude inflow on the 

upper and lower surfaces respectively. The cross-spectrum of pressure difference fluctuations 
can now be computed by taking the expected value of the pressure difference times its conjugate 

Sqq (x,x',x 3 ,x 3 ,co r ) = lim^^ ^AP(x,x 3 ,<» r ) AP(x',xj,<y r ) . (3.14) 

Noting that the only non-deterministic quantity of equation 3.13 is the upwash amplitude, \v„ 
which can be related to the upwash velocity spectrum, <j) ww as 

E \w 0 {co r ,k 3 )w 0 {(Or,k' 3 jA = —8{k 3 -k\) ( l) ww {(o r ,k 3 ) (3.15) 

L J ft 

where T = R/Uoo and R is the spatial region over which the inflow is considered (both T and R are 
normalized with the free-stream velocity), the pressure difference cross-spectrum becomes 

00 00 

spx,x\n,^Mu„pf J \(p u ,„-pP(PuJ-p L jU,M,XV % ''<lk,dk i (3.16) 

—00 —00 

where rj = x 3 - x'3. To avoid confusion it is worth noting, explicitly that equation 3.16 is for a 
double sided spectrum in frequency (-00 < co, < co) and therefore is multiplied by 2 when 
calculated to obtain a single sided spectrum for comparison purposes. 


1 The FORTRAN code developed by Reba and Kerschen (1996) was modified to remove the dependence of the 
spanwise and normal to planform wavenumber on frequency by dividing both parameters by the reduced frequency. 
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3.2.3 Application of Reba and Kerschen Surface Pressure Prediction Code 

As a check to ensure the cross-spectral density, S qq is calculated correctly a comparison 
with Amiet’s (1976a, b) theory is performed. Both theories ideally should converge to the same 
solution at zero angle of attack. Since it is not possible to integrate over the same spanwise 
wavenumber range in both theories (i.e. Kerschen’ s theory is limited to low acoustic frequencies) 
a comparison at one set of wavenumbers is performed. The following are the input parameters 
for both theories; 

Table 33: Parameters used in Kersch en/A m iet comparison 


Spanwise separation 

Turbulence intensity 

Integral length scale of turbulence 

Chordwise locations 

Figure 3.1 is a plot of the auto-spectral level for Amiet (1976a, b) and Reba and Kerschen’s 

(1996) theories as a function of chordwise location. Note that in the y-axis label G ApAp is the 
auto-spectrum of the pressure difference calculated in Hertz and is related to S* as 

<W/)=f r s G x=x '’ r >’^) 3 - 17 

^ oo 


JL 

u/Uoo 

A 

x/c 


Parameter 


Moo 

k 



Comment 


Mach No. 

Streamwise wavenumber 


Spanwise wavenumber 

Normal to planform wavenumber 
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Figure 3.1: Comparison of Amiet (1976a, b) and Reba and Kerschen (1996) theories 

Figure 3.1 shows excellent agreement between the two theories up to 4% chord at which point 
Reba and Kerschen’s theory begins to diverge slightly. The divergence of the two theories 
beyond 4% chord is likely related to the numerical integration scheme employed in Reba and 
Kerschen’s code. 
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Having confidence in the formulation of S qq the above procedure is used to investigate the 
effects of turbulence distortion on unsteady pressure and lift. The following serve as input 
parameters to Reba and Kerschen’s FORTRAN code; 

Moo = 0.087 

k = i, where i = 1, 2, 3. . . 10 

1 q[(./-I)V 3 ] 

k 3 = , where j = 1, 2, 3. . .41 and A & 3 = 0.01 

k 

k n = , where l = 1, 2, 3. . .30 and A k n = 0.06 

k 

7 = 0 

a = 0 °, 2 °, 6 ° 
x/c = 1, 2.5, 4, 6 , 9, 14% 
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4. Single Airfoil in Turbulence: Experimental Details and Data 
Quality 

This experiment sets out to measure the inviscid response of an airfoil in turbulence using 
an array of microphones mounted sub-surface. In this chapter the experimental details are 
provided and the quality of the measured unsteady pressure data examined. 

4.1 Measurement Conditions and Data Sampling 

Measurements of unsteady pressure were made in two wind tunnel configurations. A set 
of measurements was made with the Virginia Tech Stability Wind Tunnel in its solid wall 
configuration (i.e. steel plates form test section walls). A complimentary set of measurements 
was made with the steel walls replaced by acoustically treated walls. 

4.1.1 Solid Test Section Wall Test 

Measurements of unsteady surface pressure were made in solid wall configuration using 
an array of ninety-six microphones (forty-eight per side with each side a mirror image of the 
other). Measurements were made at angles of attack, a = 0°, 4°, 8°, 12°, 16°, 20° (figure 4.1 
defines positive angles of attack). The nominal flow velocity was 30 m/s corresponding to a 
chord Reynolds number of 1.17 x 10 6 . Measurements were made with the airfoil immersed in 
small and large grid turbulence. Two sampling rates were used to gather data in large grid 
turbulence; 1.6 kHz and 12.8 kHz. Data was sampled at 12.8 kHz in small grid turbulence. Five 
hundred records of data were taken with a record length of 2048. 

Simultaneous measurement of all ninety-six microphones was not possible due to the data 
acquisition systems sixty-four channel capacity. Therefore, measurements were made in three 
microphone blocks with the longest chordwise row of microphones common to each as discussed 
in chapter 2 and shown in figure 2.22. Block one measured the first two spanwise rows (row A 
& B), block two the third and forth spanwise rows (row C & D), and block three the fifth 
spanwise row (row E). Some measurements also used a forth block of microphones which was 
composed of rows B & D. 

Several subsets of unsteady pressure data were taken under alternate configurations. To 
aid in assessing noise in unsteady pressure measurements, data was taken with the turbulence 
generating grid removed. Two types of measurements were made in this configuration. One 
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matched the flow velocity of 30 m/s (nominally) which aids in assessing contamination 
associated with boundary layer pressure fluctuations. This measurement is referred to as 
‘smooth flow’ in subsequent discussion. The other matched the tunnel fan speed (RPM) 
achieved when the grid is installed which aids in determining the extent of ambient acoustic field 
contamination. Additionally, as discussed in chapter 2, measurements of unsteady pressure were 
made with a modified trailing edge and simultaneous velocity/pressure measurements were also 
performed. Table 4.1 summarizes the angle of attacks, microphone configuration, sampling 
rates, and grids used in each solid wall measurement. 


Table. 4./: Solid Wall Test Configuration Details 


Measurement 

Angle of attack 

No. of 

Mies 

Microphones Used. 
Definitions refer to 
figure 2.20 

Grid 

Sampling 

Rate 

Unsteady pressure 

0°, 4°, 8°, 12°, 
16°, 20° 

96 

Column 4 
Row A-E 

Large 

1.6kHz 

Small 

12.8kHz 

Column 4 
Row A & B 

Large 

12.8kHz 

V elocity/Pressure 

-16°, -8°, 0°, 
8°, 16° 

54 

Column 4 
Row B & D 

Large 

1.6kHz 

Small 

12.8kHz 

Modified trailing 
edge 

0°, 4°, 8°, 12°, 
16°, 20° 

54 

Column 4 
Row B & D 

Large 

1.6kHz 

Small 

12.8kHz 

Matched flow 
speed 

0°, 4°, 8°, 12°, 
16°, 20° 

54 

Column 4 
Row B & D 

No grid 

1.6kHz 

Matched RPM 

0° 

54 

Column 4 
Row B & D 

No grid 

1.6kHz 

No grid 

12.8kHz 


Mean pressure measurements were also made using the forty-eight port scani-valve 
system described in chapter 2. These measurements were made at a = 0°, 4°, 8°, 12°, 16°, 20°. 


4.1.2 Acoustically Treated Wall Test 

A complimentary set of unsteady pressure measurements were made with the Virginia 
Tech Stability Wind Tunnel solid walls replaced by acoustic treatment modules. These 
measurements were performed with an array of forty microphones configured as discussed in 
chapter 2 and shown in figure 2.21 (i.e. column 4 and row B). The instrumented airfoil was 
immersed in large grid turbulence only for this set of measurements with data taken over a range 
of angles of attack, a, from 0° to 20° in 2° increments. Two sampling rates were used for these 
measurements; 1.6 kHz and 25.6 kHz. One-thousand records were taken with a record length of 
2048. 
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Additionally, measurements were performed at a = -4° and -8°. A boundary layer trip 
was incorporated into a set of measurements covering a = -8°, -4°, and 4°. Lastly, a leading edge 
bump was installed and unsteady pressure measured at a = -8°. For this set of measurements a 
sampling rate of 1 .6 kHz with a record length of 2048 was used to take 200 records of data. 
Mean pressure measurements were made at 20% chord to ensure correct angle of attack 
alignment. 

4.1.3 Coordinate System 

The coordinate system used in the description of all single airfoil in turbulence data is 
shown in figure 4.1 and 4.2. The origin is in the plane of the leading edge with x directed 
downstream at mid-span, y directed up along the leading edge, and z completing the system by 
the right-hand rule. 

4.2 Measurement Uncertainty 

4.2.1 Uncertainty in Pressure Spectra 

An estimate of the uncertainty in unsteady pressure measurements is made based on 
repeated measurements. At any one angle of attack measurements were made with three 
microphone blocks with the longest chordwise row of microphones common to each block. By 
considering the three repeated measurements of the longest chordwise row the uncertainty in the 
pressure difference spectra can be obtained at each chordwise location. Table 4.2 presents the 
uncertainty in pressure difference spectra, s r for x/c = 1 and 4% (representative locations) in 
decibels (dB) based on repeated measurements at a = 0° in large grid flow at 20 to 1 odds (95%). 


Table 4.2: Uncertainties in unsteady pressure spectra 


(Dr 

£,. atx/c = 0.1 (±dB) 

£ r at x/c = 0.4 (±dB) 

Large Grid 



1 

0.58 

0.92 

10 

0.49 

0.40 

35 

0.81 

0.18 

Small Grid 



10 

0.58 

0.46 

100 

0.50 

0.57 

300 

0.36 

O 

U> 

oc 
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The uncertainty in auto-spectra is further reduced through spanwise averaging. The 
reduction in uncertainty based on statistical convergence is given by Bendat and Piersol (1986) 
as 


e r 

Vtv 


(4.1) 


where N = 8 (the number of microphones in a spanwise row). As such, spanwise averaging 
reduces spectral certainty sixty-five percent giving a maximum uncertainty of 0.35dB at x/c = 
4% for oj r = 1 . 


4.2.2 Uncertainty in RMS Pressure 

Mean-square pressure is calculated by integrating the pressure spectrum such that 

__ 00 

P 2 =\S qq (x = x\j] = 0,o))da> (4.2) 

o 

where S qq is the surface pressure cross-spectral density. Values of the RMS pressure, p' are 
calculated from the square root of/? 2 . Using the analysis of Kline and McClintock (1953) the 
uncertainty in the p 2 is estimated as 

1 Yj S 2 m ( OJ >) A ° j2 

£—=£ ” M (A 3t 

p 2 ° spectrum N V / 

i = 1 

where S (<s? ; ) and A co i are the spectral value and bin width at the z'th position. Here the error in 


the spectrum, s S p ectnim is taken as s c from equation (4.1). The uncertainty of the RMS pressure is 
given by Kline and McClintock (1953) as 


e 


p' 



(4.4) 


Root mean square pressure uncertainty,^, values at chordwise locations from 1% to 14% are 

evaluated at each angle of attack on the pressure and suction side of the airfoil. These values are 
found to be less than 3%. 
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4.3 Basic Flow Definition 

4.3.1 Mean and Unsteady Blockage Effects 

The combined effect of solid and wake blockage produce at most slightly more than 1% 
change in velocity in both grid flows which is expected to negligibly affect unsteady pressure 
measurements. Streamline curvature restriction produced up to a 0.38° increase in angle of 
attack occurring in the small grid flow at a = 16° and as much as a 0.042 reduction in Cl 
occurring again in the small grid flow at a = 20°. The changes implied by the streamline 
curvature restriction correction are considered small and will not adversely affect unsteady 
pressure measurements. Additional details on mean blockage effects are available in Mish (2001) 

The effects of unsteady blockage or, what is also referred to as unsteady wall 
interference, have been investigated by Karasoy (1990). Karasoy measures the unsteady surface 
pressure of an oscillating airfoil at several oscillation frequencies and at varying values of h/c, 
where h is the tunnel test section half height/width and c is the characteristic dimension of the 
body (taken to be the airfoil chord). These results show that for small values of h/c (< 0.5) a 
large rise in unsteady pressure occurs for co r < 1 . As h/c increases above 1 the unsteady surface 
pressure is minimally affected. In the case of the present experiment h/c = 1.5 and, as such the 
unsteady pressure measurements are believed to be largely unaffected. 

4.3.2 Mean Pressure 

Mean pressure measurements are summarized here (discussed in detail in Mish 2001) 
through the presentation of mean Cl versus angle of attack for both large and small grid 
measurements. Figure 4.3, adapted from Mish (2001), presents lift coefficients for a = 0°, 4°, 8°, 
12°, 16°, and 20° calculated from integrating mean pressure measurements. Also shown on 
figure 4.3 are the mean Q values for various measurements made by McKeough (1976) with a 
NACA 0015 in various scale and intensity turbulent flows. The theoretical lift curve, 2na is also 
plotted on figure 4.3. In both the large and small grid flows there is a region of unstalled flow 
which persists up to a = 16°. With the airfoil immersed in small grid flow complete stall does 
not occur at a = 20°. This is likely a result of the small scale grid generated turbulence 
interacting with the turbulent boundary layer. High levels of kinetic energy are brought into the 
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boundary layer since the scale of the grid generated turbulence in on the order of the turbulent 
boundary layer scale. In the large grid flow the airfoil stalls at a = 20°. 

4.3.3 Stagnation, Transition, and Separation Points 

The transition, stagnation, and separation locations are recapped here (discussed in detail 
in Mish 2001) through figure 4.4 (small grid) and 4.5 (large grid). These figures, which are 
adapted from Mish (2001), show a cross-section of the airfoil with the locations of transition, 
stagnation, and separation denoted by symbols along the edge of the airfoil. In the small grid 
turbulence transition, determined from high pass filtered pressure time series, was found to occur 
forward of the most leading edge microphone (1% chord) and is therefore taken to occur at the 
leading edge. Based on examination of mean pressure data, a small laminar separation bubble 
may develop over the leading edge region at a = 16° and 20°. No independent verification of this 
bubble is offered. The large grid transition locations are shown to shift down the chord on the 
pressure side as the angle of attack is increased while moving towards the leading edge on the 
suction side. Mean pressure data also suggests the occurrence of a laminar separation bubble at 
a = 16° in the large grid flow. Although not presented here, Mish (2001) shows a comparison of 
the mean pressure at a = 16° angle with panel method calculations which show a loss of pressure 
over the first 5% chord. Oil flow visualization is used to provide additional evidence on the 
presence of this bubble. 

4.3.4 Large Grid Flow Visualization 

Oil flow visualization is performed over a range of angles of attack including, a = 0°, ±4°, 
±8°, ±12°, 16°, 20° in large grid turbulence only. The oil flow visualization is performed by 
wrapping black contact paper around the airfoil as discussed in chapter 2. The contact paper will 
of course moderately alter the shape of the airfoil and smoothes out surface irregularities. 
Therefore, flow phenomena which are be triggered by local airfoil geometry may not be present 
when oil flow visualization is in progress. Hence, oil flow visualization, as with many types of 
flow visualization requires cautious interpretation of results. As such, the conclusions drawn 
here are somewhat ambiguous and, where possible, are corroborated with independent data. 

Figure 4.6 presents oil flow visualization results at a = 0° on the access hatch side of the 
airfoil (4.6a) and solid side of the airfoil (4.6b). A thinning of particles is shown to occur, 
between 1 1% and 15% chord on the hatch (pressure side when angle of attack is positive) side of 
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the airfoil and between 8% and 12% chord on the solid side. This is likely a result of boundary 
layer transition. That is, as transition occurs there is a rise in skin friction which in turn carries 
additional particles out of the transition region. The transition points noted by Mish (2001) 
based on analysis of unsteady pressure data are near 6% chord on the hatch side and 8% chord on 
the solid side. The slight disagreement may be the result of the contact paper used during flow 
visualization. As mentioned above, the contact paper will tend to smooth out small 
imperfections on the surface of the airfoil. These imperfections would otherwise act to 
destabilize the flow and cause transition to occur closer to the airfoil leading edge. 

As the angle of attack is increased from a = 4° to 16° a darker colored region persists (i.e. 
thinning of particles in the oil mixture implying an increase in skin friction) over approximately 
3% of the chord on the suction side of the airfoil that may be related to an unsteady laminar 
separation bubble (denoted in figures 4.7-4. 1 1 by the two red lines). This region, staying nearly 
the same size moves towards the leading edge with increasing angle of attack. The laminar 
separation bubble, if truly present, likely develops as follows; first the laminar boundary layer 
separates just downstream of the suction peak where the pressure is strongly adverse. Next, the 
separated laminar shear layer quickly becomes turbulent and subsequently reattaches a short 
distance downstream. The turbulent reattachment is believed to cause the thinning of oil 
particles shown in figures 4. 7-4.1 1. Analysis of mean pressure data also suggests the occurrence 
of a laminar separation bubble at a = 16°. Furthermore, the work of McKeough (1976) showed a 
similar phenomenon occurring on a NACA 0015 immersed in grid generated turbulence. 
McKeough immersed a 0.3m NACA 0015 in turbulence characterized by a length scale 40% of c 
and intensity of 6%. He observed, through tuft flow visualization, the development of a small 
leading edge laminar separation bubble as the angle of attack was increased, ultimately 
culminating in combined lead/trailing edge stall. The lack of a clear separation or reattachment 
line in the oil suggests that the laminar separation bubble is highly unsteady. This would be 
consistent with bubble formation in turbulent flow. That is, as the loading on the airfoil 
fluctuates as a consequence of interacting with the inflow, the size and position of the bubble 
will be altered blurring the separation and reattachment points. At a = 20° (figure 4.1 1) the flow 
appears to separate just before 1% chord causing the airfoil to completely stall. 

Figures 4.6-4.11 also show the approximate stagnation point location. These locations 
agree well with the mean pressure data presented in Mish (2001). Table 4.3 presents the 
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approximate locations of the stagnation point and turbulent reattachment of the conjectured 
laminar separation bubble at each angle of attack along with the theoretical stagnation point 
location (calculated based on inviscid theory using a vortex panel code). The location of 
turbulent reattachment of the speculated laminar separation bubble is approximated based on the 
observed oil thinning. 

Table 4.3: Stagnation point and laminar separation bubble locations for large grid 


measurement 


a (deg) 

Stag Pt. 
Theory (% c) 

Stag Pt. Flow 
Vis. (% c) 

Turbulent reattachment of 
speculated laminar separation 
bubble location (% c) 

0 

0 

0 

No separation 

4 

0.4 

No data 

6-9 

8 

1.6 

1.7 

2.5-6 

12 

3.6 

2.7 

2.6-6 

16 

6.7 

4.2 

0.9-2. 3 

20 

10.5 

6.0 

separation 


4.4 Quality of Unsteady Pressure Data 
4.4.1 Spanwise Homogeneity 

It is important to verify that measured unsteady pressure fluctuations are indeed 
homogeneous across the span as would be expected based on the spanwise homogeneous inflow. 
The homogeneity is checked by comparing auto-spectra from each spanwise location for any 
given chordwise location. Such comparisons are possible for chordwise positions x/c = 1, 2.5, 4, 
6, 9%. 

Small grid surface pressure auto-spectra (G pp ) from each spanwise location are presented 
in figure 4.12 at a = 0° for chordwise positions x/c = 1, 2.5, 4, 6, 9% on the suction side of the 
airfoil. The spectra at x/c = 1% are banded together within ~3 decibels (dB) for all spanwise 
locations except y/c = 33%. The calibration for the microphone at this location is incorrect and, 
as such, this microphone has been removed from subsequent analysis. Moving to x/c = 2.5%, the 
spectra across the span fall agreeably on top of each other with ~2dB of spread. Similar 
agreement is shown at x/c = 4 and 6%. The homogeneity of spectra at x/c = 9% is within 4dB 
however, microphones at locations y/c = 24 and 35% fall decades below the other locations as a 
result of bad microphone calibrations. These two microphones were also removed from 
additional analysis. 
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Small grid pressure side auto-spectra from each spanwise location are presented at a = 0° 
for chordwise positions x/c = 1, 2.5, 4, 6, 9% in figure 4.13. The spectra from different spanwise 
stations compare to within 3dB at all locations except y/c = 27%. The source of this spectrums 
rise (~3dBrelative to its spanwise companions is not clear; however may be related to flow 
asymmetries produced by small asymmetries in the leading edge. The spectra agree well at x/c = 
2.5 and 4% with ~3dB of spread across the span. Moving to x/c = 6%, spectra vary across the 
span by as much as 7dB. This poor agreement in spanwise spectra persists at x/c = 9% with 
~9dB of discrepancy. The source of this considerable variation at both x/c = 6 and 9% may be 
related to asymmetries in the flow. Inspection of microphone calibrations shows them to be 
consistent at each spanwise location. Flow asymmetries may arise from asymmetric boundary 
layer transition, wakes produced by microphone mounting holes, or leading edge asymmetries. 

Comparisons of suction side spectra at different spanwise stations for the large grid 
turbulence are presented in figure 4.14 with the airfoil at a = 0°. Spanwise data comparisons are 
made for chordwise locations x/c = 1, 2.5, 4, 6, 9%. Spectra at x/c = 1% are within ldB across 
the span. Similar agreement is shown at x/c = 2.5% with exception to the spectrum at y/c = 98%. 
The microphone here suffers from a bad calibration and is removed from remaining analysis. 
Spanwise spectra at x/c = 4 and 6% agree within ~2dB. This agreement continues at x/c = 9% 
however, two microphones at y/c = 24 and 27% suffer from bad calibrations and are also 
removed from further data analysis. 

Large grid pressure side auto-spectra from each spanwise location are presented at a = 0° 
for chordwise positions x/c = 1, 2.5, 4, 6, 9% in figure 4.15. Spectra at x/c = 1% agree within 
~3dB except for the spectrum at y/c = 46%. A flow asymmetry, likely produced by the leading 
edge, seems to be the source of spectral level reduction at this location. The agreement of 
spanwise spectra improves slight to within ~2dB at x/c = 2.5, 4, and 6%. At x/c = 9% spanwise 
spectra are within ~3dB for co r < 10. Above co r = 10 the spectra spread out with as much as KMB 
of difference occurring across the span. This high frequency variation is likely related to 
asymmetric boundary layer transition. 

In summary, the homogeneity of spanwise unsteady pressure measurements is good with 
exception to several measurement locations corrupted by bad microphone calibrations. These 
locations have been removed from further data analysis. With spanwise consistency 
demonstrated, the data is averaged across the span to aid in reducing uncertainty (a full 
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explanation of the formation of cross-spectral density of pressure fluctuations and the inherent 
averaging utilized is presented in appendix A). The remaining analysis will employ this reduced 
uncertainty data set, where appropriate. 

4.4.2 Effects of Acoustic Treatment and Data Repeatability 

Recalling that two independent measurements (solid wall and acoustic walls) of unsteady 
pressure were performed only with the airfoil in large grid turbulence, this section discusses 
repeatability and acoustic treatment effects related to the large grid data only. 

Zero Angle of Attack: Measurements of unsteady pressure were made in two wind tunnel test 
section configurations (solid and acoustically treated walls). Comparison of spectra from these 
measurements is useful for determining the repeatability of the measurement and assessing the 
effects of acoustic treatment. Of chief concern is the presence of local flow phenomena which 
may not be highly repeatable and therefore, produce differences in the data sets. The effect and 
repeatability of laminar separation bubbles is of particular interest given the evidence of such 
features in oil flow visualization. As such, comparisons of pressure spectra from both 
measurements are made in figure 4.16 and 4.17 for the suction and pressure sides of the airfoil 
respectively. Spectra are presented for a = 0° at x/c= 1, 2.5, 4, 6, 9, 14%. 

The suction side comparison (figure 4.16) consistently shows the solid wall pressure 
spectra residing above the acoustic wall spectra. At x/c = 1% the agreement is within ~1.75dB 
across the frequency range. However, the agreement deteriorates slight at x/c = 2.5% with ~3dB 
of difference occurring over all frequencies. Chordwise positions, x/c = 4, 6, and 9% show 
agreement within ~2dB across all frequencies. Spectra at 14% chord are within ~2dB for oj, < 7 
and then diverge by as much as ~7dB over higher frequencies. The divergence at high 
frequencies may result from transition occurring at slightly different chordwise locations in each 
measurement. 

Spectra measured on the airfoil pressure side for each wall configuration are presented in 
figure 4.17. As with the suction side, the solid wall spectra are consistently situated above the 
acoustic wall spectra at all chordwise locations. The leading edge measurement location, x/c = 
1% shows ~3dB of difference for co r < 10. The agreement improves to nearly perfect as the 
frequency increases. The poor agreement at this location may be the result of a local flow 
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phenomenon, such as a transient laminar separation bubble. Evidence of such a bubble is shown 
in oil flow visualization. The agreement improves to within ~ldB at x/c = 2.5 and 4% over all 
frequencies. At x/c = 6% agreement is within ~2dB for co r < 10, however, deteriorates to ~4dB 
with increasing frequency. Spectra are within ldB for co r < 10 at x/c = 9%. Due, perhaps, to 
differences in chordwise transition location, the agreement gets worse with increasing frequency. 
Similar type agreement is shown at x/c = 14%, however, for co r < 10 nearly 3dB of difference in 
spectral level occurs. 

Eight Degrees Angle of Attack: Comparisons of spectra measured in the presence of solid and 
acoustic walls are also made at a = 8° for chordwise locations, x/c = 1 and 2.5%. Figure 4.18 
and 4.19 present spectra measured on the suction and pressure side of the airfoil respectively. 
On the suction side (figure 4.18) the spectrum measured with solid walls resides above spectrum 
measured with acoustic walls. The difference is nearly 2dB at both chordwise locations, x/c = 1 
and 2.5% across the frequency range and is similar to that shown at a = 0°. On the pressure side 
(figure 4.19) at x/c = 1% the solid wall spectrum is above the acoustic wall spectrum by as much 
as 3dB across the frequency range. The agreement improves at x/c = 2.5% with no more than 
ldB of difference occurring; however, the acoustic wall spectrum has now shifted above the 
solid wall spectrum. 

Discussion on Measurement Repeatability: The repeatability of this measurement is 

reasonable; however, differences are shown on both sides of the airfoil at a = 0° and 8°. Because 
differences occur across the range of measured frequencies, it is unlikely such differences arise 
solely from the acoustically treated test section walls. The acoustic treatment was shown to 
produce significant attenuation at low frequencies and therefore, would be expect to impact 
unsteady pressure measurements in this range. The other possible sources for the observed 
differences are; 

1 . Airfoil location 

2. Turbulence generating grid 

It is possible that the airfoil was shifted upstream or downstream by nominally 5cm between the 
two measurements thereby, locating the airfoil a slightly different distance from the turbulence 
generating grid. As such, the energy in the inflow at the leading edge would change between the 
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two measurements producing pressure fluctuations of different magnitude. The data presented 
suggests (since solid wall spectra reside above acoustic wall spectra) the airfoil was located 
slightly further downstream during the acoustically treated test section measurement. The other 
possibility is the turbulence grid. The grid was removed and then reinstalled for the second 
measurement and in doing so may have been reinstalled slightly differently. These differences 
could alter the turbulence produced and cause the small differences in unsteady pressure 
measurements. 

Lastly, the large difference observed at 1% chord on the pressure side of the airfoil for 
both a = 0° and 8° at low frequency is particularly interesting. The source of this discrepancy 
may be a combined effect of that discussed above and a transient laminar separation bubble. 
Evidence of a laminar separation bubble is suggested by oil flow visualization and, if present and 
coupled with the inflow would tend to mask the pressure response over the region of occurrence. 
Such a bubble would likely not be a highly repeatable phenomenon, especially if triggered by 
some very local leading edge feature. As such, it is feasible that this bubble was present to a 
greater extent during the acoustically treated test section measurement causing a low frequency 
reduction in the pressure response at the 1% chord location. That being said, the agreement 
between the two data sets is reasonable and therefore, they will be used interchangeably 
throughout the remaining analysis. 

4.4.3 Data Symmetry 

The next step in confirming the quality of measured unsteady pressure data is a check of 
the data symmetry. Unsteady pressure spectra on the pressure and suction side of the airfoil 
should be comparable within the uncertainty of the data at a = 0°. Similarly, data should be 
comparable from both sides of the airfoil measured at positive and negative angles of attack. 
That is, pressure spectra measured on the pressure side of the foil should agree at both positive 
and negative angles of attack and similarly for the suction side. 

Zero Angle of Attack: A comparison of spectra measured on both sides of the airfoil at a = 0° 
is presented in figure 4.20. In this figure the sides of the airfoil are designated as suction and 
pressure even though the airfoil is at zero angle of attack. The agreement atx/c = 1% is marginal 
with the suction side spectrum nearly 4dB above the pressure side for low frequencies. The 
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agreement improves with increasing frequency to near perfect at o\ = 100. Chordwise locations, 
x/c = 2.5, 4, and 6% show nearly perfect agreement over the measured frequency range. Some 
difference is observed at very high frequencies (op > 200) which may be related to resonance in 
the microphone mounting hole. The agreement at x/c = 9% is very good for 1 < co, <10. Below 
(Or < 1 the suction side spectrum rises ~4dB above the pressure side spectrum. Differences above 
(Or = 10 are likely the result of asymmetric boundary layer transition. The agreement returns to 
very near perfect at x/c = 14% with exception to some disagreement at co r > 200 due possibly to 
microphone mounting hole resonance. 

Positive and Negative Eight Degrees Angle of Attack: A comparison of pressure spectra at 
a = ±8° is presented in figure 4.21 for x/c = 1, 2.5, 4, and 6% chord. In these plots spectra at a = 
±8° on the pressure side (indicated by ‘Pressure’ in the figure) should be similar and likewise for 
the suction side (indicated by ‘Suction’). At 1% chord spectra measured on the pressure side of 
the airfoil fall nearly on top of each other. On the suction side at 1% the agreement between a = 
±8° p' spectra is poor with the a = -8° spectra falling 6dB below the a = 8° spectra. Spectra at 
2.5% chord on the suction side are found to be within 2dB of agreement. On the pressure side 
the agreement is slightly worse with 3dB of difference shown. By 4% chord data on both sides 
of the airfoil are within ldB over much of the low frequency range. For (0, > 10 asymmetric 
boundary layer transition likely produces the differences in pressure spectra observed at 4% 
chord. At 6% chord the agreement remains quite good at all frequencies with less than ldB of 
asymmetry occurring across the frequency range. 

Discussion on Sources of Asymmetries: Asymmetries are shown to occur primarily in the 
leading edge region ( x/c = 1 and 2.5%) for a = 0° and ±8°. Several asymmetry sources are 
investigated, including 

1 . Airfoil angle of attack 

2. Microphone calibrations 

3. Airfoil leading edge geometry 

Airfoil angle of attack 


64 




Chapter Four 


If the airfoil were at a very small angle of attack differences in the pressure spectra from 
one side of the airfoil to the other would be visible with the greatest difference showing up at the 
leading edge. However, the airfoil angle of attack was sighted to 0° and checked with mean 
pressure measurements which confirmed the angle of attack setting. Furthermore, since spectra 
at 4, 6, 9, and 14% chord are in near perfect agreement additional doubt is cast on this argument. 

Microphone calibrations 

Data used to form these plots are spanwise averaged where possible. That is, the auto- 
spectrum at each of the 8 spanwise locations is averaged to form the auto-spectrum for 
chordwise locations. This process assists in reducing the uncertainty associated with microphone 
calibrations. Additionally, microphone calibrations are consistent with each other at each 
chordwise location where spanwise averaging is not possible. 

Airfoil Geometry 

The airfoil thickness and chordwise location of microphones at y/c = 0.33 — up to 14% 
chord — was measured and is presented in figure 4.22 along with the theoretical NACA 0015 
geometry. Variation in chordwise location of microphones in any spanwise row relative to each 
other is less than 0.1mm. The measured chordwise locations are presented in table 4.4. 

Table 4.4: Measured chordwise location of microphones 


Nominal 

x/c 

Pressure side 
measured x/c 

Suction side 
measured x/c 

0.01 

0.0115 

0.0125 

0.025 

0.0255 

0.0270 

0.04 

0.0410 

0.0430 

0.06 

0.0610 

0.0630 

0.09 

0.0900 

0.0930 

0.14 

0.1400 

0.1440 


Microphones on the suction side tend to be located 0.1% to 0.3% chord further downstream than 
those on the pressure side. This result is unexpected based on the comparison of pressure spectra 
of adjacent microphones. That is, spectra decrease by approximately 1.25dB per 1% chord with 
increasing chordwise location. Therefore, if the suction side microphones are located slightly aft 
of the pressure side microphones then the measured spectra on the suction side of the airfoil 
should fall slightly below that on the pressure side of the airfoil. Just the opposite is shown to 
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occur in the spectra of figure 4.20 at low reduced frequencies. The measured thickness agrees 
well with the airfoils theoretical shape for the measured chordwise locations. 

There are a few other clues which may help explain the leading edge asymmetries that tie 
in with the leading edge geometry. The poor repeatability combined with the asymmetries at the 
most leading edge locations suggests that some sort of local flow phenomena is present. Oil 
flow visualization implies the presence of a laminar separation bubble on the suction side of the 
airfoil (for positive angles of attack). Such a bubble would be quite sensitive to local geometry 
and configuration. Additionally, if this bubble were coupled with the inflow the unsteady 
loading could be affected locally. As such, it is possible that this bubble develops at a different 
location and possibly over a greater or less extent of the chord in each measurement (solid wall 
and acoustic wall) and, for this reason, the repeatability is poor in the leading edge region. If this 
is indeed the case, such a local flow phenomenon could produce the asymmetries shown in 
pressure spectra at a = ±8° in both experiments. It is possible that the geometry at the leading 
edge is slightly different from one side to the other and this difference causes the laminar 
separation bubble to develop differently. If the laminar separation bubble on the suction side at 
a = -8° is larger compared to that which occurs on the suction side at a = 8°, and well coupled 
with the inflow then a significant reduction of the unsteady pressure could occur on the suction 
side at a = -8°. That is, the larger bubble, if coupled to the inflow, could be shedding vorticity 
and as such, reducing the loading on the airfoil and consequently driving down the unsteady 
pressure spectral level. 

No quantitative investigation of this theory through alternate flow visualization or other 
means was undertaken. However, a small side study was performed to test the effects of slight 
changes to the leading edge geometry. A small bump made of yam and covered with electrical 
tape (as discussed in chapter 2) was attached to the leading edge and pressure measurements 
made at a = -8°. The results of these measurements are presented in figure 4.23 which compares 
the p' spectra measured at a = -8° (with leading edge bump installed) with those measured at a = 
8° (with no leading edge bump present) for the 1% chord location. The agreement of pressure 
side data remains very good with less than ldB of difference occurring. The suction side data 
remains out of agreement; however, the a = -8° spectrum rises slightly so that the difference 
between the two is ~2.2dB. This is substantial improvement from the 6dB of difference which 
occurs with no leading edge bump present. The leading edge bump may be suppressing, to some 
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extent, the formation of the leading edge laminar separation bubble and in doing so increases the 
unsteady loading on the airfoil since less vorticity is shed by the smaller bubble. Therefore, this 
result does offer additional circumstantial support for the local flow phenomenon theory 
presented above. Furthermore, it suggests that viscous effects may play a significant role in 
determining the surface pressure response. The surface pressure is assumed to be a result of the 
inviscid response of the airfoil (i.e. the non-penetration condition is imposed at the airfoil 
surface) and therefore, local geometry should not affect the response of the airfoil to the inflow 
turbulence. This leading edge bump experiment, through its influence of the 1% chord pressure 
response, suggests that local geometry does play an important role and as such, viscosity is 
significant to determining the airfoil response. Therefore, it stands to reason that the 
asymmetries shown in figure 4.20 and 4.21 are related to a local flow phenomenon which is 
influenced by leading edge geometry. 

4.4.4 Contamination and Dominance of Various Pressure Fluctuation Sources 

This experiment sets out to measure the inviscid response of an airfoil in turbulence using 
microphones. The expectation here is that the microphones will measure pressure fluctuations 
associated with the non-penetration condition. However, measurements of unsteady pressure 
made with microphones are subject to a variety of contamination sources. Contamination can 
result from the wind tunnel acoustic field, boundary layer, stall, convection of eddies, and 
electrical noise. This section sets out to determine the extent of contamination from these 
sources and reveal the dominance of the inviscid response in unsteady pressure measurements. 

Acoustic Contamination: A somewhat qualitative estimate of the extent of acoustic field 

contamination is possible through comparison of smooth flow, matched RPM, and reference 
microphone measurements. The reference microphone is expected to provide an upper bound 
estimate of the acoustic field. Also, the reference microphone spectrum falls near 30dB below 
surface pressure spectra for leading edge chordwise locations (1% < x/c < 14%). As such, the 
extent of acoustic contamination can be estimated through a comparison of reference microphone 
data with that measured on the surface of the airfoil in smooth flow and matched RPM. 

Auto-spectra of the across chord pressure difference are presented at eight chordwise 
locations (1, 2.5, 4, 6, 9, 14, 20, and 30%) for the cases of smooth flow at all angles of attack, 
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matched RPM at a = 0°, and reference microphone measurement in figures 4.24 and 4.25. 
Figures 4.26 and 4.27 present the corresponding coherence and phase of these signals. At the 
unstalled (0°, 4°, 8°, 12°) angles of attack the smooth flow and matched RPM spectra remain as 
much as 30dB below the reference microphone spectrum at all eight chordwise locations over 
much of the frequency range, varying from as high as 1 15dB at a = 12° to as low as 65dB at a = 
0°. The reference microphone, as discussed in chapter 2 may, by its existence in the flow, 
generate pressure fluctuations and thus the spectrum measured represents an upper bound on the 
tunnel acoustic field. Therefore, the fact that smooth flow data remains below the upper bound 
tunnel acoustic field indicates the low level of noise present at the airfoil surface. Also, across 
chord phasing (figures 4.26 and 4.27) remains at nearly 0° at each chordwise location for the 
unstalled angles of attack indicating the primary source of pressure fluctuation is the ambient 
acoustic field. Coherence, although reaching to 0.8 at some frequencies is quite ‘spiky’ and falls 
off at higher reduced frequencies (>10) beyond 9% chord demonstrating the weak nature of the 
dominant pressure fluctuation. 

Stalled angle of attack spectra (16° and 20°) are well above the reference microphone 
spectrum, which is expected considering the unsteadiness associated with stall. The a = 16° and 
20° spectra remain within 8dB of each other at all eight chordwise locations with, interestingly 
the a = 16° spectrum in general residing at a higher level than the a = 20° spectrum. These 
spectra achieve between 120 and 130dB at the lowest reduced frequency and fall to as low as 
80dB at the highest reduced frequency. Therefore, it can be concluded that at stalled angles of 
attack the large unsteady pressure fluctuations occurring in the presences of the separated flow 
will dominate the unsteady pressure measurement. 

Therefore, based on the low smooth flow, matched RPM, and reference microphone 
spectral levels, ambient acoustic related pressure fluctuations are not expected to significantly 
contaminate the data set. However, at stalled angles of attack, the large pressure fluctuations 
associated with separated flow are expected to dominate the unsteady pressure measurement. 

Turbulent Boundary Layer Contamination: To aid in assessing the extent of contamination 
from turbulent boundary layer pressure fluctuations a technique presented by Burley et al. (1997) 
is used. Burley et al. (1997) developed a method for extracting the desired pressure 
fluctuation — inviscid response pressure fluctuations — from microphone measurements. They 
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assume that the inviscid response will correlate across the chord of the airfoil. As such, the 
pressure fluctuations related to the inviscid response can be extracted by correlating with the 
airfoil leading edge (or at least the most forward pressure measurement location). The leading 
edge is used because it is assumed to be dominated by pressure fluctuations related to the 
inviscid response. In this way, the pressure fluctuations related to the turbulent boundary layer 
can be removed (since the boundary layer pressure fluctuations will not correlate of large 
distances) leaving the inviscid response auto-spectra. This method employs what is referred to 
as the Coherent Output Power (COP) and extracts what amounts to the inviscid response surface 
pressure auto-spectra. The COP is defined as 


J AP y AP y Inviscid 


J APAP r 


J APAP 


(4.5) 


where Gapiapi is the pressure difference auto-spectrum at 1% chord and Gapiapx is the pressure 
difference cross-spectra between microphones located at 1% chord and position x. 

Figure 4.28 presents large grid pressure difference auto-spectra measured at 1, 2.5, 4, 6, 
9, and 14% chord plotted with the inviscid response auto-spectra (COP — presented as Burley- 
Brooks in plots) at a = 0°. The measured auto-spectra compares well with COP at all chordwise 
location for co r < 10 indicating that significant acoustic and boundary layer pressure fluctuations 
are not present. Above reduced frequencies of 10 COP begins to fall below measure auto-spectra 
signifying the onset of significant acoustic and/or boundary layer pressure fluctuations. Due to 
the comparatively low acoustic level in the tunnel the measured auto-spectra are believed to be 
primarily elevated as a consequence of small scale (higher frequency) boundary layer pressure 
fluctuations. 

Now, consider the COP calculated for small grid unsteady pressure measurements. Auto- 
spectra measured in small grid turbulence are expected to contain some boundary layer pressure 
fluctuation contamination as a consequence of grid generated turbulence scale. This turbulence 
scale is on the same order as the boundary layer scale and thus the grid generated turbulence 
tends to add energy to the boundary layer. Additionally, the hydrodynamic pressure fluctuations 
associated with the convection of turbulence down the chord may contribute significantly to the 
overall pressure response. Figure 4.29 presents a comparison of small grid measured pressure 
difference auto-spectra with inviscid response auto-spectra (COP) at 1, 2.5, 4, 6, 9, and 14% 
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chord. Inviscid response auto-spectra (COP) fall below measured auto-spectra by as much as 1 
order of magnitude in spectral level at 2.5% chord and continue to fall relative to measured auto- 
spectra at the remaining chordwise locations. The large difference between measured and 
inviscid response auto-spectra results from the lack of correlation of pressure fluctuations with 
the leading edge. The inviscid response of the airfoil is, of course, still present; however, 
pressure fluctuations associated with the turbulent boundary layer and the convection of the free- 
stream turbulence along the airfoil chord are significant enough to mask the inviscid response 
and cause the COP to fall considerably relative to measured auto-spectral levels. 


Measurement System Noise: Measurement system noise is another source of contamination of 
particular concern since it will correlate across all microphones thereby, artificially driving up 
cross-spectral levels. The effects of measurement system noise can be assessed through 
comparison of large and small grid pressure spectra with pressure spectra measured in smooth 
flow. Surface pressure measurements made in smooth flow form a baseline measure of noise (as 
shown above), including measurement system noise. Comparing cross-spectra from the two 
measurements (airfoil in turbulence with smooth flow) will reveal which is more dominant. If 
cross-spectra measured with the airfoil in turbulence have a higher spectral level then they are 
not expected to be contaminated with measurement systems noise. Pressure difference cross- 
spectra are formed as 


GaP,APj — E ([^E/ { X i ) ( X i )] Pu ( X ; ) ( X j )_ 


(4.6) 


Figure 4.30 through 4.32 presents a comparison of pressure difference cross-spectra 
measured in large and small grid turbulence and smooth flow pressure difference cross-spectra 
for chordwise locations (denoted by index i of above equation) of 1%, 6%, and 14% at a = 0°. 
The cross-spectra at each location i are taken with 6 other chordwise locations (denoted by index 
j in the above equation) at 1%, 2.5%, 4%, 6%, 9%, and 14%. This data has been presented in 
pressure difference form. 

In the case of large grid turbulence, surface pressure cross-spectra are more than two 
decades in spectral level above cross-spectra measured with no grid turbulence and show 
distinctly different shape. These characteristics are visible at each chordwise location. Now, 
considering the small grid turbulence surface cross-spectra, it is shown that these spectral levels 
are as much as a decade in magnitude above smooth flow cross-spectra. These spectra also take 
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on a distinctive shape and this characteristic is present at each chordwise location. The 
difference in the shape and the increase in spectral level suggest the presence of correlated 
pressure fluctuations which are significantly more dominant than those measured in smooth 
flow. As such, this data does not appear to be contaminated by measurement system noise, 
boundary layer noise, or ambient acoustic noise. 


Dominance of Inviscid Response and Hydrodynamic Pressure Fluctuations: The correlation 
of pressure fluctuations in space and time offers an effective technique for identifying dominant 
types of pressure fluctuations occurring on the airfoil surface. Such correlations are referred to 
as pressure-pressure space-time correlation, or cross-correlation R PP {x,x', z). Investigating the 


pressure-pressure correlations will help assess the extent of viscous phenomena (or 
hydrodynamic pressure fluctuations) involved in the total airfoil response. 

The pressure-pressure correlation, R PP ■ is defined as 




(4.7a) 


where p 2 is the mean-square pressure. The quantity can also be computed from the cross-power 
spectral density of surface pressure fluctuations as 
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(4.7b) 


The space-time correlation is a three dimensional function (when 7 = 0) which, when contour 
plotted, can reveal the speeds that pressure fluctuations are moving along the airfoil chord. For 
example, consider the contour plot of R PP ’ shown in figure 4.33 with space on the vertical axis 


(x/c) and time on the horizontal axis ( zUJc ). Ridges in the correlation will have a slope 
associated with them that corresponds to a velocity. This is the velocity at which pressure 
fluctuations are moving along the airfoil chord. As such, pressure fluctuations that convect with 
the mean flow (such as eddies) have a slope of 1 (x/z = U x ) while propagating (i.e. at the speed 
of sound) pressure fluctuations will have slope of 12(sound speed + mean free-stream, x/z = 
12U„). 
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Figure 4.34 presents the correlation of pressure with the suction side 1% chord 
measurement location in space and time at a = 0°, 2°, 4°, 6°, 8°, 10°, 12°, 14°, 16°, 18°, 20° for 
the airfoil in large grid turbulence. On these plots, positive x/c indicates the pressure side of the 
airfoil and negative x/c the suction side. The measured correlation is formed from pressure data 
sampled at 25.6 kHz and is high pass filtered at 100Hz to remove uncertainties associated with 
microphone calibrations. Both sides of the airfoil are correlated with the same microphone (1% 
chord, on the suction side). At all angles of attack the correlations shows a ridge with a 12:1 
slope occurring on the pressure side of the airfoil. This region of high correlation is the airfoils 
inviscid response and clearly dominates the unsteady pressure measurement at all angles of 
attack. All angles of attack (more so with increasing angle of attack) show a ridge in the 
correlation with a slope of 1 : 1 on the airfoils suction side which likely results from convection of 
turbulence by the free-stream. The suction side also shows a strong inviscid response with a 
ridge in the correlation occurring with a 12:1 slope. It is interesting to see the break in the 
correlation function on the suction side which is evident at all angles of attack and occurs around 
the 15% chord location. As the angle of attack is increased this region of no correlation tends to 
grow along the chord of the airfoil covering as much as 20% of the chord at the higher angles of 
attack. This region may be related to the formation of a laminar separation bubble (as suggested 
by flow visualization) which masks the response of microphones in this area from the response 
of the entire airfoil. When the airfoil stalls at a=20° disturbances begin to convect at UJ2 which 
is the characteristic convection speed in a separated flow region. In summary, figure 4.34 clearly 
demonstrates the dominance of the airfoils inviscid response when interacting with large grid 
turbulence. 

For completeness sake, the correlation of pressure with the 1% chord measurement 
location on the pressure side of the airfoil in space and time is presented in figure 4.35 for large 
grid flow at a = 0°, 2°, 4°, 6°, 8°, 10°, 12°, 14°, 16°, 18°, 20°. The cross-correlation shows 
strong correlation of pressure on both sides of the airfoil along the chord for a < 8° that appears 
to be related to the airfoils inviscid response. The correlation function also has a ridge with a 
slope near 12:1 (i.e. propagation of pressure waves). A break in the correlation function is also 
observed to occur near 15% chord and convected pressure fluctuations appear to be emerging 
from this area for a < 8°. Signs of eddy convection on the suction side of the airfoil (recalling 
the correlation function is formed with the 1% chord microphone on the pressure side) are 


72 




Chapter Four 


surprising and suggest the presence of a transient laminar separation bubble that is coupled with 
the inflow. Such coupling would alter the unsteady loading and therefore, result in correlation 
between convected disturbances on the suction side of the airfoil with measurements on the 
pressure side of the airfoil. For 8°, 10°, 12°, and 14° angle of attack, significant correlation 
related to propagating pressure waves is observed along the chord on the pressure side. Eddy 
convection appears to be present on the pressure side as well, indicated by the 1:1 slope of the 
correlation ridge. Convected pressure fluctuations with a slope less than 1 : 1 also become more 
significant on the suction side; however, significantly correlated propagated pressure fluctuations 
appear to vanish. Beyond a = 14° propagated pressure fluctuations dominate the pressure side 
with significant convection occurring on the suction side. The convection on the suction side 
appears to occur at a speed close to have the mean flow velocity based on the slope of the 
correlation ridge. 

Additional support for the dominance of the inviscid response is provided by the pressure 
correlation function formed with the 1% chord measurement location of the pressure side of the 
airfoil (figure 4.35). Furthermore, this analysis has helped to substantiate the possible presence 
of a transient laminar separation bubble coupled with the inflow and occurring on the suction 
side of the airfoil. 

Figure 4.36 presents the pressure-pressure space-time correlation from small grid 
turbulence/airfoil interaction surface pressure measurements at a = 0°, 4°, 8°, 12°, 16°, and 20°. 
The most notable occurrence in figure 4.36 is the reduction in correlation level relative to the 
large grid data. The inviscid response is present up to a = 12° on both sides of the airfoil 
indicated by the correlation ridge occurring along the airfoil chord at near a 12:1 slope. At all 
angles of attack pressure fluctuations convecting along the airfoil chord (1:1 slope) are also 
visible. A standing wave pattern also is visible in the correlation at a = 4° and 8° which may be 
related to system noise. The wave has a period close to that which would corresponds to 60Hz 
power system noise. The pressure-pressure correlation seems to indicate that the inviscid 
response is not particularly strong in the small grid flow. The inflow length scale being l/80 th of 
the airfoil chord may explain why this turbulence/airfoil interaction does not excite a strong 
inviscid response. With little evidence of a significant inviscid response in these pressure 
measurements the aeroacoustic relevance of this turbulence/airfoil interaction is minimized. 
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That being said, the remaining portion of this dissertation will focus on analysis of large grid 
pressure fluctuations in which the inviscid response is significant and dominant. 

Lastly, the theoretical form of the cross-correlation function, as predicted by Amiet’s 
(1976a, b) unsteady loading theory is presented and compared with the large grid measured 
cross-correlation at a = 0° in figure 4.37. The predicted correlation drops off slightly moving 
away from the leading edge; however, the slope of the peak correlation is infinite indicating that 
the entire airfoil responses instantaneously to the inflow upwash. This, of course is much 
different from the real interaction in which the airfoil pressure field must adjust at the speed of 
sound (hence the 12:1 slope of the correlation peak). Due to hydrodynamic pressure fluctuations 
the measured correlation drops off faster moving away from the leading edge when compared 
with the predicted. The prediction also does not account for the convection of eddies along the 
airfoil chord which occurs in the real flow. Consequently, the predicted correlation does not 
show a peak occurring with a 1 to 1 slope which would be indicative of eddy convection. In 
summary, the predicted pressure-pressure correlation captures the inviscid response of the airfoil 
well in the leading edge region; however, the prediction fails to drop off as quickly moving along 
the airfoil chord as shown to occur in the measurement. 

4.5 Summary of Experimental Details and Data Quality 

Unsteady pressure measurements on the surface of an airfoil in turbulence were made 
with the goal of measuring pressure fluctuations associated with the airfoils inviscid response. 
The airfoil was immersed in two grid generated flows to assess the effect of turbulence scale on 
the airfoil response. Measurements were made at angles of attack ranging from a = 0° to 20°. 

Two separate tests were conducted in large grid flow; one measuring unsteady pressure 
with solid steel test section walls, the other with acoustically treated walls present. Comparison 
of data from these two tests reveals poor repeatability in the leading edge region. Furthermore, 
asymmetries are shown in the leading edge region. It is believed repeatability and asymmetry 
issues arise from a local flow phenomenon, such as a transient laminar separation bubble. Oil 
flow visualization and investigation of the pressure-pressure space-time correlation supports the 
presence of such bubble. The repeatability and asymmetry issues appear to be representing real 
events on the airfoil surface which impact the airfoils inviscid response. As such, the 
repeatability and asymmetry issues are not looked at as corrupting features, but rather as 
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additional intimations which can aid in understanding the physics of turbulence-airfoil 
interaction. 

Acoustic, boundary layer, and electrical noise contamination are investigated. Large grid 
unsteady pressure data is shown to be largely contamination free over much of the frequency 
range. Small grid data is shown to contain substantial hydrodynamic pressure fluctuations 
although, is largely free of acoustic and electrical noise contamination. Additional support is 
given to this conclusion based on analysis of the pressure-pressure space-time correlation. The 
space-time correlation shows the inviscid response is dominant in unsteady pressure 
measurements when the turbulence is large; however, the inviscid response is weak and pressure 
measurements contain significant hydrodynamic pressure fluctuations when the turbulence is 
small. Based on these results the remaining analysis of this data will focus on large grid 
turbulence-airfoil interaction. 
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Figure 4.1: Top down view of airfoil in wind tunnel test section defining 

positive angles of attack 



Figure 4.2: Single airfoil in turbulence 
coordinate system 
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Figure 4.3: Mean pressure calculated lift coefficient shown with McKeough’s (1976) data 
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Figure 4.4: Stagnation and separation locations of 
small scale turbulence flow 



x/c 

Figure 4.5: Stagnation, transition, and separation locations of large 

scale turbulence flow 
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Figure 4.6: a = 0° oil flow visualization on the hatch-pressure-(a) and so\\d-suction-(b) side of 
airfoil. The numbers in black correspond to the location of microphones in percent chord 
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Figure 4.7: a = 4° oil flow visualization on the hatch-pressure-(a) and so\\d-suction-(b) side of 
airfoil. The numbers in black correspond to the location of microphones in percent chord 
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Figure 4.8: a = 8° oil flow visualization on the hatch-pressure-(a) and so\\d-suction-(b) side of 
airfoil. The numbers in black correspond to the location of microphones in percent chord 
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Figure 4.9: a = 12° oil flow visualization on the hatch-pressure-(a) and so\id-suction-(b) side of 
airfoil. The numbers in black correspond to the location of microphones in percent chord 
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Figure 4.10: a = 16° oil flow visualization on the hatch-pressure-(a) and so\\d-suction-(b) side of 
airfoil. The numbers in black correspond to the location of microphones in percent chord 
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Figure 4.11: a = 20° oil flow visualization on the hatch-/>rewwre-(a) and so\\d-suction-(b) side of 
airfoil. The numbers in black correspond to the location of microphones in percent chord 
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Suction side--x/c=0.01 



Figure 4.12: Spanwise homogeneity of small grid unsteady pressure data at a = 0° on the suction 

side of airfoil for chordwise positions x/c = 1, 2.5, 4, 6, 9%. First colu mn of numbers in legend 

indicates spanwise position and second column is unique number given to each microphone. 
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Figure 4.12 continued: Spanwise homogeneity of small grid unsteady pressure data at a = 0° on the 
suction side of airfoil for chordwise positions x/c = 1, 2.5, 4, 6, 9%. First column of numbers in 
legend indicates spanwise position and second column is unique number given to each microphone. 
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Figure 4.12 continued: Spanwise homogeneity of small grid unsteady pressure data at a = 0° on the 
suction side of airfoil for chordwise positions x/c = 1, 2.5, 4, 6, 9%. First column of numbers in 
legend indicates spanwise position and second column is unique number given to each microphone. 
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Figure 4.13: Spanwise homogeneity of small grid unsteady pressure data at a = 0° on the pressure 
side of airfoil for chordwise positions x/c = 1, 2.5, 4, 6, 9%. First column of numbers in legend 
indicates spanwise position and second column is unique number given to each microphone. 
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Figure 4.13 continued: Spanwise homogeneity of small grid unsteady pressure data at a = 0° on the 
pressure side of airfoil for chordwise positions x/c = 1, 2.5, 4, 6, 9%. First column of numbers in 
legend indicates spanwise position and second column is unique number given to each microphone. 
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Figure 4.13 continued: Spanwise homogeneity of small grid unsteady pressure data at a = 0° on the 
pressure side of airfoil for chordwise positions x/c = 1, 2.5, 4, 6, 9%. First column of numbers in 
legend indicates spanwise position and second column is unique number given to each microphone. 
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co r 


CO 


r 


Figure 4.14: Spanwise homogeneity of large grid unsteady pressure data at a = 0° on the suction 
side of airfoil for chordwise positions x/c = 1, 2.5, 4, 6, 9%. First column of numbers in legend 
indicates spanwise position and second column is unique number given to each microphone. 
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Suction side--x/c=0.06 



Figure 4.14 continued: Spanwise homogeneity of large grid unsteady pressure data at a = 0° on the 
suction side of airfoil for chordwise positions x/c = 1, 2.5, 4, 6, 9%. First column of numbers in 
legend indicates spanwise position and second column is unique number given to each microphone. 
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Figure 4.14 continued: Spanwise homogeneity of large grid unsteady pressure data at a = 0° on the 
suction side of airfoil for chordwise positions x/c = 1, 2.5, 4, 6, 9%. First column of numbers in 
legend indicates spanwise position and second column is unique number given to each microphone. 
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Pressure side--x/c=0.025 
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Figure 4.15: Spanwise homogeneity of large grid unsteady pressure data at a = 0° on the pressure 
side of airfoil for chordwise positions x/c = 1, 2.5, 4, 6, 9%. First column of numbers in legend 
indicates spanwise position and second column is unique number given to each microphone. 
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Figure 4.15 continued: Spanwise homogeneity of large grid unsteady pressure data at a = 0° on the 
pressure side of airfoil for chordwise positions x/c = 1, 2.5, 4, 6, 9%. First column of numbers in 
legend indicates spanwise position and second column is unique number given to each microphone. 
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Figure 4.15 continued: Spanwise homogeneity of large grid unsteady pressure data at a = 0° on the 
pressure side of airfoil for chordwise positions x/c = 1, 2.5, 4, 6, 9%. First column of numbers in 
legend indicates spanwise position and second column is unique number given to each microphone. 
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Figure 4.16: Comparison of large grid p' spectra measured with solid and acoustically treated 
wind tunnel test section walls. Spectra are from suction side of airfoil (x/c = 1, 2.5, 4, 6, 9, 
14%) at a = 0° and are normalized on the free-stream dynamic pressure. 
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Figure 4.16 continued: Comparison of large grid p' spectra measured with solid and 
acoustically treated wind tunnel test section walls. Spectra are from suction side of airfoil 
(x/c = 1, 2.5, 4, 6, 9, 14%) at a = 0° and are normalized on the free-stream dynamic pressure. 
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Figure 4.16 continued: Comparison of large grid p' spectra measured with solid and 
acoustically treated wind tunnel test section walls. Spectra are from suction side of airfoil 
(x/c = 1, 2.5, 4, 6, 9, 14%) at a = 0° and are normalized on the free-stream dynamic pressure. 


99 


(coVq 


Pressure side--x/c=0.01 Pressure side--x/c=0.025 



Figure 4.17: Comparison of large grid p' spectra measured with solid and acoustically treated 
wind tunnel test section walls. Spectra are from pressure side of airfoil (x/c = 1, 2.5, 4, 6, 9, 
14%) at a = 0° and are normalized on the free-stream dynamic pressure. 
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Figure 4.17 continued: Comparison of large grid p' spectra measured with solid and 
acoustically treated wind tunnel test section walls. Spectra are from pressure side of airfoil 
( x/c = 1, 2.5, 4, 6, 9, 14%) at a = 0° and are normalized on the free-stream dynamic pressure. 
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Figure 4.17 continued: Comparison of large grid p' spectra measured with solid and 
acoustically treated wind tunnel test section walls. Spectra are from pressure side of airfoil 
( x/c = 1, 2.5, 4, 6, 9, 14%) at a = 0° and are normalized on the free-stream dynamic pressure. 


102 


b/(®) 



Figure 4.18: Comparison of large grid p' spectra measured with solid and acoustically treated 
wind tunnel test section walls. Spectra are from suction side of airfoil ( x/c = 1, 2.5%) at a = 
8° and are normalized on the free-stream dynamic pressure. 
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Figure 4.19: Comparison of large grid p ' spectra measured with solid and acoustically treated 
wind tunnel test section walls. Spectra are from pressure side of airfoil (x/c = 1, 2.5%) at a = 
8° and are normalized on the free-stream dynamic pressure. 
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Figure 4.20: Comparison of p' spectra on the pressure and suction side of airfoil at a 
Comparison made at x/c = 1, 2.5, 4, 6, 9, 14%. 
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Figure 4.20 continued: Comparison of p' spectra on the pressure and suction side of airfoil 
at a = 0°. Comparison made at x/c = 1, 2.5, 4, 6, 9, 14%. 
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Figure 4.20 continued: Comparison of p' spectra on the pressure and suction side of airfoil 
at a = 0°. Comparison made at x/c = 1, 2.5, 4, 6, 9, 14%. 
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Figure 4.21: Comparison of pressure and suction side p‘ spectra at a = 8° and -8° for x/c = 

0.01,0.025,0.04, and 0.06 
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Figure 4.21 continued: Comparison of pressure and suction side p‘ spectra at a = 8° and -8° 

for x/c = 0.01, 0.025, 0.04, and 0.06 
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Figure 4.22: Measured geometry at each microphone location 
versus actual wing geometry 
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Figure 4.23: Pressure spectra comparison at a = 8° and -8° for x/c 
0.01 . Data at a = -8° is taken with leading edge bump installed. 
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Figure 4.24: Pressure difference smooth flow data plotted at all angles of attack with matched 
RPM data at a = 0° and reference (B&K) microphone data at x/c = 0.01, 0.025, 0.04, 0.06 
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Figure 4.25: Pressure difference matched flow speed data plotted at all angles of attack with matched 
RPM data at a = 0° and reference (B&K) microphone data at x/c = 0.09, 0.14, 0.2, 0.3 
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Figure 4.26: Coherence and phase of matched flow speed data plotted at all angles of attack 
with matched RPM data at a = 0° and reference (B&K) microphone data at x/c = 0.01, 0.025, 

0.04, 0.06 
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Figure 4.27: Coherence and phase of matched flow speed data plotted at all angles of attack with 
matched RPM data at a = 0° and reference (B&K) microphone data at x/c = 0.09, 0.14, 0.2, 0.3 
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Figure 4.28: Large grid auto-spectra plotted with COP (Burley-Brooks) at a = 0° for x/c 

0.025, 0.04, 0.06, 0.09, and 0. 14 
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Figure 4. 29: Small grid auto- 
spectra plotted with COP (Burley- 
Brooks) at a = 0° for x/c = 0.025, 
0.04, 0.06, 0.09, and 0.14 
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Figure 4.30: Pressure difference cross-spectra at x/c = 0.01, a = 0° measured in 
smooth flow and in large and small grid turbulence 
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Figure 4.31: Pressure difference cross-spectra at x/c = 0.06, a = 0° measured in 
smooth flow and in large and small grid turbulence 
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Figure 4.32: Pressure difference cross-spectra at x/c = 0.14, a = 0° measured in 
smooth flow and in large and small grid turbulence 
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Figure 4.34: Contour plots of pressure-pressure correlation, 
R pp (x=\%,x r) for large grid flow. The correlation is formed with 
the 1 % chord microphone on the suction side of the airfoil 
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Figure 4.34 continued: Contour plots of pressure-pressure 
correlation, R pp (x=\ %,x r) for large grid flow. The correlation is 
formed with the 1% chord microphone on the suction side of the airfoil 
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Figure 4.34 continued: Contour plots of pressure-pressure 
correlation, R pp (x=\ %,x r) for large grid flow. The correlation is 
formed with the 1% chord microphone on the suction side of the airfoil 
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Figure 4.34 continued: Contour plots of pressure-pressure 
correlation, R pp (x=\ %,x r) for large grid flow. The correlation is 
formed with the 1% chord microphone on the suction side of the airfoil 
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Figure 4.34 continued: Contour plots of pressure-pressure 
correlation, R pp (x=\ %,x r) for large grid flow. The correlation is 
formed with the 1% chord microphone on the suction side of the airfoil 
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Figure 4.34 continued: Contour plots of pressure-pressure 
correlation, R pp (x=\ %,x r) for large grid flow. The correlation is 
formed with the 1% chord microphone on the suction side of the airfoil 
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Figure 4.35: Contour plots of pressure-pressure correlation, 
R pp (x=\%jc t) for large grid flow. The correlation is formed with the 
1% chord microphone on the pressure side of the airfoil 
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Figure 4.35 continued: Contour plots of pressure-pressure correlation, 
R pp (x=l% r K f) for large grid flow. The correlation is formed with the 1% 
chord microphone on the pressure side of the airfoil 
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Figure 4.35 continued: Contour plots of pressure-pressure correlation, 
R pp (x=l% r K f) for large grid flow. The correlation is formed with the 1% 
chord microphone on the pressure side of the airfoil 
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Figure 4.35 continued: Contour plots of pressure-pressure correlation, 
R pp (x=l% r K f) for large grid flow. The correlation is formed with the 1% 
chord microphone on the pressure side of the airfoil 
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Figure 4.35 continued: Contour plots of pressure-pressure correlation, 
R pp (x=l% r K f) for large grid flow. The correlation is formed with the 1% 
chord microphone on the pressure side of the airfoil 
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Figure 4.35 continued: Contour plots of pressure-pressure correlation, 
R pp (x=l% r K f) for large grid flow. The correlation is formed with the 1% 
chord microphone on the pressure side of the airfoil 
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Figure 4.36: Contour plots of R pp (x=l%,x r) for small grid flow. 
The correlation is formed with the 1% chord microphone on the 
suction side of the airfoil 
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Figure 4.36 continued: Contour plots of R pp (x=\%,x r) for small 
grid flow. The correlation is formed with the 1% chord microphone 
on the suction side of the airfoil 


135 




faeoe-B o o" o 


Amiet Prediction 



Large Grid , a = 0° 



tU /c 


Figure 4.37: Comparison of measured and predicted pressure-pressure 
correlation, R pv {x=l%jc r) for large grid flow. The correlation is formed 
with the 1% chord microphone on the suction side of the airfoil. 
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5. Single Airfoil in Turbulence: Analysis of Unsteady Pressure Data 

This chapter presents analysis of unsteady surface pressure measured on the NACA 0015 
airfoil in the large grid turbulence. Measurements of unsteady pressure were made at angles of 
attack, a = 0° to 20°. Comparisons between Amiet’s (1976a, b) flat plate, zero angle of attack 
theory and measured pressure spectra are presented. Effects of angle of attack on pressure 
spectra are examined and comparisons made with Reba and Kerschen’s (1996) unsteady loading 
theory. Unsteady pressure measurements are then used to calculate lift and the effects of angle 
of attack examined. Measured lift spectra are compared with Amiet’s (1976a, b) theory for a = 
0° and Reba and Kerschen’s (1996) at non-zero angle of attack. The sources of observed angle 
of attack effect are examined through investigation of pressure correlation length scales, a 
modified trailing edge, and scaling of unsteady lift. Based on this analysis, a prediction scheme 
is presented which accurately captures the effects of angle of attack at low co r within the 
parameters of this experiment. Lastly, this chapter demonstrates how surface pressure 
measurements can be used to make estimates of far-field noise. 

5.1 Pressure Spectra and Mean Square Pressure 

5.1.1 Pressure Difference Spectra at Zero Angle of Attack 

Measured pressure spectra are compared with Amiet’s (1976a, b) theory in this section. 
Recall from chapter 3, Amiet (1976a, b) theoretically computes the unsteady response of a flat 
plate airfoil at zero angle of attack. Figure 5.1 presents a comparison of Amiet’s predicted 
pressure difference spectra (see appendix A for definition of pressure difference) with measured 
pressure difference spectra at a = 0° for chordwise locations, x/c = 1, 2.5, 4, 6, 9, and 14% 
normalized on free-stream dynamic pressure, q. Both the measurement and prediction fall 
monotonically moving away from the leading edge. At all chordwise locations Amiet’s theory 
agrees well with the measured shape of Gapap spectra. However, the theory tends to over 
predict spectra levels, particularly at 1% chord with as much as 7dB of difference occurring here. 
The over prediction of surface pressure in the leading edge region is likely the consequence of 
the leading edge singularity which occurs in Amiet’s theory. The agreement in spectral level 
improves to within ldB moving away from the leading edge. The measured data tends to roll off 
slightly faster than the prediction over 5 < op < 20. The disagreement in slope may be a result of 
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the airfoil thickness. Other researchers have shown similar results, although unpublished 
(private communication, Gershfeld 2002). The slope of the measured data decreases for co r > 20. 
The change is subtle for the first 4 chordwise locations and may be related to convection of 
eddies by the free-stream. At x/c = 9 and 14% the slope changes substantially more for co r > 20 
as significant boundary layer pressure fluctuations begin to affect the high frequency response. 

At some chordwise locations a high frequency spike is present in measured Gapap- This 
spike, occurring over a range of 300 < co r < 800 appears to be related to resonance within the 
microphone mounting hole that was not accounted for in the microphone calibration. 

5.1.2 Pressure Difference Spectra at Non-Zero Angle of Attack 

Measured pressure spectra are compared with Reba and Kerschen’s (1996) linearized 
RDT predictions for similar chordwise locations and angles of attack. Such a comparison will 
aid in assessing the effects of turbulence distortion on surface pressures. Recall from chapter 3, 
Reba and Kerschen (1996) developed a solution for the response of a flat plate at small angle of 
attack using Rapid Distortion Theory (RDT). The FORTRAN code produced by Reba and 
Kerschen (1996) is used to calculate the airfoil response from which the pressure difference 
cross-spectral density is computed (equation 3.16). However, due to the restriction in Reba and 
Kerschen’s theory of low acoustic wavenumber it is not possible to integrate over a large range 
of spanwise wavenumbers. Consequently, the absolute level of predictions from Reba and 
Kerschen’s theory are somewhat questionable. As such, a direct comparison between predicted 
and measured unsteady pressure level is not appropriate, although a comparison of spectral shape 
and variation with angle of attack and chordwise position seems reasonable. 

Figure 5.2 presents the pressure difference spectra as a function of frequency at six 
chordwise ( x/c = 1, 2.5, 4, 6, 9, 14%) locations and three angles of attack (a = 0°, 2°, 6°) 
normalized on free-stream dynamic pressure, q. This figure presents the RDT predictions of 
Reba and Kerschen (1996) and the measured pressure difference spectra. The predicted spectra 
of Reba and Kerschen fall monotonically moving away from the leading edge at all angles of 
attack, as shown to occur in measured data. Additionally, the spectra show a somewhat flat 
range for co r < 2 and then begin to fall off with a slope similar to that shown in measured data. 
The agreement in level between the predicted a = 0° spectrum and the measured a = 0° spectrum 
at 1% chord is believed to be a coincidence for the reasons stated above. 
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The interesting occurrence shown in figure 5.2 is the significant rise in predicted spectral 
level with increasing angle of attack. This rise is particularly dominant at low reduced 
frequencies, with up to a 16dB increase and occurs at all chordwise locations. The spectral 
levels tend to converge at higher reduced frequencies. The measured data shows a slight 
increase (2.5dB) at 1% chord; however tends to show little variation as the angle of attack is 
increased. The reason for the large change in spectral level predicted by Reba and Kerschen 
(1996) with increasing angle of attack is likely related to the flat plate model they assume. The 
potential field associated with a flat plate at angle of attack will distort the inflow substantially 
more than an actual airfoil with thickness due to the very large velocity gradients that stretch and 
tilt the vorticity vectors of the inflow turbulence in the vicinity of the leading edge singularity. 
Clearly, though the main result shown here is that a flat plate RDT representation of turbulence- 
airfoil interaction is not sufficient for correctly predicting the effects of angle of attack on the 
surface pressure given the relative inflow scale and airfoil chord used in this experiment. 

5.1.3 Mean Square Pressure at Zero Angle of Attack 

It is now interesting to examine Amiet’s prediction of the root mean square (RMS) 
pressure and its agreement with measured RMS pressures. The unsteady pressure spectrum at 
each chordwise location on the pressure and suction side of the airfoil is integrated to obtain root 
mean square p' values as 

/max 

p 2 (x)= | S qq (x,x',r 1 ,f)df (5.1) 

/min 

where S qq is the surface pressure cross-spectral density, x = x\ and r/ = 0. The integration is 
performed from a minimum frequency, f min = 6Hz (co r = 0.4) to a maximum frequency, f max = 
1570Hz (co r = 100). These values are compared with those calculated from Amiet’s (1976a, b) 
theory in figure 5.3 at a = 0° as a function of chordwise position (up to 14% chord) and are 
normalized on the ffee-stream dynamic pressure, q. The pressure spectra presented in figure 5.1 
show very good agreement with Amiet’s theory at chordwise locations beyond 2.5%. Therefore, 
not surprisingly, the agreement of p' values compare within 3% across the chord. The 1% chord 
location shows a drop in measured p' on the pressure side which may be related to a local flow 
phenomenon as discussed in chapter 4. In the leading edge region Amiet’s prediction begins to 
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rise rapidly as a consequence of the leading edge singularity in this theory, producing large 
differences between measured and predicted RMS pressure. 

5.1.4 Mean Square Pressure at non-Zero Angle of Attack 

The effects of angle of attack on p' values are investigated in figures 5.4 and 5.5. The 
spectrum from chordwise locations 1% to 14% are each integrated to obtain p 2 values for both 
pressure and suction sides of the airfoil at all angles of attack. The integration is performed over 
frequencies ranges co r = 0.4 to 100. Root mean square pressures are normalized on free-stream 
dynamic pressure. 

Figure 5.4 presents the root mean square pressure calculated from surface pressure 
spectra as a function of angle of attack for x/c = 0.01, 0.025, 0.04, 0.06, 0.09, and 0.14. The 
leading edge positions (1% to 6% chord) on the suction and pressure side are most affected by 
angle of attack. The pressure reaches a peak on the suction side (at the unstalled angles of 
attack) at a = 12° at 1% chord. At a = 20° the airfoil is stalled causing p' to rise substantially on 
the suction side at all locations. On the pressure side of the airfoil p' peaks at a = 18°. Again, 
the leading edge locations tend to most affected by angle of attack. 

Root mean square pressures, -^Ap' 2 calculated from pressure difference spectra are 
presented in figure 5.5. It is interesting to see that RMS Ap' 2 remain essentially flat up to a = 
14° at all chordwise locations. A slight rise occurs at x/c = 1% between 0° and 4° angle of attack. 
The ^Ap’ 2 lack of dependence on angle of attack suggests the increased eddy stretching (by the 

mean velocity field) that likely takes place with increasing angle of attack does not substantially 
increase the intensity of pressure fluctuations in the leading edge region. 

5.2 Lift Spectra and Mean Square Lift 

5.2.1 Motivation for Forming Unsteady Lift 

The effects of angle of attack on the surface pressure response can be well summarized 
by investigating the unsteady lift calculated from unsteady pressure measurements. Such a 
calculation inherently contains surface pressure and pressure difference auto- and cross-spectra 
and provides an integrated view of the airfoil response. 
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5.2.2 Unsteady Lift Formulation 


This section presents the equations for calculating lift from unsteady pressure 
measurements and addresses issues associated with such calculations. Integration of stochastic 
quantities is not a trivial matter and as such, considerable effort is spent in identifying sources of 
error and verifying results. 

Unsteady lift is calculated from pressure measurements made over the first 14% of the 
airfoil chord. This corresponds to 6 microphones per airfoil side used in the calculation. The 
microphones used in the lift calculation are denoted in figure 5.6 by red triangles and their 
chordwise and spanwise location presented in table 5.1. 


Table 5.1: Location of microphones, used in unsteady lift calculation 


y/c 

x/c 

0.33 

0.01 

0.33 

0.025 

0.33 

0.04 

0.33 

0.06 

0.33 

0.09 

0.33 

0.14 


The unsteady pressure measurements must first be resolved into normal-to-chord and tangential- 
to-chord force components as shown in figure 5.7. The relationship between these forces and lift 
is a simple coordinate transformation from the airfoil to flow aligned coordinates given as 

L' = N a cos(#)-r a sin(a) (5.2) 

where L is lift and a is angle of attack. Since the lift if a random quantity the expected value is 
taken to give the lift spectrum 

g u =e[l'-l'] 


= E[N a "N a cos 2 (a)+TjT u sin 2 (a)-T u ‘N u sin(a)cos(a)-Ai I 'r > sin(a)cos(a)] (5.3) 


= G nn cos 2 {a)+G TT s in 2 (a ) - 2 Re ( G TN ) sin (a ) co s (a ) 

Gnn and Gtt are the normal and tangential force spectrum and Gtn is the cross-spectrum between 
these forces. These spectra can be computed from unsteady pressure measurements as 


G nn (%)=£ £ £' AP(x, co r )* A p(x' , co r \lxdx' = G apap . (co r , x, xf 


lx • dx' 


(5.4) 
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G TT (a> r ) = E j o Z ”“ W P u (z(x),(O r ) + P L (z(x),co r ) dz(x)-\^ {X \p u (z(x'),co r ) + P L (z(x),(O r )^z(x y ) 


= J ( ^ * Gp [r p u '(o) r ,x,x'^+Gpp . ^(0 r ,x,x'^+2R.e(G Pu f x,x')) dz(x)-dz(x') (5-5) 


G tn M = E y * P L! (z (x) ,co r ) +P L (z (x) ,a> r ) J Jz (x) • JJ™ [P v (x \co r ) - P L (x co r )}/x ' 

= j Q Z(X) \_G Pu p v •(», ,x,x')-G PlPl '(&>. , x, x ')+G^ '(a> r ,x,x ')-G P[jPl <(a> r ,x,x ')]dz (x) • dx ' (5.6) 

where Gapap is given as 

Gapap- ( ' G>r > *> * ') = E [( P C/ (• X, 0) r ) - P L (x, 0) r )) [p v (x \<a T ) - . P L (x \<a r )) 

= E^AP(x,co r ) A P(x',<y r )J 


(5.7) 


Of course the unsteady pressure measurements are not defined spatially in a continuous manner 
and therefore, the above equations must be discretized. The final discretized form of these 
equations is presented below. 


6 6 


G,vv (p) r ) Y G\ P \p {co r ,x t ,x . )ax ; Ax . (5-8) 

i = 1 j = 1 

A yXj )+Gp L p L (oj r , X,. , x ; )+2Re(G f; ^ (© r , x, , x ; ))]Az (x f )Az(x ; ) 


6 6 


i=l 7=1 


6 6 


G Pjv(®r) = XX 

i=l 7=1 


Gp a P u ( 0J , >A yXj )-G Pl p l (0J r , X,. , Xj )+ 
Gp a P L (o\ , a , Xj )-G P u p l ( co r , X , , Xj ) 


Az(x ; .)Ax y 


(5.9) 


(5.10) 


The indices i and j run from 1 to 6 and correspond to the six microphone locations used in the 
unsteady lift calculation. 


5.2.3 Verification of Lift Formulation 

The formulation of the unsteady lift in this manner raises two important questions which 
are addressed prior to presenting the calculated lift spectra: 

1 . What is the effect of the discrete microphone spacing used to form the unsteady lift 
spectra? 
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2. What is the effect of the airfoil to flow-aligned coordinate transformation? 

Effects of Discrete Microphone Spacing on Lift: The equations used in the unsteady lift 
formulation are presented above in definite integral form. For practical implementation with a 
discrete measurement system these equations must be discretized. Consider, for example the 
discretization of the normal-to-chord force spectra Gnn 

G m («.,)=['“ P“ G„ («., , x, ■ *■ * £ £ G„„. , x, , x, K , Ax, (5.11) 

(=1 y=l 

Gapap is defined only at microphone measurement locations and therefore must be approximated 
at all other locations. Error can result from this discrete computation if Gapap is not adequately 
defined in rapidly changing regions. Consider figure 5.8 which is a contour plot of large grid 
Gapap as a function of x/c and xVc at co, = 3 for a = 0°. In the leading edge region Gapap reaches 
a peak and drops off quickly moving down the chord. This function must be defined with 
sufficient resolution for the lift computation to be accurate. 

The effects of discretization can be assessed using Amiet’s theory. The assessment can 
be done by computing the unsteady lift discretely and continuously from Amiet’s theory. That is 
Amief s theory can be used to compute the unsteady pressure at the same chordwise locations 
used in the experiment from which the unsteady lift can be computed. This lift spectrum can 
then be compared from a continuous integration of Amiet’s solution. Figure 5.9 presents a visual 
depiction of the exact and discrete techniques used with Amiet’s theory. 

The continuous unsteady lift spectrum calculation is performed in a quasi-continuous 
manner by taking a very small Ax ( = 0.14% chord). The discrete lift is calculated using the 
same chordwise microphone distribution shown in figure 5.6 (chordwise location presented in 
table 5.1). The von Karman turbulence spectrum is used in this calculation with the integral 
scale and turbulence intensity taken as the values in the large grid turbulence measurement. 

Figure 5.10 presents a comparison of Amiet’s continuous versus discrete lift calculations 
in large scale turbulence (the angle of attack is zero; i.e. Amief s solution does not allow for 
airfoil mean loading). Notice the discrete lift falls below the continuous by 0.7dB. The decrease 
of the simulated discrete lift is a consequence of cutting off the pressure peak at the leading edge. 
That is, consider figure 5.11, a plot of Amiet’s Gapap at x/c = 2.5% chord for co r = 10 as a 
function of x/c-x/c ’ (where x/c ’ is every chordwise position) for both the discrete and continuous 
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case. In the leading edge region, Gapap rises up at a high rate as a consequence of the singularity 
occurring there. This portion of the pressure distribution is not included in the discrete lift 
calculation causing the overall level to fall below the continuous calculation. Thus, it can be 
concluded that in the presence of a noise free measurement calculating lift discretely from 
pressure measurements will result in a slight reduction (< ldB) of unsteady lift. 


Effect of Airfoil to Flow-Aligned Coordinate Transformation: The next item to consider is 
the effect of the coordinate transformation implemented in the calculation of unsteady lift to 
rotate from airfoil aligned to flow aligned coordinates. When calculated unsteady lift the 
pressure difference, A P is taken across the chord (as depicted in figure 5.12a and 5.12b) such 
that, 

A/ 3 (x,)= P^x^)- P L {x^) (5.12) 

where Pu(x \ ) and Pl(x\ ) are the pressure on the upper and lower side of the airfoil, respectively 

at chordwise location x\ in airfoil aligned coordinates. However, it is desired to form the 

pressure difference in the direction of lift (also shown in figure 5.12b) such that, 

a p( x ;)= p u (x;)~ p l (x;) (5.i3> 

where Pu(x t ' \ ) and P/ix a \) are the pressure on the upper and lower side of the airfoil, respectively 

at location x a \ in flow aligned coordinates. Since the actual microphone distribution is fixed and 
thus, does not allow for pressure differencing in the direction of lift, a coordinate rotation must 
be implemented to obtain the unsteady lift. This transformation is taken as, 

G ll (co r ) = G nn cos 2 (a)+ G TT sin 2 (a)- 2 Re(G m )sin(a)cos(a) (5-14) 

and is derived above. 

The effects of this coordinate transformation can be simulated using Amiet’s thin airfoil 
theory. For some angle of attack, a the actual microphone locations (as denoted in table 5.1) on 
the top and bottom of the wing are projected into a flow aligned coordinate system as shown in 
figure 5.13. The pressure difference between these projected top and bottom points is then 
calculated. The lift spectrum, Gll is then computed by numerically integrating these values for 
each angle of attack and compared with each other. Gll is computed as follows (using the 
notation of figure 5.13), 


G„w=yy g i Ai a a** A* ■ “+vvg ( a( a 

LL\ rJ P{x nU a ]p{x mU a ) »U mU Z_i La p[x nU a )p[x mL ° ) 


^nU^mL 


+ 


£ £' g h.. i -ihw) Ax -"‘ + 2 £' g ^m-a 


(5.15) 
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Figure 5.14 presents unsteady lift calculated in this manner for each angle of attack. The 
lift spectra calculated using Amiet’s simulation show little variation (<1 dB) as a consequence of 
the coordinate transformation. Therefore, it is concluded that lift spectra are not significantly 
effected by the transformation from airfoil to flow aligned coordinates. 

Another perspective on the coordinate transformation can be taken which considers the 
phasing of the pressure response across the chord. When the airfoil is at angle of attack pressure 
sensors on the suction side of the airfoil are further away from the stagnation point than their 
companions on the pressure side. This results in a slight phase difference between the pressures 
measured on each side of the airfoil as shown in figure 5.15. The phase difference can be 
modeled by considering the pressure fluctuations to be harmonic. Referring to the notation of 
figure 5.16 the pressure at a = 0° is given as 

J'M 

AP(a = 0) = 2Pe 1 u "> (5.16) 

and the pressure at non-zero angle of attack as 

f d„ +x n \ ( d„ —x n \ 

A P(a) = Pe 1 * >+Pe ' u ° K (5.17) 

Now, by taking the ratio of eq. 5.16 and eq. 5.17 the effect of the phasing difference is revealed 
as 


AP(a) 
AP(« = 0) 


Pe 


• i d.+x n 
l(0\ Z-fA s " 

U rri 


+ Pe 


■ i d. -x n 
ico\ t + - 1 — ^ 
U m 


2Pe 


= cos 


f (OX, ^ 




(5.18) 


In forming the power spectral density of pressure fluctuations the pressure difference is 
multiplied by its conjugate and, as such, the cosine term of eq. 5.18 is squared. If x 0 is assumed 
constant then integrating along the chord (as performed to obtain unsteady lift) will not affect the 
result of equation 5.18. A plot of the phase shift model (equation 5.18 squared) as a function of 
frequency for each angle of attack is presented in figure 5.17. The phase shift model shows no 
reduction in unsteady lift for co r < 10; however, reduction in unsteady lift above co, = 10 is 
possible based on this model. 


Unsteady Lift Spectra: Lift spectra, Gll calculated from unsteady pressure measurements are 
presented in figure 5.18 with spectral level plotted against ox. The results are very surprising in 
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light of previous mathematical and experimental work on this problem. Two frequency 
dependent interaction regions are shown in this figure. At low frequencies (ny < 10) a 
suppression of unsteady lift with increasing angle of attack occurs with a maximum lift reduction 
of 5dB. The opposite effect occurs at high frequencies ( co, >10) with a rise in unsteady lift of up 
to lOdB occurring with increasing angle of attack. The reduction in unsteady lift at low co, is a 
never before shown result and unexpected based on previous experimental studies, in addition to 
not being predicted by theoretical formulations. McKeough (1976) showed, through unsteady 
lift measurements on a NACA 0015 in turbulence an increase in unsteady lift with an increase 
from a = 0° to 10°. Patterson and Amiet (1976a) discussed, although did not present supporting 
data, an increase in surface pressure spectral magnitude with increasing angle of attack. This 
result is also not supported by theoretical formulations which account for airfoil mean loading, 
such as the theories developed by Reba and Kerschen (1997) and Atassi (1984). It should be 
noted, the lift spectrum at a = 20° shows significant unsteadiness at low co r which results from 
large scale pressure fluctuations associated with stall. Airfoil stall at this angle of attack is 
corroborated by mean pressure measurements and flow visualization. 

Effects of Forming the Spectra over the First Fourteen Percent Chord: With the unsteady 
lift presented one last issue remains to be examined. The effects of limiting the integration of 
unsteady pressure to the first 6 microphones (i.e. up to 14% chord) must be examined to ensure 
the trends with angle of attack are not influences by this bounding. Effects which may result 
from this bounding are examined through comparison of unsteady lift spectra calculated using all 
thirteen chordwise microphones with lift calculated from the first 6 chordwise microphones. 
Using all thirteen microphones will result in integration along 85% of the airfoil chord. 

Figure 5.19 shows lift spectra calculated from partial and full chord integration. Lift 
spectra calculated over 85% chord show the same angle of attack observations revealed in 
spectra formed with the first 14% chord. Lift spectra calculated over the first 85% chord are not 
as smooth as the 14% chord lift spectra due to turbulent boundary layer pressure fluctuations and 
trailing edge separation. Furthermore, the inclusion of more measurement locations into the 85% 
chord lift calculation causes the overall lift spectra level to rise. Large scale, low frequency 
pressure fluctuations associated with stall are clearly revealed by the increase in unsteady lift 
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level at low frequency at a = 20° when the unsteady lift is calculated from all microphones (85% 
chord). 

5.2.4 Lift Spectra at Zero Angle of Attack 

The a = 0° lift spectrum, Gll is presented in figure 5.20 along with that calculated from 
Amiet’s (1976a, b) theory. Over low co r (< 5) the measurement and prediction are within 3dB 
and the agreement in shape is very good. The measured Gll begins to fall off with a greater 
slope than the predicted for 5 < co r < 20. The spectral levels are within 13dB in this region. The 
difference in slope of the predicted Gll is likely related to the zero thickness flat plate model 
used in this theory. The predicted and measured Gll crossover at co r = 120 as a consequence of 
the change in slope of the measured Gll ■ The change in slope which occurs at higher co r is 
related to turbulent boundary layer pressure fluctuations which become significant at higher 
frequencies and tend to dominant the unsteady lift over this frequency range. The peak in Gll 
which occurs at = 600 is related to resonance of the microphone-mounting hole system. 

5.2.5 Comparison of Predicted and Measured Unsteady Lift at Angle of Attack 

The lift spectra calculated from pressure measurements are compared with lift spectra 
computed from Reba and Kerschen’s (1996) Rapid Distortion Theory in figure 5.21. Lift is 
calculated for a = 0°, 2°, and 6° using Reba and Kerschen’s theory. A direct comparison of 
measured and theoretical unsteady lift levels is not appropriate given the spanwise wavenumber 
integration restriction (discussed in section above); however, a comparison of spectral shape and 
trends with angle of attack are still of substantial value. For small angles of attack (a < 6°), the 
measured Gll falls by less than 1.5dB. Over the remaining angles of attack the lift falls up to 
5dB (note, at a = 20° the airfoil is stalled and, as such, a significant rise in unsteady lift occurs at 
low o\). The RDT prediction of Reba and Kerschen shows a substantial increase in unsteady lift 
over a similar range of co r . The large increase (as much as 14dB) appears to occur as a 
consequence of the flat plate model assumed in the theory. The velocity gradients associated 
with the leading edge singularity of a flat plate at some a in inviscid flow are much larger than 
those of an airfoil with thickness. As such, the inflow is not distorted as it would be when 
encountering a real airfoil with thickness. 
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The measured lift spectra begin to converge and show no effect of angle of attack over 5 
< co r <10. Over the remaining oj, range Gll remains unchanged for small a’s (a < 6°). The lift 
spectra’s lack of dependence on angle of attack over this a> r range agrees well with the 
measurements made by McKeough (1976). The Gll predicted by Reba and Kerschen over this 
co r range increase with increasing angle of attack by as much as 6dB; however, the spectra appear 
to be converging at higher ny. Again, in the presence of an airfoil with thickness it is likely that 
the theory of Reba and Kerschen would not predict as large an increase in Gll with increasing 
angle of attack. In fact, the convergence of the lift spectra over this co r range may be associated 
with the airfoils finite thickness. 

For co r > 10 the slope of the measured Gll begins to decrease as pressure fluctuations 
associated with the turbulent boundary layer and hydrodynamic pressure fluctuations begin to 
affect the unsteady lift. Comparisons with Reba and Kerschen’ s theory were not possible over 
this frequency range due to an unreliable complex error routine used in their FOTRAN code. 

5.2.6 Mean Square Lift at non-Zero Angle of Attack 

The lift spectrum calculated from the first 14% chord is integrated at each angle of attack 
to obtain z 2 values as 

L 2 (a) = \G Ll (a,f)lf (5.19) 

The mean squared lift if obtained by integrating over = 0.35 - 10. Mean square L 
values normalized on ( qc ) , where q is the dynamic pressure and c the chord are presented in 
figure 5.22 as a function of angle of attack. As pointed out by McKeough and Graham (1980) 
angle of attack effects should appear first in terms of order a . Therefore, figure 5.22 also shows 
a curve representing an a variation for comparison with the data. As the angle of attack is 
increased the mean square L decreases nearly as a up to a = 14°. Agreement with the a curve 
is within 6% up to a = 14°. The coefficients, A and B are calculated to be -0.00137 and 0.5045 
respectively. Unsteady lift associated with separation becomes significant at a = 16° causing the 
rate of decrease in z 2 to fall. At a = 20° the unsteadiness in the lift is due primarily to stall 
causing an apparent increase ini 2 . 
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5.3 Admittance Function 

An airfoil interacting with turbulence can be treated as a linear system (as is done with 
most theoretical representations) and as such the response (unsteady lift) can be written in terms 
of the excitation (inflow turbulence) multiplying a transfer function (airfoil response). If the 
output is taken as the unsteady lift, the transfer function is often referred to as an admittance 
function and can be written as 


®c,M = 








(5.20) 


where O ww (<u) is the upwash velocity spectrum, |d(&»)| 2 is the admittance function and <J> Ci (co) 
is the normalized lift spectrum defined as 


® c, (®) ('it (®) 


pUlb 


(5.21) 


Based on velocity measurements presented in Bereketab et al. (2000) the energy spectrum of the 
inflow turbulence can be represented over the entire wavenumber range (including the inertial 
subrange) using the interpolation of Panchev and Kesich (Hinze, 1975, pp. 241) giving 
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where f a is given as 
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The upwash is then determined as (Hinze, 1975) 
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(5.24) 
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Integrating over k z and k y gives the one-dimensional velocity spectrum used in equation 5.20. 
The wavenumber, k x can be related to angular frequency, co by Taylors hypothesis, k x = 
UJco. Figure 5.23 presents the von Karman and modified von Karman spectrum. Note the von 
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Karman spectrum is referred to as Amiets von Karman in figure 5.23 and is calculated from 
Patterson and Amiet (1976a, pp 57). The modified von Karman is calculated from equation 
5.22. 

The admittance is extracted using equation 5.20 and compared with the experimental 
results of McKeough (1976) for large scale turbulence interaction in figure 5.24. Recall, 
McKeough made unsteady lift measurements on a NACA 0015 in turbulence at a = 0° and 10°. 
The scale of the turbulence was 40% of the chord with intensity of 5.8%. For all angles of attack 
the admittance forms a ‘V’ shape, dropping from oj, = 1 0 to oj, =25 and then rising up. The rise 
in admittance at high frequency may be related to eddy convection or boundary layer pressure 
fluctuations. That is, the unsteady lift spectra slope decreases at the higher frequencies (as a 
result of eddy convection and boundary layer pressure fluctuations) while the turbulence 
spectrum continues to drop off at a constant slope as shown in figure 5.23. 

The angle of attack variation is preserved in the admittance function with as much as 5dB 
of decrease observed at low reduced frequencies with increasing angle of attack. At a = 0° the 
slope of the flat region (co r = 2 to 20) is in agreement with McKeough’s data; however, the data 
of the present measurement is shifted to the right as a consequence of the smaller scale 
turbulence (A/c = 13%). Surprisingly, the magnitude of the admittance of the present 
measurement is above McKeoughs. McKeoughs data also shows an increase in admittance 
magnitude with increasing angle of attack. At a = 10° the slope of McKeoughs data agrees well 
with the present data at a = 8°. 

5.4 Surface Pressure Correlation Length Scales 

The spanwise and chordwise correlation length scales are presented to offer additional 
insight into the source of the angle of attack effect observed in unsteady lift spectra and 
admittance function. Furthermore, comparisons are made with Amiet’ s (1976a, b) theory to 
determine the ability of his solution to accurately predict correlation length scales at zero angle 
of attack. 

5.4.1 Spanwise Correlation Length Scale 

The spanwise correlation length scale is computed at the five most leading edge 
measurement chordwise locations (x/c = 1, 2.5, 4, 6, 9%). The spanwise correlation length scales 
are computed from the pressure cross-spectrum as 
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(5.25) 


J 0 S 9g (x,x',Tj = 0,f) 

where L q is normalized on the airfoil chord, Jj max = 96% chord, and is the pressure difference 
cross-spectral density. This calculation is carried out at each chordwise location from 1% to 9% 
chord with x = x' for all frequency. 


Zero Angle of Attack: Figure 5.25 presents L n calculated from measured unsteady pressure and 
as predicted by Amiet’s (1976a, b) theory as a function of frequency at x/c = 1, 2.5, 4, 6, and 9% 
for a = 0° (note: high frequency data was only taken in spanwise rows at 1 and 2.5% chord). At 
1% chord the spanwise correlation length of the inflow turbulence is also presented and 
calculated from the von Karman spectrum (Amiet, 1975) as 


(5.26) 
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Reasonable agreement between L t and the measured 


L n at 1% chord occurs for 1 0 < co r < 30. This result is seems logical since the surface pressure 
length scale in the leading edge arises from the turbulence-airfoil leading edge interaction. 
Moving away from the leading edge, the measured is shown to increase over the entire 
frequency range. This is an expected result since the change in airfoil lift occurs by propagation 
of acoustic waves. That is, the acoustic wave that forms at the leading edge propagates in all 
directions with a spherical wave front. Amiet’s theory under predicts the spanwise length scale 
below co r = 10 for all chordwise positions which is likely related to the flat plate model used in 
this theory. At 1 and 2.5% chord agreement within 2% chord is shown to occur for co r > 20. The 
data begins to fall below Amiet’s prediction for co r > 20 beyond 2.5% chord. At these higher 
frequencies boundary layer pressure fluctuations (high frequency), which correlate over very 
small spanwise separations (< 1% chord) begin to contaminate L q by driving up S qq at rj = 0 and 
thus, causing L n to fall. In summary, Amiet’s (1976a, b) theory is limited in its ability to 
accurately predict the spanwise correlation when the airfoil is not well modeled by a flat plate. 
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Effects of Angle of Attack: Investigating the variation of the spanwise correlation with angle of 
attack offers additional hints on the source of the angle of attack effect. Figure 5.26 presents the 
spanwise correlation length scale as a function of reduced frequency for each angle of attack (a 
= 0°, 4°, 8 °, 12°, 16°, 20°) at x/c = 1 and 2.5% chord. At both chordwise locations a strong 
dependence on angle of attack is shown in L, r As the angle of attack is increased L n falls over 
much of the frequency range; however, most dramatically over < 10. The reduction in the L n 
in the leading edge region with increasing angle of attack suggests that eddies are being 
increasingly stretched by the mean field producing turbulence with a smaller spanwise length 
scale. This type of phenomenon may be the mechanism producing the reduction in unsteady 
loading at low reduced frequencies. Figure 5.27a, b depicts how such a process could reduce the 
unsteady loading. At small a’s the inflow eddy (labeled ‘A’) is convected towards the leading 
edge of the airfoil and moderately stretched by the mean field velocity gradients (eddy labeled 
‘B’) as shown if figure 5.27a. Now, at a higher angle of attack, the same eddy ‘A’ is stretched 
considerably more as a consequence of the increased rate of strain in the mean field which 
reduces the spanwise length scale to that shown by eddy ‘C’ (figure 5.27b). Since the spanwise 
length scale of the inflow is related to the spanwise correlation length scale of the surface 
pressure, as suggested in figure 5.25, a reduction in the surface pressure L v occurs with 
increasing angle of attack. The reduction in L n means that the pressure has less of a correlated 
area to act over. Since the intensity of pressure fluctuations appears to remain largely unaffected 
by angle of attack, the overall unsteady loading decreases with increasing angle of attack as a 
consequence of the reduction in correlated area. 

5.4.2 Chordwise Correlation Length Scale 

Correlation length scales are typically defined for quantities which are homogeneous (i.e. 
ergodic) in space and time. The cross-spectrum of surface pressure fluctuations is indeed 
homogeneous in the spanwise direction, and, as such the distance between two spanwise points 
( 77 ) becomes the defining parameter. However, in the chordwise direction pressure fluctuations 
are inhomogeneous (for example, this is apparent through comparison of chordwise auto-spectra) 
and thus the absolute chordwise position relative to any other chordwise position must be 
considered. Therefore, it is not necessarily appropriate to compute the chordwise correlation 
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length scale. Nevertheless, this quantity is essential for computing far-field noise from surface 
pressure fluctuations and is thus presented here. 

The chordwise correlation length scale is calculated much like the spanwise correlation as 


y s„(x,x\'i=o,f) 
c ~ l s, s (w=o,/) 

with L c normalized on the airfoil chord, x = x', x max = 85% chord, and S qq 
difference cross-spectral density. 


(5.27) 


is the pressure 


Zero Angle of Attack: Figure 5.28 presents L c as function of frequency calculated from 

unsteady pressure measurements and that predicted by Amiet’s theory at x/c = 1, 2.5, 4, 6, and 
9% for a = 0°. With increasing chordwise distance from the leading edge L c increases nearly the 
same amount at all frequencies up to 6% chord. At 9 and 14% chord L c continues to increase for 
co r < 10 while beginning to decrease at higher frequencies. The increase in L c is, as with the 
spanwise correlation length scales an expected phenomenon due to the propagation of acoustics 
waves as a consequence of lift adjustment. The reduction shown to occur over co r > 10 at 9 and 
14% chord is likely due to contamination from boundary layer pressure fluctuations much like 
that which affects the spanwise correlations. Recalling the length scale of the inflow turbulence, 
X/c = 0. 13, it is interesting to note that L c is in close agreement with the inflow length scale at 1% 
chord for low reduced frequencies. 

Amiet’s theory predicts L c to within 5% c at all frequencies up to 6% chord. The theory 
begins to over predict L c for frequencies above 10 at 9 and 14% chord. This may be a combined 
result of the flat plate model assumed in Amiet’s theory and contamination of the data through 
boundary layer pressure fluctuations. 


Effects of Angle of Attack: Angle of attack effects on the chordwise correlation length scale 
are presented in figure 5.29. It is interesting to see that variation in the airfoils mean loading 
does not affect L c nearly as significantly as shown in the spanwise correlation. Furthermore, at 
1% chord L c increases slightly for co r < 10. At 2.5% chord slight reduction in L c does occur as 
the angle of attack is increased to a = 16°. When a = 20° L c rise up significantly since the airfoil 
is stalled and large scale pressure fluctuations develop at low frequencies. 
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5.5 Correlation of Velocity and Pressure 

Examining the correlation of pressure and velocity will aid in identifying both 
hydrodynamic and acoustic pressure fluctuations and may offer additional insight into the 
laminar separation bubble theory presented in chapter 4. Recalling that for some angle of attack 
simultaneous measurements of velocity and pressure were performed, correlation of these two 
quantities is presented here. Details on the instrumentation associated with the velocity-pressure 
measurement can be found in chapter 2; although briefly, hotwire measurements were made 
(simultaneously with pressure measurements) at 14% chord on the suction side (suction side 
defined for positive a’s) of the airfoil. The hotwire was located Y" off the airfoil surface. The 
velocity-pressure space-time correlation is formed as 


, s p(x,t)u(t + T) 

R x,r (5.28) 

pu v ’ nTTi 
V P \ x ) u 

where u is the fluctuating velocity measured by a single hotwire located at 14% chord. 
Investigating the correlation of velocity and pressure assists in determining whether pressure 
fluctuations arise from convected or propagated disturbances. In this way, the correlation 
function will help determine if vortex shedding is occurring as a consequence of an inflow 
coupled laminar separation bubble and thereby altering the airfoil loading. 

The correlation function is presented in figure 5.18 for a = 0°, ±8°, ±16°. Superimposed 
on the pressure-velocity correlation contours are gradient vectors computed by taking the 
derivative of the correlation with respect to time and space. These vectors assist in identifying 
ridges in correlation. At a = 0° on the pressure side of the airfoil a correlation ridge forms 
(indicated by the red dashed line) at a 12:1 slope with significant pressure-velocity correlation up 
to 85% chord indicating propagation of sound speed pressure fluctuations. Bearing in mind the 
location of the hotwire on the airfoil suction side, the correlation will not show convection on the 
pressure side. However, on the suction side of the airfoil both convected and propagated pressure 
fluctuations are shown in the correlation function. The convect disturbances appear to be 
moving at a velocity slightly less than the mean ffee-stream. At a = 8° the correlation looks 
much like that at a = 0°. As the angle of attack is increased to 16° a significant inviscid response 
is shown on the pressure side of the airfoil with convection occurring near the mean free-stream 
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velocity on the suction side. The convection may be related to vorticity shedding from a leading 
edge transient laminar separation bubble. 

The correlation is changed dramatically at a = -8° and -16° with the hotwire now located 
on the pressure side of the airfoil. At these angles of attack the correlation shows little evidence 
of a strong inviscid response. Some possible convection is evident (indicated by the solid red 
line) and appears to be slightly faster than the mean free-stream velocity. Eddies impacting the 
leading edge and then convecting along the pressure side of the airfoil appear to either be 1) 
substantially distorted such that they pass below the V" gap between the airfoil surface and 
hotwire or 2) have no significant component which survives after impacting the leading edge. If 
in fact there is no significant component of the eddy surviving, then a Rapid Distortion Theory 
(RDT) model could be quite simplified such that the distortion would only need to be considered 
in neighborhood of the leading edge. 

5.6 Effects of a Modified Trailing Edge 

To examine the effects of trailing edge geometry, measurements were performed with a 
modified trailing edge. The modification transforms the NACA 0015’s rounded trailing edge to 
one that is sharp (see chapter 2 for geometry details) and extends the airfoil chord by 3%. The 
rounded trailing edge allows for infinite number of velocity solutions to the inviscid flow around 
the airfoil. Of course in steady, inviscid flow only one solution occurs in reality; however, in 
unsteady flow little evidence of the applicability of the classical Kutta condition is available. 
Furthermore, the experimental work of Poling and Telionis (1986) suggests that the “classical 
Kutta condition” is never satisfied in unsteady flow. Indeed, they show that in unsteady flow at 
low a>r the shedding of vorticity results in a finite velocity at the trailing edge. As such, the 
unsteady circulation may be slightly altered relative to that which would occur if the “classical 
Kutta condition” were met and would clearly impact the unsteady loading of an airfoil 
encountering turbulence. 

5.6.1 Correlation of Pressure in Space and Time 

The space-time correlation is investigated with a modified trailing edge in Figure 5.31 for 
each angle of attack (a = 0°, 4°, 8°, 12°, 16°, 20°). The correlation is formed by taking x = 1% 
chord on the suction side (negative x/c values) and varying x' with 7=0. The data is high pass 
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filtered at 100Hz to eliminate uncertainties associated with microphone calibrations below this 
range (i.e. this data is processed the same way as the data presented without modified trailing 
edge in chapter 4). The inviscid response dominates on the pressure side of the airfoil at all 
angles of attack with the ridge of the correlation nearing a 12:1 slope. On the suction side the 
inviscid response is dominant up to a = 8°; however, convection (1:1 slope) of eddies begins to 
become significant at higher angles of attack. Significant separation occurs at a = 16° causing 
substantial convection on the suction side of the airfoil. The airfoil is stalled at a = 20° and 
shows little correlation on the suction side. 

The space-time correlation measured with the modified trailing edge compares very 
closely with that measured with an unmodified trailing edge (figure 4.34) at all angles of attack. 
This finding indicates that small changes in the chord length impact the inviscid response very 
little. This occurrence makes sense given that the dominant energy in the inflow turbulence is in 
wavenumbers much smaller than the chord length. Additionally, the sharpness of the trailing 
edge appears to impact the inviscid response minimally suggesting that the “classical Kutta 
condition” can be reliably applied in unsteady flow (as assumed in unsteady thin airfoil theory). 

5.6.2 Unsteady Lift 

A comparison of unsteady lift calculated from measurements with both types of trailing 
edges will assist in understanding the effects of trailing edge geometry. The unsteady lift spectra 
calculated from the modified trailing edge pressure measurements are compared with those 
calculated from unmodified trailing edge measurements in figure 5.32. Little variation (within 
the uncertainty of the data) between the spectra calculated with and without the modified trailing 
edge occurs across the frequency range. This result confirms that indeed, a rounded trailing edge 
(within the parameters of this experiment) behaves much like a sharp trailing edge in unsteady 
flow. Additionally, the slight increase in chord length associated with the sharp trailing edge 
seems to affect the airfoils response negligibly. 

5.7 Sources of Low Frequency Angle of Attack Effect 

The analysis of unsteady pressure measurements presented here and in chapter 4 seems to 
suggest two reasonable sources for the observed reduction in unsteady lift at low co,-. The 
transient laminar separation bubble suggested by oil flow visualization and leading edge pressure 
asymmetries is one means for the relaxation in unsteady loading. Another possibility is the 
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reduction in spanwise correlation length scale as a result of inflow distortion. The likelihood of 
each is examined in this section. 

5.7.1 Laminar Separation Bubble 

A possible source of the reduction in unsteady lift at low reduced frequency may the 
presence of a laminar separation bubble in the leading edge region of the airfoil. The flow 
visualization presented revealed the possibility of a laminar separation bubble occurring in the 
leading edge region. Some additional support of such a phenomenon is offered by the 
asymmetries shown in pressure spectra. Also, the reduction in the pressure-pressure correlation 
shown to occur on the suction side in the leading edge region may be related to the formation of 
such a bubble. This bubble, if indeed present could serve as an explanation for the reduction of 
unsteady lift with increasing angle of attack. At low reduced frequencies the bubble could be 
synchronizing with the unsteady lift such that the bubble sheds vorticity and thereby relaxes the 
unsteadiness in the lift. There is little evidence in the RMS pressure to indicate the presence of a 
laminar separation bubble; however, this alone is not evidence enough to eliminate such a 
possibility but, also does not aid in the support of this theory. Careful measurement of the 
velocity field surrounding the airfoil coupled with improved flow visualization (such as helium 
bubble flow visualization) would allow for unambiguous proof (or disproof) of such a scenario. 

5.7.2 Distortion of Inflow 

Analysis of the spanwise correlation length scale suggests distortion of turbulence as a 
consequence of the airfoils presence in the flow is a significant factor in reducing unsteady lift 
with increasing angle of attack. The spanwise extent of correlated pressure fluctuations in the 
leading edge region appears to be produced by the inflow turbulence as suggested through 
comparison of the undistorted inflow spanwise turbulence length scale ( L t ) with the spanwise 
pressure correlation length scale (L n ). As such, with the airfoil at angle of attack, stretched 
eddies (consequently, with a reduced L t ) impact the leading edge thereby reducing L n . The 
reduction in L n reduces the correlated area over which the pressure acts. Since the intensity of 
pressure fluctuations appears to remain largely unaffected by angle of attack (as shown in figure 
5.5), the overall unsteady loading decreases with increasing angle of attack as a consequence of 
the reduction in correlated area. This theory may be the most conceivable of the two presented 
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based on the observed variation of L v and \]Ap ' 2 with angle of attack which offer particularly 
strong support. 

Furthermore, while the distortion of the inflow appears to be critical in this experiment, 
wherein the integral scale (A/c) is 13% of the chord, this effect does not seem to be significant in 
large integral scale (relative to the chord) flow. Recall McKeough’s measurement of unsteady 
lift on a NACA 0015 which showed a 3dB (for co r < 1 ) increase in lift with an increase from a = 
0° to 10° and had an inflow integral scale of 40% of the chord. As such, it appears the relative 
scale of the inflow to chord length is important in determining the airfoils response when subject 
to mean loading. Unsteady loading models like those developed by Graham (1970a), Atassi 
(1984), and Reba and Kerschen (1996) seem to encapsulate mean loading effects well for this 
large scale turbulence case. In fact, Graham’s (1970a) theory agrees well with McKeough’s 
unsteady lift measurements based on comparisons he presented. However, for the case of small 
scale turbulence (relative to the airfoil chord), these models do not depict the dominant physics 
of turbulence-airfoil interaction well. The reason seems to be related to the flat plate with 
leading edge stagnation point assumed in these models that does not account for the inflow 
distortion correctly. 

5.8 A Model for Predicting Angle of Attack Effects 

A new model is suggested here that will correctly account for the inflow distortion and 
therefore, the dominant factor impacting the airfoil response for the parameters of this 
experiment ( A/c = 13%). This new model combines accurate modeling of the turbulence 
distortion using Rapid Distortion Theory (RDT) with Amiet’s (1976a, b) flat plate unsteady 
loading theory. This type of modeling scheme will require as input the distorted turbulence 
energy spectrum to account for the effects of angle of attack. 

Rapid Distortion Theory Background: Rapid Distortion Theory (RDT) is a scheme for 

predicting the evolution of turbulence as it is stretched and compressed by the mean velocity 
field. The foundation for this theory was laid by Batchelor and Proudman (1954) and uses the 
linearized form of the vorticity transport equation to predict distortion effects. This theory works 
under the assumption that 

1 . Distorting effects of turbulence on itself are negligible 
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2. The scale of the turbulence is small compared to the scale of the distortion 

3. The duration of the distortion is small, such that the influence of viscous dissipation can be 
neglected. 

In Lagrangian form the vorticity transport equation, or what is often referred to as 
Cauchy’s equation is given as (in tensor notation) 

dx, 


Q, = 


q; 

da, ] 


(5.29) 


where Q ' is the vorticity at the beginning of the distortion (time f) with position determined by 

Ox- 

a.j and Q ; is the vorticity at the end of the distortion (time t) with position xj. The term — L is 

dcij 

referred to as the distortion tensor and will depend on the mean velocity field. Batchelor and 
Proudman (1954) define the distortion tensor as 


_ dx 
ij da 


L = s y + 


rdU, 


l da , 


L dt 


(5.30) 


where aj is the initial position of the particle at time t\ x, is the position at time t, U ) is the mean 
velocity component and Sy is the Kronecker delta. Using Goldstein’s (1979) model the distortion 
tensor can be computed from a two dimensional mean velocity field as 


S 9 = 


dX x 
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dX x 

ds 

dn 
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0 

dX \ 
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dX 3 

ds 

dn 


-i-i 


(5.31) 


Note that in eq. 5.31, the mean field does not vary in the y direction (i.e. along the span of the 
airfoil) and hence the zero and unity terms. In this equation 5 and n form a curvilinear coordinate 
system with 5 along a streamline and n perpendicular to the streamline. The quantity X3 is 
proportional to the stream function with X\ defined as the drift function and given as 


u„ 


- = — + p 

u J -°° 


1 


1 


U\x v X,) U_ r _ 


\dx. 


(5.32) 


where here U ’ is the mean velocity component parallel to the ffee-stream velocity (14, ) and xi is 
the component of the fluid particle position vector parallel to U«>. Through manipulation of the 
drift function, the distortion simplifies to 
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s „ = 


U, 
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— 0 
U x dn 

0 1 0 
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0 


u. 


(5.33) 


where 


U, 


dX l 


is the magnitude of the mean velocity vector. The term — - is difficult to compute 

dn 


because it implies integration along a streamline (eq. 5.32) followed by differentiation 
perpendicular to the streamline; however, this calculation is feasible with careful application of 
numerical procedures. 

Batchelor and Proudman (1954) show that the distorted turbulence energy spectrum can 
be calculated using the distortion tensor and the undistorted energy spectrum as 


®’ii-disXX)dX = £ikl S rnn P S iab £ u W ,XlXbX^KK S kJau® P S k ) dk 


(5.34) 


with 27 = tiS.y' is the non-dimensional distorted wave number vector, 5, y is the distortion 
tensor, Syk is the unit alternating tensor and 0 ;; (£) is the non-dimensional undistorted energy 
spectrum tensor. 


New Airfoil Response Model: A new unsteady response model is developed which uses a 
circular cylinder to predict the inflow distortion. Ideally, the distortion tensor is calculated based 
on the actual mean field produced by the airfoil; however, using a cylinder has the distinct 
advantage of significantly simplifying the distortion tensor calculation. The circular cylinder 
radius is determined based on a fit to the airfoil curvature at the location of the mean stagnation 
point for any given angle of attack. Figure 5.33 shows a circular cylinder fit to the NACA 0015 
at a location that corresponds to the a = 4° stagnation point. The location of the mean stagnation 
point is easily calculated using a vortex panel method. 

The stagnation streamline associated with this circular cylinder is then used to predict the 
distortion of the inflow via RDT (eq. 5.34) as depicted in figure 5.34. The first step in this 
calculation is determining the stagnation streamline velocity which is obtained for a circular 
cylinder from potential flow theory as 
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U(x) 


U„ 


= 1 - 


r r \ 2 


yxj 


(5.35) 


where r is the cylinder radius. The distortion tensor is then readily computed from the potential 
flow solution using eq. 5.33 




U(x) 
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0 
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0 
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t/ 


U(x) 


-1 


(5.36) 


dX, 


Notice the most difficult term to compute in the distortion tensor, — L is zero when calculated 

dn 

along the stagnation streamline of a circular cylinder and hence, the real convenience of using 
the cylinder approximation. 

This calculation is carried out for a = 0°, 4°, 8°, and 12°. The distortion is compute at a 
location, x (referring to the notation of figure 5.34) that is some percentage of the integral scale 
(A/c = 13%) upstream of the cylinder such that 

(5.37) 

The constant, A 2 is taken as 1 for the present calculation resulting in calculation of the distorted 
turbulence spectrum one integral scale upstream from the circular cylinder. The ratio of 
distortion location, x to cylinder radius, r for each angle of attack is given in table 5.2. 

Table 5.2: Ratio of distortion location to cylinder radius 


x = r + A 2 / l . 


Angle of 
Attack 

x/r 

0° 

-9.81 

4° 

-3.14 

8° 

-2.11 

12° 

-1.71 


With the distorted turbulence spectrum computed, Amiet’s (1976a, b) theory is used to predict 
the airfoils response to the distorted inflow at each simulated angle of attack. 

The results of this calculation are presented in figure 5.35 with lift spectra, Gll plotted 
against co r for each angle of attack. The lift spectra are calculated discretely using the same six 
chordwise locations used to calculate lift from measured unsteady pressure. This new unsteady 
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loading prediction scheme shows a 6dB reduction is lift for co r < 10 as a is increased from 0° to 
12°. This is excellent prediction of the angle of attack trend and surprisingly good prediction of 
amount of spectral level reduction. In light of the approximations used in the formulation (i.e. 
calculation of distortion one integral scale upstream of a circular cylinder), prediction of the 
correct change in spectral level with angle of attack was not expected. This model also predicts 
some reduction in unsteady lift for uy > 10 which is not shown in the measurement. Clearly 
though, this model provides solid evidence that, in fact, the distortion of the inflow is crucial in 
producing the reduction of unsteady lift with increasing angle of attack. 

To shed additional light on the relative importance of the inflow distortion for large scale 
turbulence (large relative to chord length), this model is used to predict the unsteady lift based on 
McKeough’s flow conditions (£/«, = 15m/s, X/c = 0.4, 77= 6.5%). This calculation is carried out 
as outlined above with x/r = 3.1 and 27.2 for a = 0° and 12° respectively. Figure 5.36 presents 
the predicted unsteady lift spectra as a function of co r . Less than ldB of reduction in lift spectral 
level is predicted for co r < 10 with the spectra falling nearly on top of each other for higher 
frequencies. As expected, this theory does not predict the increase in unsteady lift shown to 
occur in McKeough’s measurements. However, this prediction does prove the point that the 
inflow distortion is not an important factor in determining the unsteady loading (for X/c = 0.4) 
and, in fact, those factors which are correctly modeled in Graham’s (1970a) theory are likely 
dominant. 

5.9 High Frequency Angle of Attack Effects 

The unsteady lift calculated from pressure measurements (Figure 5.18) shows an increase 
in spectral with increasing angle of attack for co r > 10. This increase in unsteady loading can be 
explained by scaling lift spectra. Lift spectra are scaled on the square of the mean lift as 

g u(«./)= G -frf ) (5.38) 

L (<z) 

—2 

where G LL is the unsteady lift as calculated from eq. (4.3) and L is the mean lift squared. 

Figure 5.37 presents the scaled lift spectra (note, data for a = 0° is not presented due to 
singularity which results from zero mean lift) as a function of reduced frequency. As expected, 
this scaling amplifies the angle of attack effect at low frequencies (co r < 10) causing a greater 
reduction in unsteady lift with increasing angle of attack. However, above uy = 10 the lift 
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spectra are collapsed to within 5dB by the mean lift scaling for a > 4°. This is a considerable 
reduction in the lOdB of variation with angle of attack shown over this frequency range to occur 
in the unsealed lift spectra. The collapsing of data over this frequency range suggests the 
variation in unsteady lift here is closely related to the mean loading on the airfoil. 

Furthermore, this scaling suggests that effects of angle of attack are frequency dependent 
with the distortion of the inflow defining the low frequency effects and the airfoils mean loading 
producing the high frequency effect. In the present experiment the crossover appears to occur 
near co, = 10 . 


5.10 Computing Far-field Noise from Surface Pressure 


Lastly, because it is the ultimate goal of any surface response theory to generate 
predictions of the far-field noise produced during turbulence-airfoil interaction the usefulness of 
this data for such schemes is presented. Casper and Farassat (2002) present a technique for 
computing far-field loading noise (i.e. inviscid response noise) from the unsteady surface 
pressure distribution. The solution depends on the time derivative and surface gradient of the 
unsteady pressure. The required input for their formulation is obtained by expressing the surface 
pressure as the cross-correlation of the substantial derivative of p, such that 


p(x,y,t) 

and the cross-correlation is given as, 


dp(x,y,t) u dp(x,y,t) 


dt 


dx 


(5.39) 


R i , p (x,x\jj,r) = p(x,y,t)p(x',y',t + T) . 


(5.40) 


The pressure difference inverse Fourier transform of eq. (5.38) — the convective derivative 
spectrum — is considered here for zero spanwise separation (i.e. rj= 0 ) and is obtain as 


UU 

Got (*> x 7, /) = | R Am (x,x\ r, 7 , r )e" ! ' 2 * /r dr 


(5.41) 


where the left hand side of eq. (5.39) is computed from time series data as 
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G Dl (x,x’,f) = 


n-s d 
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(5.42) 


with 
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(5.43) 


The spatial and temporal derivatives of eq. (5.41) are evaluated numerically such that 

dAp(x,t) ^ Ap(x i} t)-Ap(x i+2 ,t) 
dx x i+2 -Xi 

dAp(x,t) ^ Ap(x, t i ) - Ap(x,t i+2 ) 


(5.44) 


The time derivative, Gdt and space derivative, G)/ v , both components of the convective derivative 
spectrum, are also considered and can be directly calculated from eq. (5.41) as 


G dl (x,x\f) 
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(5.45) 


and 
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(5.46) 


The above quantities are calculated for a = 0° at x/c = 1, 2.5, 4, 6, 9, 14% chord and 
presented in figure 5.38 as a function of reduced frequency along with the pressure difference 
spectral density, Gapap . The spatial derivative spectrum, Gd x (eq. 5.46) closely resembles the 
pressure difference spectrum in shape for all chordwise positions. This indicates the acoustic 
wavelength is so large that the pressure field appears to be the same over very small chordwise 
separations. The time derivative, Gdt (eq. 5.45) rolls off with nearly the same slope as Gapap for 
op > 3 at all chordwise locations. Gdt begins to fall off for co r < 3 which is expected since the 
time derivative of a slowly varying signal will tend to zero. The convective derivative spectrum, 
Got (eq. 5.42) bears close resemblance to the time derivative spectrum at all chordwise locations 
for op > 3. The rise in G ( / x tends to offset the fall in Gj t resulting in leveling off of Gut at co r < 3 . 
Beyond 9% chord the accuracy of Gdx becomes questionable due to the large chordwise 
separations. In addition, boundary layer pressure fluctuations begin to contaminate the data for 
co, > 20 at 9% chord and beyond. 
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5.11 Summary 

Unsteady pressure measurements made on the surface of a NACA 0015 immersed in 
large grid turbulence (A/c = 13%) are presented and the effects of angle of attack examined. 
Pressure measurements are used to calculate unsteady lift. With increasing angle of attack, lift 
spectra show a decrease by as much as 5dB in spectral level for co, < 10 and an increase of as 
much as lOdB in spectral level for co r > 10. 

The increase in unsteady lift observed for co r > 10 is found to be related to the mean lift of 
the airfoil based on lift spectrum scaling using the mean lift squared. Based on analysis of the 
spanwise pressure correlation length scale, the reduction in unsteady lift is believed to result 
from distortion of the inflow by the mean velocity field. However, predictions of unsteady lift 
made with Reba and Kerschen’s (1996) unsteady loading theory, which accounts for effects of 
mean loading through Rapid Distortion Theory, show a significant rise in lift spectra level with 
increasing angle of attack. Their theory does not capture the correct angle of attack effect at low 
co r because the distortion of the inflow is not accurately predicted with the flat plate model they 
assume. 

As such, a model is developed that accurately predicts the effects of angle of attack on 
unsteady lift spectra. This model uses a circular cylinder fit to the airfoil to compute effects of 
distortion on the inflow turbulence. The distorted inflow velocity spectrum is then used with 
Amiet’s (1976a, b) theory to predict the unsteady loading. This model successfully captures the 
reduction observed in unsteady lift spectra calculated from pressure measurements. 
Furthermore, it is shown that angle of attack effects arising from inflow distortion are significant 
only when the relative scale of the inflow turbulence to airfoil chord is sufficiently small (A/c = 
13% for present experiment). Factors which are well predicted by theories such as Graham 
(1970a), Atassi (1984), Reba and Kerschen (1996) become dominant when the scale of the 
inflow is large compared to the airfoil chord. 

Lastly, it is shown that unsteady pressure measurements can likely be used to make far- 
field noise predictions using the method presented by Casper and Farassat (2002). 
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Figure 5.1: Comparison of measured and predicted pressure difference 
spectra for a = 0° at x/c= 1, 2.5, 4, 6, 9, and 14%. Prediction is Amiet’s 
(1976a, b) unsteady thin airfoil theory. 
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Figure 5.1 continued: Comparison of measured and predicted pressure 
difference spectra for a = 0° at x/c= 1, 2.5, 4, 6, 9, and 14%. Prediction is 
Amiet’s (1976a, b) unsteady thin airfoil theory. 
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Figure 5.1 continued: Comparison of measured and predicted pressure 
difference spectra for a = 0° at x/c= 1, 2.5, 4, 6, 9, and 14%. Prediction is 
Amiet’s (1976a, b) unsteady thin airfoil theory. 
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Figure 5.2: Predictions of pressure difference spectra from the Rapid 
Distortion Theory of Reba and Kerschen (1996) plotted with measured data 
for a = 0°, 2°, and 6° at chordwise locations, x/c = 1, 2.5, 4, 6, 9, and 14% 
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Figure 5.2 continued: Predictions of pressure difference spectra from the Rapid 
Distortion Theory of Reba and Kerschen (1996) plotted with measured data for 
a = 0°, 2°, and 6° at chordwise locations, x/c = 1, 2.5, 4, 6, 9, and 14% 
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Figure 5.2 continued: Predictions of pressure difference spectra from the Rapid 
Distortion Theory of Reba and Kerschen (1996) plotted with measured data for 
a = 0°, 2°, and 6° at chordwise locations, x/c = 1, 2.5, 4, 6, 9, and 14% 
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Figure 5.3: Pressure and suction side root mean squared of/? as a 
function of x/c at a = 0°, formed from integrating pressure spectrum over 
co r = 0.4 - 100. Comparison made with Amiet (1976a, b). 
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Figure 5.4: Pressure and suction side root mean squared of p as a function 
of angle of attack for x/c = 0.01, 0.025, 0.04, 0.06, 0.09, and 0.14. 



Figure 5.5: Large grid root mean square p computed from pressure difference spectra 
as a function of angle of attack for x/c = 0.01, 0.025, 0.04, 0.06, 0.09, and 0.14 
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Figure 5.6: Microphones used in unsteady lift calculation indicated by triangle symbol 
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Figure 5.8: G APAP (oj r , x ,x ’) = G APAP (3, x ,x ’) showing the variation 
of the pressure difference cross-spectra in x 



Figure 5.9: Depiction illustrating how Amiet’s Thin Airfoil Theory is 
used to assess errors related to discretizing lift spectra equations 
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Figure 5.11: G APAP at x/c = 2.5%, co r = 10 from Amiet’s 
theory computed continuously and discretely, plotted 
against all other chordwise locations 
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Figure 5.12b: Depiction of coordinate transformation used to form unsteady lift 
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Figure 5.13: Projection of microphone locations into flow aligned coordinate system for Amiet 

thin airfoil theory simulation 
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Figure 5.14: Unsteady lift computed from Amiet simulation shown for each 

angle of attack 
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Figure 5.15: Drawing showing phase shift in pressure when 
airfoil is at non-zero angle of attack 



Figure 5.17: Ratio of unsteady lift at non-zero angle of attack to zero angle of attack 

unsteady lift based on phasing model 
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Figure 5.18: Unsteady lift calculate from measured unsteady pressure in large grid 

turbulence 
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Figure 5.19: Large grid unsteady lift calculated using pressure measurements from 

the first 14% and the first 85% chord 
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Figure 5.20: Unsteady lift calculated from high angle of attack resolution data for large 
grid turbulence at a = 0° compared with that predicted from Amiet’s (1976a, b) theory. 
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Figure 5.21: Unsteady lift calculated from pressure measurements compared with lift 
computed from Reba and Kerschens (1996) unsteady response theory. 


185 


n n 


10000L" i /(qC) 


Large Grid 



Figure 5.22: Measured L ' 2 as a function of angle of attack compared 
with a curve representing a 2 
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Figure 5.23: Upwash velocity spectrum for large grid turbulence 



Figure 5.24: Admittance function calculated from unsteady pressure measurements 
compared with McKeough’s (1976) measurements 
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Figure 5.25: Comparison of measured and predicted spanwise correlation length scales, a = 0° 


188 


x/c =0.01 



x/c =0.025 



Figure 5.26: Effects of angle of attack on spanwise correlation length scales at x/c = 1 and 2.5% chord 
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Figure 5.27a: Top down view of airfoil at a = 0° showing the reduction in cross-sectional area 
of an eddy encountering an airfoil as a result of eddy stretching by the mean potential field. 
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Figure 5.27b: Top down view of airfoil at a > 0° showing the significant reduction in cross-sectional 
area of an eddy encountering an airfoil as a result of eddy stretching by the mean potential field. 


190 




L c /C L / c L /c 


x/c =0.01 x/c =0.025 



x/c =0.04 x/c =0.06 



x/c =0.09 x/c =0.14 



Figure 5.28: Comparison of measured and predicted chordwise correlation length scales, a = 0° 
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Figure 5.29: Effects of angle of attack on chordwise correlation length scales at x/c = 1 and 2.5% chord 
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Figure 5.30: Contour plots of pressure-velocity correlation function, R pv (pc= 14%, x r) 
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Figure 5.31: Contour plots of R (x=l% 9 x\ f) with trailing edge modification installed 
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Figure 5.32: Unsteady lift calculate with and without the modified trailing edge 
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Figure 5.33: Diagram showing a circular cylinder fit to the airfoils radius of curvature at 

the a = 4° stagnation point location 
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Figure 5.34: Diagram showing the streamline along which the turbulence distortion is calculated 
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Figure 5.35: Unsteady lift spectra calculated using the distorted inflow turbulence energy spectra 
with Amiet’s (1976a, b) solution. Distorted turbulence energy spectra are computed using the 
distortion tensor calculated along the stagnation streamline of a circular cylinder. 
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Figure 5.36: Unsteady lift spectra calculated using the distorted inflow turbulence energy spectra 
with Amiet’s (1976a, b) solution for McKeough’s (1976) flow conditions. Distorted turbulence 
energy spectra are computed using the distortion tensor calculated along the stagnation streamline 

of a circular cylinder. 
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Figure 5.37: Lift spectra scaled on the mean lift squared 
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Figure 5.38: Convective derivative of pressure difference cross-spectral 
density for a = 0° at chordwise locations, x/c = 1, 2.5, 4, 6, 9, and 14% 
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Figure 5.38 continued: Convective derivative of pressure difference 
cross-spectral density for a = 0° at chordwise locations, x/c = 1, 2.5, 4, 6, 

9, and 14% 
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6. Conclusions 

6.1 Primary Contributions 

Measurements of fluctuating surface pressure were made on a NACA 0015 airfoil 
immersed in grid generated turbulence. The airfoil model has a 2' chord and spans the 6' 
Virginia Tech Stability Wind Tunnel test section. Two grids were used to investigate the effects 
of turbulence length scale on the surface pressure response. A large grid which produced 
turbulence with an integral scale 13% of the chord and a smaller grid which produced turbulence 
with an integral scale 1.3% of the chord. Measurements were performed at angles of attack, a 
from 0° to 20°. An array of microphones mounted subsurface was used to measure the unsteady 
surface pressure. The goal of this measurement was to characterize the effects of angle of attack 
on the surface pressure response. Three major contributions have resulted from this work. 

The first significant result is shown in large grid lift spectra calculated from unsteady 
surface pressure measurements. Lift spectra were calculated for each angle of attack and 
comparisons made which revealed two distinct interaction regions; for co r < 10 a reduction in 
unsteady lift of up to 7 decibels (dB) occurs while the unsteady lift increases for co r > 10 as the 
angle of attack is increased. The reduction in unsteady lift at low co r with increasing angle of 
attack is a result that has never before been shown either experimentally or theoretically and as 
such, is a significant and new finding. 

The source of the reduction in lift spectral level appears to be closely related to the 
distortion of inflow turbulence based on analysis of surface pressure spanwise correlation length 
scales. This result suggest eddies are stretched as they approach the airfoil thereby reducing their 
spanwise integral scale which in turn reduces the surface pressure spanwise correlations length 
scale. The effects of distortion are more pronounced as the angle of attack is increased due to the 
larger velocity gradients associated with the mean field which results in a reduction in spanwise 
correlation length scale with increasing angle of attack. The unsteady pressure therefore, has a 
smaller correlated area over which to act and thus, a reduction in unsteady lift occurs with 
increasing angle of attack. Furthermore, while the distortion of the inflow appears to be critical 
in this experiment, wherein the integral scale (A/c) is 13% of the chord, this effect does not seem 
to be significant in larger integral scale (relative to the chord) flows based on the previous 
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experimental work of McKeough (1976). This leads to the second major result of this work; the 
airfoils size relative to the inflow integral scale is critical in defining how the airfoil will respond 
under variation of angle of attack. 

The flat plat theoretical formulations developed by Reba and Kerschen (1996), Graham 
(1970a), and Atassi (1984) appear to accurately predict the airfoil response when the inflow 
integral scale is large relative to the airfoil chord. However, such flat plat predictions schemes 
are not appropriate for the case of smaller scale inflow because they do not model the distortion 
of the inflow correctly. The third major contribution of this work is the development of a 
prediction scheme which correctly accounts for the effects of distortion when the inflow integral 
scale is small relative to the airfoil chord. This scheme utilizes Rapid Distortion Theory to 
account for the distortion of the inflow with the distortion field modeled using a circular 
cylinder. 

6.2 Suggestions for Future Work 

The relative size of the airfoil to inflow integral scale appears to be a critical parameter in 
determining the airfoil response subject to variation of angle of attack. McKeough’s (1976) 
measurements of unsteady lift showed an increase in lift spectral level with increasing angle of 
attack for an inflow integral scale 40% of the airfoil chord. The current experiment shows a 
reduction in unsteady lift when the angle of attack is increased for an inflow integral scale 13% 
of the airfoil chord. This finding suggests the existence of an integral scale to airfoil chord ratio 
which results in an unsteady loading that is independent of angle of attack. Therefore, 
examination of the airfoil response under variation of both angle of attack and A/c will provide 
insight on how the dominant physics vary with A/c. Such an experiment would offer valuable 
details necessary for developing a unified model of the airfoil response. With an instrumented 
airfoil this measurement could be carried out with an active grid for generation of turbulence of 
varying integral scale. 

Analysis of oil flow visualization and unsteady surface pressure data suggests the 
presence of a transient laminar separation bubble in the airfoil leading edge region. It is not clear 
how this bubble, if in fact it really exists, affects the unsteady loading of the airfoil. This bubble, 
if present, may be coupled to the inflow and as such, shedding vorticity which provides means 
for relaxing the unsteady loading on the airfoil. Careful examination of the flow field 
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surrounding the airfoil through both velocity measurements and flow visualization together with 
surface pressure measurements could offer the evidence required to determine this bubbles 
presence and its effect on the airfoil response. 
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Appendix A 

A.l Formulation of Pressure Cross-spectral Density from Measured Pressure 

The first step in producing the surface pressure cross-spectral density from measured 
pressure data is the application of microphone calibrations. This begins by applying the Fourier 
transform to each microphone voltage time series, v m such that 

oo 

V m ( X > y>f) = J V m {x,y,ty i2nft dt (A.l) 


which is numerically implemented as 


N S — (7-l)(*-l) 

v m {x,yJ k ) = Y J v m (x,y,t i y> N ° (A.2) 

j = i 

where j and k are indices for each time and frequency increment, respectively and N s = 2048 (i.e. 
the record length). The microphone calibration, M s [volts/Pa] is applied to each microphone 
voltage spectrum to obtain pressure, p' [Pa] as 


P’(x,y,f) 


V m( X >y>f) 

M s (m n ,fY 


(A.3) 


where m n is a unique number assigned to each microphone. The surface pressure cross-spectral 
density is then formed as 


S m { x n x pVj) = 


N-S„ 


p'(xj,yjj) 


(A.4) 


where E\\ is the expected value, * indicates the complex conjugate, Sr is the sampling rate, and rj 
= / - y. Also, i and j denote different measurement locations. The factor of two in equation 
(A.4) accounts for the fact that (A.l) is only evaluated over the positive time domain. The 
pressure difference cross-spectral density is calculated as 



x »x, 



N-S 


AP(x l ,y„f)AP(x i ,y J ,f)' 


R 


(A.5) 


where A P is the difference between the lower and upper surface pressure measurements 

A P(x,y,f) = P L (x,y,f)-P u (x,y,f). (A.6) 

Throughout this work S qq is used in several different forms which are identified here. 
The surface pressure auto- spectrum, G pp is related to S qq as 
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G pp ( x > /) = S qq ( x i >x t ,i! = 0 , /) (A.7) 

The pressure difference auto-spectrum, Gapap is related to the cross-spectral density as 

G^(x,f) = S‘ (x„x„ n = 0,/) . (A.8) 


A.2 Formulation of Coherence and Pressure Phasing 

The coherence of the unsteady pressure signal across the airfoil is chord is found as 


r 


2 

P 


G p p 

r U r L 

G p u p u G p l p l 


(A.9) 


where G P Pl is the cross-spectrum between the upper and lower surface of the airfoil, G p p is the 

auto-spectrum on the upper surface, and G P f , is the auto-spectrum on the lower surface. The 

corresponding phasing between pressure signals on the upper and lower surface of the airfoil is 
found as 


(j) = arctan 


W 

y Re ( G P L :P L ) j 


(A. 10) 
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